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1.0  SUMMARY 


1.1  Contract  Description 

This  contract  focused  on  the  physics  and  applications  of  coherently  coupled  phase 
conjugators.  Of  particular  interest  is  the  physics  of  mutually  pumped  phase  conjugation  (MPPC) 
and  its  application  to  wavefront  matched  heterodyne  communications.  These  studies  involve 
theoretical  investigations  of  Brillouin  induced  MPPC  in  both  transmission  and  reflection 
geometries,  and  experimental  investigations  of  high  gain  nondegenerate  two-wave  mixing  in 
Cr:YA103. 

1.2  Scientific  Problem 

One  limitation  for  MPPC  is  the  limited  dynamic  range  over  which  photorefractive  materials 
will  operate  in  MPPC  mode.  This  is  a  problem,  for  example,  in  a  wavefront  matched  heterodyne 
receiver  application,  where  the  power  in  an  aberated  incoming  signal  wave  may  be  very  weak  as 
compared  with  the  power  required  for  the  local  oscillator  plane  wave.  To  alleviate  this  problem  we 
studied  alternatives  such  as  using  a  Kerr  media  as  a  replacement  for  photorefractive  media  in 
MPPC. 


It  would  be  desirable  to  do  MPPC  experiments,  and  many  other  nonlinear  optical 
experiments,  using  a  convenient  solid  state  Kerr  media  with  low  power  cw  lasers.  CW  lasers  are 
often  easier  to  work  with  than  high  power  pulse  lasers,  and  they  often  provide  better  quantitative 
data  due  to  the  lack  of  any  need  to  characterize  pulses  and  their  pulse-to-pulse  fluctuations.  When 
this  is  possible,  it  often  brings  about  rapid  progress  in  experimental  research,  as  was  the  case  with 
photorefractives.  Unfortunately,  to  our  knowledge,  no  one  has  demonstrated  the  required  gain  of 
e4  ~  55  in  a  solid  state  bulk  (as  opposed  to  guided  wave)  Kerr  material  using  a  cw  source.  In  this 
program  we  demonstrated  the  largest  cw  nonlinear  optical  gain  in  a  bulk  solid-state 
nonphotorefractive  material.  Although  the  net  gain  was  less  than  that  required  for  MPPC, 
theoretical  calculations  indicate  that  larger  gain  is  possible. 

The  same  optical  nonlinearity  that  gives  rise  to  nonlinear  optical  gain  in  Cr:YAlC>3  also 
gives  rise  to  self-focusing.  To  obtain  an  accurate  description  of  nonlinear  optical  gain  in  Cr:YA103, 
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it  is  necessary  to  include  the  effects  of  self-focusing  in  the  theoretical  model.  One  theoretical  effort 

was  focused  on  an  analysis  of  the  effects  of  self-focusing  on  nonlinear  optical  gain  in  Cr:YA103. 

1.3  Progress  Summary 

There  were  many  areas  of  significant  progress  achieved  under  this  contract  in  the 

development  of  a  solid  state  Kerr  media  for  MPPC  and  in  the  theory  of  MPPC.  These  include: 

•  Completion  of  a  study  of  metal-ion-doped  insulators  as  nonlinear  optical  materials,  and 
identification  of  the  potential  of  Cr:  YAIO3  to  achieve  high  nonlinear  optical  gain  using  cw 
lasers. 

•  Conclusively  demonstrated  the  potential  for  large  dynamic  range  in  Brillouin  induced 
mutually  pumped  phase  conjugation  in  reflection  geometry  by  obtaining  the  first  solutions 
to  the  transcendental  equation  giving  the  phase-conjugate  reflectivity. 

•  First  accurate  determination  of  the  nonlinear  refractive  index  and  saturation  intensity  of 
Cr:YA103. 

•  Procurement,  processing  and  fabrication  of  several  cubes  from  a  highly  doped  boule  of 
Cr:YA103  to  be  used  in  experiments  for  obtaining  large  nonlinear  optical  gain. 

•  A  gain  of  22  times  was  obtained  by  nondegenerate  two-wave  mixing  in  Cr:YA103.  This  is 
currently  the  largest  cw  two-wave  mixing  gain  obtained  in  a  bulk  solid  state 
nonphotorefractive  material,  and  within  a  factor  of  two  of  that  required  for  MPPC. 

•  Observation  of  a  nonlinear  optical  beam  break-up  phenomena  that  is  the  result  of  the 
interaction  between  spatial  inhomogeneities  and  the  optical  nonlinearity  in  a  crystal  of 
Cr:YA103. 

•  Theoretical  prediction  that  significantly  higher  gain,  in  excess  of  200,  is  possible  in  a 
homogeneous  crystal  of  Cr:YA103. 

•  Conclusive  demonstration  of  the  near  equal  strengths  of  transmission  and  reflection 
gratings  in  metal-ion-doped  insulators  like  Cr:YA103,  for  which  there  currently  exists  a 
disagreement  in  the  literature.  This  indicates  that  there  is  very  little  spatial  migration  of  the 
excitation  on  time  scales  corresponding  to  the  excited  state  lifetime  of  33  ms. 
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Analysis  of  transverse  effects  on  nonlinear  optical  gain  in  Kerr  media  indicating  that  a 
significant  enhancement  in  nondegenerate  two-wave  mixing  gain  is  possible  due  to  self- 
focusing  phenomena. 


Details  of  this  progress  are  presented  in  Section  2.0,  and  in  the  publications  listed  below 

that  are  included  in  this  report  as  Appendices. 

1.4  Publications 

1.  I.  McMichael,  T.  Chang,  M.  Noginov,  M.  Curley,  P.  Venkateswarlu  and  H.  Tuller, 
“Application  of  Laser  and  Related  Materials  to  Demonstrate  Large  Nonlinear  Optical  Effects 
and  Diffraction  Efficiency,”  to  be  published  by  the  Optical  Society  of  America  in  the  Trends 
in  Optics  and  Photonics  Series  volume  on  Advanced  Solid  State  Lasers,  1996. 

2.  I.  McMichael  and  T.  Chang,  “Self-Oscillation  and  Self-Pumped  Phase  Conjugation  in 
Cr:YAL03,”  paper  CFF7  to  be  presented  at  the  Conference  on  Lasers  and  Electro-Optics, 
Anaheim,  CA,  June  2-7,  1996. 

3.  I.  McMichael,  T.  Chang,  M.  Noginov  and  H.  Tuller,  “Application  of  Laser  and  Related 
Materials  to  Demonstrate  Large  Nonlinear  Optical  Effects  and  Diffraction  Efficiency,” 
presented  at  the  conference  on  Advanced  Solid  State  Lasers,  San  Francisco,  CA,  January 
31 -February  3,  1996. 

4.  I.  McMichael,  M.  Ewbank  and  F.  Vachss,  “Efficiency  of  Phase  Conjugation  for  Highly 
Scattered  Light,”  Optics  Comm.  1 19, 13  (1995). 

5.  I.  McMichael,  T.  Chang  and  H.  Tuller,  “Optical  Nonlinearity  of  Cr:YA103  and  Other 
Metal-Ion  Doped  Insulators,”  presented  at  the  Annual  Meeting  of  the  Optical  Society  of 
America,  Portland,  OR,  September  10-15,  1995. 

6.  I.  McMichael,  R.  Saxena,  T.  Chang,  Q.  Shu,  S.  Rand,  J.  Chen  and  H.  Tuller,  “High  Gain 
Nondegenerate  Two-Wave  Mixing  in  Cr:YA103,”  Opt.  Lett.  19,  1511  (1994). 

7.  I.  McMichael,  M.  Ewbank  and  F.  Vachss,  “Efficiency  of  Phase  Conjugation  for  Highly 
Scattered  Light,”  presented  at  the  Annual  Meeting  of  the  Optical  Society  of  America, 
Dallas,  TX,  October  2  -  7,  1994. 

8.  R.  Saxena,  T.  Chang  and  I.  McMichael,  “Effects  of  Self-Focusing  on  Nondegenerate 
Two-Wave  Mixing  in  Kerr  Media,”  presented  at  the  Annual  Meeting  of  the  Optical  Society 
of  America,  Dallas,  TX  (1994). 

9.  R.  Saxena  and  I.  McMichael,  “Brillouin  Induced  Mutually  Pumped  Phase  Conjugation  in 
Reflection  Geometry,”  presented  at  the  Conference  on  Nonlinear  Optics,  Waikoloa,  Hawaii 
(1994). 
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10.  I.  McMichael,  R.  Saxena,  T.  Chang,  Q.  Shu,  S.  Rand,  J.  Chen  and  H.  Tuller,  “High  Gain 
Nondegenerate  Two-Wave  Mixing  in  Cr:YA103,”  presented  at  the  Conference  on 
Nonlinear  Optics,  Waikoloa,  Hawaii  (1994). 

11.  I.  McMichael,  R.  Saxena,  T.  Chang,  R.  Neurgoankar,  and  S.  Rand,  High  Gain 
Nondegenerate  Two-Wave  Mixing  in  Cr:YA103,”  in  Conference  on  Lasers  and  Electro- 
Optics,  Vol.  11  of  1993  OSA  Technical  Digest  Series  (Optical  Society  of  America, 
Washington,  D.C.,  1993),  paper  CThS34,  pp.  514-515. 

12.  R.  Saxena  and  I.  McMichael,  “Brillouin  Induced  Mutually  Pumped  Phase  Conjugation  in 
Reflection  Geometry,”  submitted  to  the  Journal  of  the  Optical  Society  of  America  B. 

13.  R.  Saxena  and  P.  Yeh,  “Nonlinear  Optical  Phase  Conjugation  by  Four-Wave  Mixing  in 

Kerr  Media,”  in  Phase  Conjugation,  M.  Gower  and  D.  Proch,  eds.,  (Springer- Verlag 
Berlin,  1994),  Chap.  2.  5 

14.  R.  Saxena,  Diffraction  Properties  of  Photorefractive  Gratings,”  in  Photorefractive 
Materials,  Effects  and  Applications;  Critical  Reviews,  P.  Yeh  and  C.  Gu,  eds.,  (Society  of 
Photo-Optical  Instrumentation  Engineers,  Bellingham,  Washington,  1993). 
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2.0  PROGRESS 


2.1  Identification  of  the  Potential  of  Cr:YAL03 

As  noted  in  the  summary,  a  convenient  solid  state  Kerr  media  with  large  cw  nonlinear 
optical  gain  is  required  in  experiments  on  MPPC  for  applications  such  as  wavefront  matched 
heterodyne  communications.  We  had  previously  demonstrated  that  nonlinear  optical  gain  could  be 
achieved  in  ruby  by  the  process  of  nondegenerate  two-wave  mixing  of  two  beams  from  a  cw  argon 
laser.1  However,  calculations  indicated  that  it  would  be  difficult  to  achieve  the  gain  required  for 
MPPC  using  ruby.  We  therefore  initiated  a  search  of  the  literature  to  see  if  some  other  metal-ion- 
doped  insulator  might  work,  and  found  that,  of  the  materials  previously  investigated,  Cr:YAlC>3 
was  the  most  promising  for  achieving  the  high  gain  required  in  MPPC.2 

2.2  First  Demonstration  of  High  Nonlinear  Optical  Gain  in  Cr:YAI03 

After  identifying  Cr:YA103  as  a  promising  candidate  material  for  future  investigations  of 
MPPC  in  Kerr  media,  we  then  focused  our  attention  on  demonstrating  the  potential  for  obtaining 
high  nonlinear  optical  gain.  Using  a  borrowed  sample  of  Cr:YA103,  we  demonstrated  a  beam 
coupling  gain  of  6  times.3  At  that  time,  this  was  the  largest  cw  nonlinear  optical  gain  obtained  in  a 
bulk  solid  state  nonphotorefractive  material.  Details  of  this  experiment  are  presented  in  the  paper 
entitled  “High  Gain  Nondegenerate  Two-Wave  Mixing  in  Cr:YAlC>3,”  that  was  presented  at  the 
Conference  on  Lasers  and  Electro-Optics  and  is  included  in  the  Appendix.  The  important  data  are 
shown  in  Fig.  1,  where  gain  for  a  weak  probe  beam  is  plotted  as  a  function  of  the  pump  intensity. 
A  fit  to  this  data  yields  the  nonlinear  refractive  index  112=  1.4  xlO'7  cm2/W  and  the  saturation 
intensity  Is  =  1.5  kW/cm2.  For  comparison,  we  find  that  112/a  =  3  x  10'7  cm2/W  for  Cr:YA103,  is 
more  than  one  order  of  magnitude  larger  than  the  value  for  ruby,  ri2/(X  =  1.7  x  10' 8  cm2/W.  It  is 
also  satisfying  that  the  value  of  the  saturation  intensity  is  consistent  with  calculations  using  the 
expression  Is  =  hv/GT  together  with  values  of  the  cross  section  a  and  the  decay  time  x  reported  in 
the  literature. 
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Fig.  1  Nonlinear  optical  gain  in  Cr:YA103  as  a  function  of  the  pump  intensity. 

2.3  Procurement,  Processing  and  Fabrication  of  Highly  Doped  CriYAIQa 

We  recognized  that  it  should  be  possible  to  obtain  higher  gain  than  reported  above  by  using 
higher  Cr  doping  and  a  longer  interaction  length.  We  procured  a  crystal  boule  with  which  we  later 
obtained  significantly  higher  gain,  and  allowed  us  to  continue  our  claim  of  demonstrating  the 
largest  cw  nonlinear  optical  coupling  in  a  non-photorefractive  media.  A  survey  of  several  potential 
vendors  and  several  individuals  indicates  that  this  crystal,  sitting  on  a  shelf  to  be  sold  as  “gem” 
material,  may  have  been  the  only  crystal  available  in  the  US  at  that  time. 

As  grown  Cr:YA103  has  a  color  center  absorption  that  competes  with  the  absorption 
responsible  for  the  optical  nonlinearity.  A  collaboration  was  established  with  Harry  Tuller  at  MIT 
who  removed  the  color  centers  by  heating  the  as  grown  material  in  a  forming  gas  atmosphere.4 
Figure  2  shows  absorption  spectra  for  the  12  mm  thick  end  of  the  boule  of  Cr:YAlC>3,  as  grown 
(dotted  line),  after  a  first  anneal  (dashed  line)  in  95%  N2/5%  H2  at  1050°  C  for  8  hours  with  a  flow 
rate  of  250  cc/min  and  a  ramp  rate  of  1  °C/min,  and  after  a  second  anneal  (solid  line)  at  1300°  C 
for  24  hours  with  the  same  atmosphere,  flow  and  ramp  rates  as  used  in  the  first  anneal.  After  the 
second  anneal,  the  boule  was  fabricated  into  several  crystallographic  cut  cubes,  one  of  which  was 
used  in  the  following  experiments.  Microprobe  measurements  indicated  that  the  Cr  concentration 
was  0. 12  weight  percent. 
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Fig.  2  Cr:YA103  absorption  spectra  showing  the  removal  of  color  centers  by  processing. 

Following  the  heat  treatment  a  strategy  was  developed  to  obtain  one  large  cube  and  several 
smaller  cubes  from  the  boule,  specifications  were  established,  and  the  boule  was  sent  out  for  x-ray 
orientation  and  fabrication.  Estimates  of  the  potential  gain  in  the  larger  cube  are  of  the  order  of  400 
times,  when  one  ignores  transverse  effects  such  as  self-focusing.  The  total  cost  of  materials, 
processing  and  fabrication  of  the  cubes  will  be  less  than  the  cost  of  a  single  BaTi03  crystal. 

2.4  Demonstration  of  High  INonlinear  Optical  Gain  in  Highly  Doped 
Cr:YA103 

Based  on  the  measured  value  of  n2/oc,  calculations  indicated  that  one  could  obtain 
significantly  higher  gain  by  increasing  the  doping  and  interaction  length.  A  crystal  with  higher 
doping  was  procured,  processed  and  fabricated,  and  then  used  in  experiments  to  demonstrate 
higher  gain.  Details  of  these  experiment  are  presented  in  the  papers  entitled,  “High  Gain 
Nondegenerate  Two-Wave  Mixing  in  Cr:YA103,”  that  were  published  in  Optics  Letters5  and 
presented  at  the  Conference  on  Nonlinear  Optics,6  and  are  included  in  the  Appendix.  The  important 
data  are  shown  in  Fig.  3,  where  gain  for  a  weak  probe  beam  is  plotted  as  a  function  of  the  pump 
intensity.  Figure  3  (a)  shows  a  plot  of  the  gain  as  a  function  of  the  intensity  of  the  pump  for  the 
range  from  0  to  80  W/cm2.  The  circles  are  the  experimental  measurements,  and  the  line  is  a  fit  of 
the  theoretical  form, 
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with  a L  -  0.69,  L  =  9. 1  mm,  X  =  5 14.5  nm  and  Is  =  hv/ax  =  1 .27  kW/cm2,  where  h  is  Planck’s 
constant,  v  is  the  optical  frequency,  a  is  the  absorption  cross-section  and  x  is  the  metastable  state 
lifetime.  Using  n2  as  the  only  adjustable  parameter,  the  fit  yields  n2  -  3.2  x  10-7  cm2/W.  For 
comparison,  we  find  that  n2/a  =  4.2  x  1(F7  cm2/W  for  Cr:YA103,  is  more  than  one  order  of 
magnitude  larger  than  the  value  for  ruby,  n2/a  =  1.7  x  1(F8  cm2/W. 


(a)  (b) 

Fig.  3  Nonlinear  optical  gain  in  Cr:YAlC>3  as  a  function  of  the  pump  intensity,  (a)  from  0 
to  70  W/cm2,  and  (b)  from  0  to  400  W/cm2. 


Figure  3  (b)  shows  the  gain  for  the  range  from  0  to  400  W/cm2.  Using  the  previously 
stated  parameters  and  n2  as  determined  from  the  fit  for  the  lower  intensity  range  of  Fig.  3(a),  the 
line  in  Fig.  3(b)  represents  projections  based  on  the  theoretical  form  expressed  by  Eq.  (1).  A  gain 
of  22  times  was  obtained  at  approximately  400  W/cm2,  but  the  measured  gain  does  not  increase 
with  intensity  as  much  as  theory  predicts.  (Theory  predicts  a  gain  of  200  at  400  W/cm2.)  As 
described  in  the  following  section,  at  first  we  believed  this  was  due  to  “beam  break-up”  resulting 
from  spatial  nonuniformities  of  the  nonlinear  refractive  index  change  in  the  interaction  region. 
However,  as  we  discuss  in  later  in  Section  2.6,  measurements  of  the  luminescence  as  a  function  of 
pumping  intensity  indicate  that  the  reduction  in  beam  coupling  gain  is  due  mostly  to  a  reduction  of 
the  metastable  state  population,  which  is  greater  than  that  predicted  by  depletion  of  the  ground 
state,  as  expressed  in  Eq.  (1),  when  using  the  usual  expression  for  the  saturation  intensity,  I§  - 
hv/ax. 
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2.5  Nonlinear  Optical  Beam  Breakup 

In  our  experiments  on  high  gain  nondegenerate  two-wave  mixing  in  Cr:YA103,  we 
observed  a  nonlinear  optical  beam  break-up  phenomena  that  is  the  result  of  the  interaction  between 
spatial  inhomogeneities  and  the  optical  nonlinearity  in  the  crystal.  This  phenomena  is  illustrated  in 
Fig.  4.  The  figure  shows  a  transmitted  pump  beam,  (a)  at  low  intensity  and  normal  exposure, 
indicating  the  transmitted  beam  size,  (b)  at  low  intensity  and  long  exposure,  showing  an  “X” 
shaped  pattern  of  linear  scattering  that  is  believed  to  be  the  result  of  scattering  by  crystal  striae  and 
(c)  at  high  intensity  and  exposure  equivalent  to  case  (b),  showing  the  nonlinear  scattering 
phenomena  due  to  nonuniformities  in  the  nonlinear  refractive  index  change. 

Spatial  nonuniformities  of  the  nonlinear  refractive  index  change  in  the  interaction  region  can 
be  the  result  of  shadows  or  distortions  from  surface  or  bulk  defects,  or  from  spatial  variations  in 
the  Cr-ion  concentration.  Photographs  of  the  crystal  illuminated  by  a  low  intensity  collimated  beam 
show  radial  and  concentric  patterns  are  thought  to  be  Cr  striae.  Images  of  a  low  intensity  pump 
beam  inside  the  crystal  also  indicate  that  the  striae  have  a  significant  effect  on  the  beam  profile. 
Beam  break-up  can  result  from  both  variations  in  Cr  concentration,  yielding  a  variation  in  n2,  and 
also  from  intensity  variations,  yielding  variations  in  the  induced  nonlinear  index  change. 


(a)  ■'  (b)  (c) 


Fig.  4  Photographs  of  a  beam  transmitted  through  Cr:YA103  illustrating  the  nonlinear  optical 
beam  breakup  phenomena,  (a)  Gaussian  beam  size,  (b)  “X”  shaped  pattern  resulting 
from  linear  scattering,  and  (c)  nonlinear  scattering. 
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Qualitative  behavior  indicated  that  the  scattering  at  high  intensity  was  nonlinear,  but  to 
confirm  this  we  measured  the  intensity  of  a  small  region  of  the  scattering  as  a  function  of  the  pump 
intensity.  An  aperture  was  closed  down  to  a  diameter  such  that,  at  low  intensity,  it  would  pass  one 
half  of  the  pump  that  transmitted  through  the  crystal,  and  then  it  was  placed  in  the  “speckle”  of  the 
scatter  indicated  by  the  arrow  in  Fig.  2  (c).  Figure  5  shows  the  scattered  power  divided  by  the 
pump  power  as  a  function  of  the  pump  power.  The  circles  are  the  data,  and  the  line  is  a  quadratic 
fit,  indicating  that  the  nonlinear  portion  of  the  scatter  has  a  cubic  dependence  on  the  pump 
intensity. 


Pump  Intensity  (kW/cm2) 


Fig.  5  Measurements  showing  the  nonlinear  scattering  in  Cr:YA103. 

We  made  the  following  observations  concerning  the  scattering:  first,  the  strength  of  the 
nonlinear  scattering  depended  on  feedback  into  the  crystal,  and  second,  a  threshold  existed  above 
which  the  speckle  pattern  of  the  scattering  was  dynamic.  An  experiment  looking  at  the  interference 
between  a  portion  of  the  transmitted  pump  and  a  “speckle”  did  not  detect  any  frequency  shift.  This 
indicated  that  the  dynamic  scattering  at  high  intensities  was  not  due  to  beating  between  frequency- 
shifted  stimulated  scattering  resulting  from  NDTWM  gain,  and  some  linear  scattering  from  crystal 
inhomogeneities.  When  the  crystal  was  oriented  to  retroreflect  the  pump,  or  when  an  external 
mirror  was  used,  the  strength  of  the  nonlinear  scattering  increased  and  the  threshold  of  dynamic 
scattering  decreased.  To  obtain  the  highest  gains  in  the  NDTWM  measurements,  we  rotated  the 
crystal  away  from  normal  to  reduce  the  nonlinear  scattering.  An  angle  of  incidence  for  the  pump  of 
0.03  radians  yielded  the  largest  gains.  Further  reductions  in  the  feedback  resulting  from  the  rear 
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surface  reflection  might  improve  the  gain.  This  could  be  achieved  by  an  antireflection  coating  01  by 
index  matching,  but  increasing  the  angle  further  may  not  help  if  it  reduces  phase  matching. 

2.6  Luminescence  Measurements 

Figure  6  shows  measurements  of  the  luminescence  intensity  from  the  highly  doped  crystal 
ofCr:YA103  as  a  function  of  the  pump  beam  intensity.  The  crosses  are  the  measurements,  and  the 
solid  line  is  a  fit  of 


This  functional  form  describes  the  saturation  of  the  luminescence  due  to  depletion  of  the  ground 
state  for  pumping  by  a  beam  with  a  Gaussian  spatial  profile.  The  fit  gives  ls  =  870  W/cm2,  which 
is  less  than  the  value  of  1,270  W/cm2  calculated  using  Is  =  hv/crc.  The  nonlinearity  of  the 
luminescence  vs.  pumping  can  be  due  to  a  number  of  factors  including  not  only  depopulation  of 
the  ground  state,  but  also  up-conversion,  excited  state  absorption  (ESA),  and  possibly  other 
effects.  However,  for  the  purpose  of  subsequent  calculations  these  effects  aie  conveniently 
expressed  in  terms  of  an  effective  saturation  intensity  that  results  from  the  fit  of  the  above  equation. 
The  subsequent  calculations  using  this  effective  saturation  intensity  should  be  correct  to  hist  oidei. 


0  200  400  600  800  1000 


Pump  Intensity  1 1  (0)  [W/cm2] 


Fig.  6  Luminescence  vs.  pump  intensity  in  the  highly  doped  crystal  of  Cr.YAlOy. 
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Using  the  value  for  Is  from  the  fit  in  Fig.  6,  we  can  refit  the  beam  coupling  gain  as  shown 
in  Fig.  3.  The  line  is  a  fit  to  the  range  from  0  to  80  W/cm2  that  gives  n2  Id  =  4.0  x  10~7  cm3/W. 
The  new  fit  is  much  better  than  the  previous  one  shown  in  Fig.  3.  Most  of  the  remaining 
discrepancy  at  high  pump  intensity  can  be  explained  by  the  nonexponential  decay  of  the 
luminescence  that  takes  place  at  high  intensity  (see  section  on  spectroscopy  and  kinetics).  The  fit  in 
Fig.  3  assumed  that  the  pump  and  probe  beams  were  plane  waves  of  infinite  extent.  The  fit  in  Fig. 
5  includes  the  effects  of  a  Gaussian  beam  profile.  Although  the  effect  is  relatively  small,  it  does 
yield  a  slightly  better  fit  to  the  data. 


Fig.  7  New  fit  to  the  beam  coupling  measurements  of  Fig.  3,  using  Is  as  determined  from 
a  fit  to  Fig.  6. 

2.7  A  New  Sample  of  Cr:YAI03  Exhibits  Different  Saturation 

We  had  a  new  crystal  of  Cr:YAlC>3  grown  by  Preciosa  Crytur,  Ltd.,  which  we  refer  to  as 
sample  2  (we  refer  to  the  previous  crystal  as  sample  1)  Since  the  Cr  concentration  in  this  crystal 
(0.028  wt.%)  was  less  than  we  specified  (0.1  wt.%),  it  did  not  exhibit  as  high  a  gain  as  the 
previous  crystal  (0.12  wt.%),  but  it  also  did  not  exhibit  the  premature  saturation  behavior.  Figure  8 
shows  measurements  of  the  luminescence  vs.  pump  intensity  for  sample  2.  The  fit  gives  1$  =  1640 
W/cm2,  which  is  greater  than  870  W/cm2  obtained  for  sample  1.  Thus,  sample  2  does  not  show 
premature  saturation,  and  in  fact  the  saturation  appears  to  be  less  than  that  predicted  by 
Is  =  hv/ax.  One  possible  explanation  for  the  measured  value  being  greater  than  the  predicted 
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value  is  that  the  cross  section  may  contain  terms  due  to  residual  color  center  absorption  that  do  not 
contribute  to  the  metastable  state  population. 


Fig.  8  Luminescence  vs.  pump  intensity  in  sample  2. 

Figure  9  shows  measurements  of  the  beam  coupling  gain  in  sample  2.  The  dashed  line  is  a 
fit  using  the  value  of  n?  /a  from  sample  1 ,  and  Ig=l  ,270  W/cm2  It  fits  well  at  low  intensity,  but  it 
greatly  underestimates  the  gain  obtained  at  high  intensity.  The  solid  line  is  a  fit  that  gives  n2/ot  = 
4. 1  x  I  O'7  cnvVW  and  =  1 ,630  W/cm2.  These  values  are  in  excellent  agreement  with  the  value  of 
n2  /cc  from  sample  1,  and  the  value  of  1$  from  the  lit  of  the  luminescence  for  sample  2. 


Fig.  9  Gain  vs.  pump  intensity  in  the  new  sample  of  Cr:YAlC>3. 
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To  conclude  this  section,  by  fitting  the  luminescence  vs.  pump  intensity,  we  obtain  an 
effective  saturation  intensity  that  lumps  together  the  effects  of  ground  state  depletion,  up- 
conversion,  ESA,  etc.  Using  this  effective  saturation  intensity,  we  obtain  a  better  fit  to  the  beam 
coupling  gain  than  is  possible  using  the  expression  Is  =  hv/ax.  Thus,  the  processes  that  give  rise  to 
saturation  of  the  luminescence  and  the  metastable  state  population,  rather  than  the  beam  breakup 
that  we  had  previously  proposed,  are  the  major  factors  responsible  for  the  premature  saturation  of 
the  beam  coupling  gain  in  sample  1. 

2.8  Self-Oscillation  in  Cr:YAI03 

At  normal  incidence,  the  reflectivity  of  Cr:YA103  is  1 1%,  and  the  absorption  loss  in  sample 
2  is  17%.  Thus,  the  round-trip  gain  (122  =  144)  is  sufficient  to  overcome  the  losses  ( (0.1 1x0. 83)2 
=  1/120)  of  a  resonator  formed  by  the  opposing  surfaces  of  the  crystal.  (We  previously 
demonstrated  that  the  gain  for  a  reflection  grating  is  comparable  to  that  for  a  transmission  grating.) 
With  the  probe  off  and  the  crystal  normal  oriented  at  a  small  angle  with  respect  to  the  incident 
pump  beam,  an  oscillation  builds  up  between  the  opposing  surfaces  of  the  crystal.  Figure  10 
shows  far  field  patterns  of  the  oscillating  mode.  As  shown  in  the  figure,  it  is  possible  to  change  the 
oscillating  mode  by  translating  the  crystal  perpendicular  to  the  pumping  beam.  Using  a 
microscope,  we  observe  a  moving  fringe  pattern  in  the  fluorescence  from  the  interaction  region  of 
the  crystal  that  is  formed  by  the  interference  between  the  beam  incident  on  the  crystal  and  the 
frequency  shifted  self-oscillation.  By  beating  the  oscillation  with  a  portion  of  the  incident  beam,  we 
measure  a  frequency  shift  of  6  Hz.  This  is  consistent  both  with  our  measurements  of  the  frequency 
shift  at  which  the  nondegenerate  two-wave  mixing  gain  was  maximum,  and  with  the  theoretical 
expression,  l/27tx,  where  x  is  the  response  time  of  the  nonlinearity  (the  lifetime  of  the  metastable 
state,  33  ms). 
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Fig.  10  Far  field  patterns  of  the  oscillating  mode  from  a  self-pumped  linear  oscillator  that 
is  formed  in  a  crystal  of  Cr:YA103. 

The  bidirectional  oscillation  acts  as  pumping  beams  and  the  beam  incident  on  the  crystal 
acts  as  a  probe  beam  in  a  four-wave  mixing  process  to  generate  what  we  believe  is  the  self-pumped 
phase-conjugate  of  the  beam  incident  on  the  crystal.  Figure  1 1  shows  the  resulting  phase-conjugate 
return.  We  verified  that  the  return  was  along  a  direction  consistent  with  a  retroreflection  of  the 
incident  beam  and  that  it  built  up  with  the  oscillation.  This  is  consistent  with  the  interpretation  of  a 
phase-conjugate  return.  However,  we  have  not  yet  verified  the  phase-conjugate  nature  of  the 
wavefront  by,  for  example,  demonstrating  distortion  correction.  Note  that  the  return  is  elongated  in 
the  direction  perpendicular  to  the  plane  of  incidence  defined  by  the  pump  and  probe  beams.  A 
similar  effect  was  previously  observed  in  photorefractive  ring  conjugators. 
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Fig.  1 1  Phase-conjugate  return  from  the  self-pumped  linear  oscillator  that  is  formed  in  a 
crystal  of  Cr:YA103. 
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A.1 

Application  of  Laser  and  Related  Materials  to  Demonstrate 
Large  Nonlinear  Optical  Effects  and  Diffraction  Efficiency 
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Abstract 

Long  metastable  lifetimes  in  laser  and  related  materials 
combined  with  a  large  index  change  makes  possible 
large  nonlinear  optical  effects  using  cw  lasers.  We 
describe  experiments  demonstrating  such  effects  using 
the  solid  state  laser  material  Cr:YA103,  and  we  show 
how  studies  of  the  spectroscopy  and  kinetics  can  be 
used  to  understand  the  results.  Since  the  same  optical 
nonlinearity  will  affect  the  propagation  of  beams,  it 
should  be  considered  in  the  design  of  lasers. 

Keywords 

nonlinear  optics,  nonlinear  optical  materials,  laser 
materials,  rare  earth  and  transition  metal  solid  state 
lasers 


Introduction 

Laser  and  related  materials  can  be  used  to  demonstrate 
large  nonlinear  optical  effects  using  cw  lasers.  For 
example,  using  Cr:YA103  we  have  demonstrated  an 
optical  gain  of  22  for  a  weak  probe  beam  by  using  a 
moving  grating  technique  [1].  Furthermore,  we  have 
used  this  gain  to  demonstrate  "self-oscillation"  in  which 
an  oscillation  spontaneously  builds  up  between  the 
opposing  surfaces  of  a  crystal,  when  a  single  beam  is 
incident  [2].  Since  these  materials  are  solid  state,  they 
are  much  easier  to  use  than  the  metal-vapors  that  also 
have  large  optical  nonlinearities.  Photorefractives  are 
solid  state  materials  that  can  be  used  to  demonstrate 
large  nonlinear  optical  effects  using  cw  lasers,  but  their 
nonlinear  optical  response  is  typically  nonlocal  and  not 
proportional  to  intensity.  Thus,  laser  and  related 
materials  may  play  an  important  role  in  nonlinear  optics 
by  providing  convenient  solid  state  materials,  with  a 


local  nonlinear  optical  response  proportional  to 
intensity,  that  will  allow  researchers  to  demonstrate 
large  nonlinear  optical  effects  using  cw  lasers. 
Furthermore,  since  the  same  optical  nonlinearity  will 
affect  the  propagation  of  beams,  it  should  be  considered 
in  the  design  of  lasers. 

In  this  paper  we  first  present  results  of  a  simple 
model  of  the  optical  nonlinearity  in  laser  materials. 
These  results  can  be  used  to  find  rules  of  thumb  for 
developing  materials  with  large  refractive  index 
changes.  As  an  example  of  such  a  material,  the  second 
section  describes  our  experiments  demonstrating  large 
nonlinear  optical  gain  using  chromium-doped  yttrium 
aluminate  (Cr:YA103).  We  study  two  crystals  from 
different  growers  that  we  refer  to  as  sample  1  and 
sample  2,  with  Cr  concentrations  of  0.12  and  0.028 
wt.%  respectively.  From  measurements  of  the 
luminescence  vs.  pumping  intensity,  we  obtain  an 
effective  saturation  intensity  that  lumps  together  the 
effects  of  various  mechanisms  causing  saturation  of  the 
metastable  state  population.  Using  this  effective 
saturation  intensity,  we  are  able  to  account  for  most  of 
the  premature  saturation  of  the  beam  coupling  gain  that 
is  observed  in  sample  1 .  The  final  section  of  our  paper 
describes  studies  of  the  spectroscopy  and  kinetics  that 
yield  clues  as  to  the  mechanisms  behind  the  effective 
saturation  intensity.  We  hope  that  this  work  will 
stimulate  others  to  help  us  find  better  materials  and 
greater  understanding. 

Results  of  a  Simple  Model  of  the  Optical 
Nonlinearity  of  Laser  Materials 

The  optical  nonlinearity  in  laser  materials  results  from 
the  light  induced  population  of  a  metastable  state  and 
the  accompanying  change  of  the  refractive  index.  Long 
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metastable  state  lifetimes  make  it  possible  to  achieve  a 
large  population  of  the  metastable  state,  and  thus  a  large 
change  in  refractive  index,  using  relatively  low  power 
cw  lasers.  At  least  as  long  ago  as  1977,  the  dominant 
contribution  to  the  light  induced  change  of  the 
refractive  index  in  ruby  was  explained  [3]  by 
considering  just  the  strong  charge-transfer  (CT) 
transitions  [4-6].  In  ruby,  the  CT  transitions  are 
allowed  transitions  in  which  an  electron  is  transferred 
from  a  nonbonding  orbital,  localized  predominantly  on 
the  O'*  ligands,  to  an  antibonding  orbital  on  the  metal 
ion,  and  it  results  in  a  very  strong  absorption  in  the  UV 
at  approximately  180  nm.  This  suggestion  continues  to 
be  used  in  the  literature  to  explain  the  nonlinear 
refractive  index  of  laser  materials.  If  we  follow  this 
suggestion,  then  we  can  show  that  the  maximum  index 
change  Anmax  is  given  by, 
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where  n2  is  the  nonlinear  refractive  index,  Is  =  hv/ox  is 
the  saturation  intensity,  h  is  Planck's  constant,  v  is  the 
frequency  of  the  light  interacting  with  the  material,  c  is 
the  absorption  cross-section,  x  is  the  metastable  state 
lifetime,  n0  is  the  linear  refractive  index,  N  is  the 
number  density  of  the  active  species  (i.e.,  Cr  for  the 
case  of  ruby),  e  is  the  elementary  charge,  e0  is  the 
permit-tivity  of  vacuum,  m  is  the  electron  mass,  f  j  and 
(Oj  are  the  oscillator  strength  and  angular  (i.e.  co  =  2  v) 
frequency  for  the  transition  from  the  metastable  state  to 
the  CT  state,  respectively;  f0  and  co0  are  the  oscillator 
strength  and  angular  frequency  for  the  transition  from 
the  ground  state  to  the  CT  state,  respectively;  co  is  the 
angular  frequency  of  the  light  interacting  with  the 
material.  This  expression  is  closely  related  to  that  for 
the  polarizability  used  by  Powell  and  Payne  [7]. 

Assuming  we  can  detune  the  light  frequency  to 
reduce  absorption,  and  that  x  is  long  enough  that  we  can 
reach  Is  with  a  cw  laser,  then  this  equation  gives  the 
figure-of-merit  for  obtaining  a  large  refractive  index 
change.  It  indicates  the  desirability  of  a  large  linear 
refractive  index,  a  large  doping,  large  oscillator  strength 
for  the  transition  from  the  metastable  state  to  the  CT 
state,  small  oscillator  strength  for  the  transition  from  the 
ground  state  to  the  CT  state,  a  CT  state  that  is  closer  to 
the  metastable  state,  and  a  metastable  state  energy  that 
is  close  to  the  laser  photon  energy.  However,  when  a 
CT  state  closer  to  the  metastable  state  is  considered,  at 
some  point  the  resulting  absorption  must  be  included. 
Also  note  that  under  these  conditions  the  two  frequency 
dependent  terms  in  this  equation  can  have  the  same 


negative  sign  and  combine  to  yield  an  even  larger 
effect. 

Although  we  have  discussed  only  the  CT  transition, 
since  historically  this  has  been  considered  as  the 
dominant  contribution  in  ruby,  the  above  equation  for 
the  maximum  index  change  has  a  more  general  validity. 
That  is,  the  above  discussion  also  applies  as  an 
approximation  if  the  CT  transition  is  simply  replaced 
with  the  appropriate  dominant  transition  for  a  given 
material,  and  to  a  better  approximation  if  a  sum  is  taken 
over  all  possible  transitions. 


Demonstration  of  Large  Nonlinearity 

Figure  1  is  a  simplified  schematic  of  a  setup  we  have 
used  to  measure  gain  for  a  weak  probe  in  a  moving 
grating  experiment  [1].  This  powerful  technique  can  be 
used  to  determine  both  the  real  and  imaginary  parts  of 
the  nonlinear  refractive  index  and  the  decay  time  for  the 
metastable  state  [8].  Light  from  an  argon  laser  is  split 
into  a  weak  probe  beam  I2(0)  and  a  strong  pump  7/0), 
with  a  fixed  ratio  12(0)11  j(0)  =  1/1000.  The  frequency 
of  the  pump  is  shifted  by  approximately  5  Hz  with 
respect  to  the  probe  by  reflecting  it  from  a  mirror 
mounted  to  a  piezoelectric  transducer  and  driven  by  a 
triangle-wave  voltage  source.  The  grating  resulting 
from  the  interaction  of  the  pump  and  probe  in  the 
crystal  of  Cr:YA103  results  in  amplification  for  the 
probe  when  the  mirror  is  moving  in  one  direction,  and 
attenuation  when  the  mirror  moves  in  the  opposite 
direction. 

For  sample  1,  Fig.  2  shows  the  gain  for  the  probe 
beam  as  a  function  of  the  pump  beam  intensity  for  the 
range  from  0  to  80  W/cm2.  The  crosses  are  the 
measurements,  and  the  line  is  a  fit  of  a  theoretical 
equation  for  the  beam  coupling  gain  [1], 
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with  aL  =  0.69,  L  =  9.1  mm,  X  =  514.5  nm  and  Is  = 
hv/ox  =  1.27  kW/cm2.  Using  n2  as  the  only  adjustable 
parameter,  the  fit  yields  n2/a  =  4.2  x  10'7  cm3/W. 
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Figure  1.  Schematic  of  the  setup  to  measure  gain  for  a  weak 
probe  in  a  moving  grating  experiment. 
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Luminescence  Measurements 


¥ 


Pump  Intensity  1^0)  [kW/cm2] 


Figure  2.  Gain  vs.  pump  intensity  from  0  to  80  W/cm2. 

Figure  3  shows  the  gain  for  the  range  from  0  to  400 
W/cm2.  Using  n2  as  determined  from  the  fit  for  the 
lower  intensity  range  of  Fig.  2,  the  line  in  Fig.  3 
represents  projections  based  on  Eq.  (2).  A  gain  of  22 
was  obtained  at  approximately  400  W/cm2,  but  the 
measured  gain  does  not  increase  with  intensity  as  much 
as  the  simple  theory  predicts  (a  gain  of  200  at  400 
W/cm2).  Previously  [1]  we  suggested  that  this  might  be 
due  to  a  "beam  break-up"  phenomenon  we  observed 
and  attributed  to  spatial  nonuniformities  of  the 
nonlinear  refractive  index  change  in  the  interaction 
region.  However,  as  we  discuss  in  the  following 
section,  measurements  of  the  luminescence  as  a 
function  of  pumping  intensity  indicate  that  the 
reduction  in  beam  coupling  gain  is  due  mostly  to  a 
reduction  of  the  metastable  state  population,  which  is 
greater  than  that  predicted  by  depletion  of  the  ground 
state,  as  expressed  in  Eq.  (1),  when  using  the  usual 
expression  for  the  saturation  intensity,  Is  =  hv/crc. 


Fig.  3.  Gain  vs.  pump  intensity  from  0  to  400  W/cm2. 


Figure  4  shows  measurements  of  the  luminescence 
intensity  from  sample  1  as  a  function  of  the  pump  beam 
intensity.  The  crosses  are  the  measurements,  and  the 
solid  line  is  a  fit  of 
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This  functional  form  describes  the  saturation  of  the 
luminescence  due  to  depletion  of  the  ground  state  for 
pumping  by  a  beam  with  a  Gaussian  spatial  profile. 
The  fit  gives  Is  =  870  W/cm2,  which  is  less  than  the 
value  of  1,270  W/cm2  calculated  using  Is  =  hv/or.  The 
nonlinearity  of  the  luminescence  vs.  pumping  can  be 
due  to  a  number  of  factors  including  not  only 
depopulation  of  the  ground  state,  but  also  upconversion, 
excited  state  absorption  (ESA),  and  possibly  other 
effects.  However,  for  the  purpose  of  subsequent 
calculations  these  effects  are  conveniently  expressed  in 
terms  of  an  effective  saturation  intensity  that  results 
from  the  fit  of  the  above  equation.  The  subsequent 
calculations  using  this  effective  saturation  intensity 
should  be  correct  to  first  order. 


Figure  4.  Luminescence  vs.  pump  intensity  in  sample  1. 


Using  the  value  for  Is  from  the  fit  in  Fig.  4,  we  can  refit 
the  beam  coupling  gain  as  shown  in  Fig.  5.  The  line  is 
a  fit  to  the  range  from  0  to  80  W/cm2  that  gives  n2  /a  = 
4.0  x  10'7  cm2/W.  The  new  fit  is  much  better  than  the 
previous  one  shown  in  Fig.  3.  Most  of  the  remaining 
discrepancy  at  high  pump  intensity  can  be  explained  by 
the  nonexponential  decay  of  the  luminescence  that  takes 
place  at  high  intensity  (see  section  on  spectroscopy  and 
kinetics).  The  fit  in  Figs.  2  and  3  assumed  that  the 
pump  and  probe  beams  were  plane  waves  of  infinite 
extent.  The  fit  in  Fig.  5  includes  the  effects  of  a 
Gaussian  beam  profile.  Although  the  effect  is  relatively 
small,  it  does  yield  a  slightly  better  fit  to  the  data. 
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Figure  5.  New  fit  to  the  beam  coupling  measurements  of  Fig. 
3,  using  Is  as  determined  from  a  fit  to  Fig.  4. 

Sample  2  Exhibits  Different  Saturation 

Figure  6  shows  measurements  of  the  luminescence  vs. 
pump  intensity  for  sample  2.  The  fit  gives  Is  =  1640 
W/cm2,  which  is  greater  than  870  W/cm2  obtained  for 
sample  1.  Thus,  sample  2  does  not  show  premature 
saturation,  and  in  fact  the  saturation  appears  to  be  less 
than  that  predicted  by  Is  =  h v/gt.  One  possible 
explanation  for  the  measured  value  being  greater  than 
the  predicted  value  is  that  the  cross  section  may  contain 
terms  due  to  residual  color  center  absorption  that  do  not 
contribute  to  the  metastable  state  population. 


Figure  6.  Luminescence  vs.  pump  intensity  in  sample  2. 


Figure  7  shows  measurements  of  the  beam  coupling 
gain  in  sample  2.  The  dashed  line  is  a  fit  using  the 
value  of  n2  la  from  the  previous  crystal,  and  I  s=  1,270 
W/cm2.  It  fits  well  at  low  intensity,  but  it  greatly 
underestimates  the  gain  obtained  at  high  intensity.  The 
solid  line  is  a  fit  that  gives  n2! a  =  4.1  x  10’7  cm3/W 
and  Is  =  1,630  W/cm2.  These  values  are  in  excellent 
agreement  with  the  value  of  n2  fa  from  the  previous 
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Figure  7.  Gain  vs.  pump  intensity  in  sample  2. 

crystal,  and  the  value  of  Is  from  the  fit  of  the 
luminescence  for  sample  2. 

To  conclude  this  first  section,  by  fitting  the 
luminescence  vs.  pump  intensity,  we  obtain  an  effective 
saturation  intensity  that  lumps  together  the  effects  of 
ground  state  depletion,  upconversion,  ESA,  etc.  Using 
this  effective  saturation  intensity,  we  obtain  a  better  fit 
to  the  beam  coupling  gain  than  is  possible  using  the 
expression  Is  =  hv/ax.  Thus,  the  processes  that  give 
rise  to  saturation  of  the  luminescence  and  the 
metastable  state  population,  rather  than  the  beam 
breakup  that  we  had  previously  proposed  [1],  are  the 
major  factors  responsible  for  the  premature  saturation 
of  the  beam  coupling  gain  in  sample  1 . 

Spectroscopy  and  Kinetics 

We  studied  the  spectroscopic  and  kinetics  properties  of 
sample  1  (the  sample  exhibiting  premature  saturation 
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Figure  8.  Luminescence  spectra  for  strong  pumping. 
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Figure  9.  Luminescence  kinetics  taken  at  A,  =750  nm.  1  -  weak 
pumping,  2  -  strong  pumping  (409  W/cm2). 


behavior).  The  room  temperature  luminescence  spectra 
taken  at  high  (  10  4  W/cm2)  pumping  density  is 
presented  in  Fig.  8  for  a  pump  wavelength  of  514.5  nm. 
According  to  Ref.  [9],  two  narrow  lines  at  -723  nm  and 
725.5  nm  are  attributed  to  the  transition  from  the  2E 
level  to  the  ground  state  4A2  (R  lines),  and  the  irregular 
shaped  broad  luminescence  band  centered  at  -740-750 
nm  is  due  to  vibronic  transitions. 

Using  long  (200  ms)  pulsed  excitation  of  the 
sample  with  488  nm  or  514.5  nm  light  from  an  Ar  laser, 
we  studied  the  luminescence  decay  kinetics  for  each 
wavelength  at  which  there  was  a  significant  peak  in  the 
luminescence  spectrum.  For  long  pulse  excitation,  the 
decay-kinetics  taken  at  most  of  the  peaks  in  the 
luminescence  spectra  were  nearly  exponential  with  a 
time  constant  of  36  ms  (Fig.  9  trace  1).  This  value  is 
equal  to  the  room  temperature  lifetime  of  Cr  ions  in 
YAIO3  reported  in  Ref.  [9]. 

As  shown  in  Fig.  10,  the  decay-kinetics  taken  at  732  nm 
were  strongly  nonexponential,  and  the  time  taken  to 
reach  the  1/e  level  was  equal  to  24  ms. 

Although  the  decay  kinetics  at  725.5  nm  and  750 
nm  were  very  similar  for  long  pulsed  excitation,  they 
were  different  for  short  (2  ms)  pulse  excitation.  For 
short  pulsed  excitation,  all  the  kinetics  studied  were 
nonexponential  to  some  degree.  This  nonexponential 
behavior  was  almost  negligible  at  750  nm,  but 
significant  at  725.5  nm. 

These  spectroscopic  and  kinetics  studies  imply  the 
presence  of  at  least  three  Cr  luminescence  centers  in 
Cr:YA103,  that  are  responsible  for  the  dominant 


Figure  10.  Luminescence  kinetics  taken  at  A=732  nm.  1  - 
weak  pumping,  2  -  strong  pumping  (409  W/cm2). 


Figure  11.  Luminescence  intensity  at  A=725.5  nm  (squares) 
and  A=750  nm  (triangles)  as  a  function  of  pumping  density; 
Apump=488  nm.  The  solid  line  represents  the  slope  of  the 
725.5  nm  luminescence  curve  obtained  from  separate 
measurements  taken  at  low  pumping,  and  the  dashed  line 
represents  the  same  for  750  nm. 

features  of  the  luminescence  spectra  at  725.5  nm,  732 
nm,  and  750  nm.  The  presence  of  three  types  of  Cr 
centers  in  YAIO3  was  also  reported  in  Ref.  [10]. 
According  to  Ref.  [9],  Cr  decay  curves  are 
nonexponential  due  to  the  quenching  of  the 
luminescence  of  single  Cr  ions  on  Cr  pairs  [11]. 
However,  in  the  presence  of  several  types  of  Cr  centers, 
we  cannot  exclude  the  possibility  of  energy  transfer 
between  them,  which  can  also  lead  to  nonexponential 
luminescence  kinetics. 

The  dependence  of  the  luminescence  intensity  on 
the  pumping  density  is  shown  in  Fig.  11.  Similar 
nonlinear  dependencies  were  obtained  at  725.5  nm  and 
750  nm,  as  well  as  for  the  spectrally  integrated  emission 
that  was  shown  in  Fig.  4.  In  general,  the  nonlinear 
dependence  can  be  explained  by  depletion  of  the  ground 
state,  excited  state  absorption  (ESA),  or  energy  transfer 
upconversion  that  shortens  the  effective  life-time  of  the 
Cr  metastable  level.  In  Fig.  1 1,  as  was  the  case  for  the 
spectrally  integrated  emission  shown  in  Fig.  4,  the 
nonlinearity  of  the  luminescence  intensity  is  stronger 
than  predicted  by  a  model  neglecting  ESA  and 
upconversion. 

The  study  of  Cr  luminescence  kinetics  have  shown  that 
at  A  750  nm  and  many  other  wavelengths  except  725.5 
nm  and  732  nm,  upconversion  is  weak  (the  effective 
decay-time  changes  by  no  more  than  for  5%)  for 
pumping  densities  as  high  as  I/S  5x10  3  W/cm2.  At 
A  725.5  nm,  upconversion  is  significant  for  I/S  5x10  3 
W/cm2,  but  negligible  for  I/S  10  3  W/cm2.  The 
strongest  upconversion  was  found  at  732  nm.  As  shown 
in  Fig.  10,  at  this  wavelength,  the  initial  decay-time  of 
the  luminescence  changed  from  9  ms  to  6.5  ms  at  I/S 
=410  W/cm2.  At  all  wavelengths  studied,  a  higher  rate 
of  quenching  at  weak  pumping  corresponded  to 
stronger  upconversion  at  intense  pumping. 
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Assuming  that  different  types  of  Cr  centers  in  the 
crystal  behave  independently,  the  strong  nonlinearity  of 
the  luminescence  intensity  (at  725.5  nm  and  X  750  nm 
as  well  as  the  spectral  integrated  luminescence 
intensity)  vs .  pumping  density  at  I/S<1000  W/cm2 
cannot  be  explained  by  upcon version. 

We  obtained  an  unexpected  result  from  the 
luminescence  spectral  measurements  at  high  pumping 
intensity.  Accounting  for  the  upconversion  kinetics 
measurements,  we  expected  the  relative  intensity  of  the 
luminescence  peak  at  732  nm  to  be  smaller  at  high 
pumping  density  than  at  weak  pumping.  On  the 
contrary,  in  the  experiment,  the  relative  height  of  the 
732  nm  peak  was  greater  at  strong  pumping  than  at 
weak  excitation.  This  result  cannot  be  explained  under 
the  assumption  of  independent  behavior  of  different 
types  of  centers  as  well. 

In  principle,  the  last  experimental  result  can  be 
explained  if  we  assume  that  one  type  of  centers 
(manifesting  itself  at  732  nm)  absorbs  most  pump 
power  and  then  transfers  it  to  the  other  types  of  centers. 
In  this  case,  the  reduction  in  the  excited  state 
concentration  of  732  nm  centers  due  to  upconversion 
will  cause  the  reduction  in  the  energy  transferred  to  the 
other  centers.  To  explain  the  growth  of  the  relative 
height  of  732  nm  peak  at  intense  pumping,  we  should 
also  assume  that  some  additional  upconversion  (that  can 
be  rather  weak)  reduces  the  effective  decay  times  of 
luminescence  at  725.5  nm,  750  nm,  etc.  If  we  accept 
this  model,  then  at  I/S  400  W/cm2,  upconversion  can 
explain  up  to  a  30%  reduction  of  Cr  excited  state 
concentration.  This  model  also  is  consistent  with  the 
strong  quenching  of  732  nm  luminescence  at  weak 
pumping.  However,  at  the  present  moment  we  do  not 
have  enough  experimental  data  to  prove  or  reject  this 
hypothesis. 

Conclusion 

Laser  and  related  materials  play  an  important  role  in 
nonlinear  optics  by  providing  convenient  solid  state 
materials  with  a  local  response  proportional  to  intensity, 
which  allows  researchers  to  demonstrate  large  local 
nonlinear  optical  effects  using  cw  lasers.  For  example, 
using  Cr:YAlC>3  we  demonstrated  a  beam  coupling  gain 
of  22,  and  self-oscillation  based  on  this  gain  [2].  Since 
the  same  optical  nonlinearity  will  affect  the  propagation 
of  beams,  it  should  be  considered  in  the  design  of 
lasers. 

One  of  the  two  crystals  of  Cr:YA103  we  studied 
shows  a  saturation  in  both  beam  coupling  gain  and 
luminescence  as  a  function  of  pump  intensity  that  is 
greater  than  that  predicted  using  Is  =  hv/ax  and  known 


values  of  the  parameters.  In  this  crystal,  although  the 
luminescence  decay  kinetics  at  most  wavelengths  were 
nearly  exponential,  they  were  strongly  non-exponential 
and  intensity  dependent  at  732  nm.  The  second  crystal 
we  studied  did  not  exhibit  the  premature  saturation 
behavior.  If  the  premature  saturation  behavior  of  the 
first  crystal  is  not  linked  to  its  higher  concentration, 
then  crystals  like  the  second  one,  but  of  higher 
concentration,  may  provide  very  large  beam  coupling 
gain. 

This  research  was  sponsored  by  the  U.S.  Air  Force 
Office  of  Scientific  Research  under  contract  F49620- 
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Abstract 

Continuous  self-oscillation  and  self-pumped  phase  conjugation  via  nondegenerate  two- 
wave  mixing  gain  are  demonstrated  in  Cr:YA103.  To  our  knowledge,  this  is  the  first 
demonstration  of  these  effects  in  a  solid  state,  non-photorefractive  material. 
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Summary 

When  two  beams  of  different  frequency  interact  in  a  medium  with  a  local  nonlinear 
optical  response,  one  beam  experiences  gain  at  the  expense  of  the  other.  We  refer  to  this  process 
as  nondegenerate  two-wave  mixing.  Figure  1  shows  measurements  of  the  nondegenerate  two- 
wave  mixing  gain  for  a  weak  probe  beam  interacting  with  a  strong  pump  beam  in  a  crystal  of 
Cr:YA103.  The  laser  wavelength  is  515  nm  and  the  probe  is  downshifted  with  respect  to  the 
pump  by  approximately  5  Hz.  Gain  in  excess  of  12  times  is  obtained. 

At  normal  incidence,  the  reflectivity  of  Cr:YA103  is  11%,  and  the  absorption  loss  for  this 
crystal  is  17%.  Thus,  the  roundtrip  gain  (122  =  144)  is  sufficient  to  overcome  the  losses  ( 
(O.llxO.83)2  =  1/120)  of  a  resonator  formed  by  the  opposing  surfaces  of  the  crystal.  (We 
previously  demonstrated  that  the  gain  for  a  reflection  grating  is  comparable  to  that  for  a 
transmission  grating.)  With  the  probe  off  and  the  crystal  normal  oriented  at  a  small  angle  with 
respect  to  the  incident  pump  beam,  an  oscillation  builds  up  between  the  opposing  surfaces  of  the 
crystal.  Figure  2  shows  far  field  patterns  of  the  oscillating  mode.  As  shown  in  Fig.  2,  it  is 
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possible  to  change  the  oscillating  mode  by  translating  the  crystal  perpendicular  to  the  pumping 
beam.  Using  a  microscope,  we  observe  a  moving  fringe  pattern  in  the  fluorescence  from  the 
interaction  region  of  the  crystal  that  is  formed  by  the  interference  between  the  beam  incident  on 
the  crystal  and  the  frequency  shifted  self-oscillation.  By  beating  the  oscillation  with  a  portion  of 
the  incident  beam,  we  measure  a  frequency  shift  of  6  Hz.  This  is  consistent  both  with  our 
measurements  of  the  frequency  shift  at  which  the  nondegenerate  two-wave  mixing  gain  was 
maximum,  and  with  the  theoretical  expression,  1/2  x,  where  x  is  the  response  time  of  the 
nonlinearity  (the  lifetime  of  the  metastable  state,  33  ms). 

The  bidirectional  oscillation  acts  as  pumping  beams  and  the  beam  incident  on  the  crystal 
acts  as  a  probe  beam  in  a  four-wave  mixing  process  to  generate  what  we  believe  is  the  self- 
pumped  phase-conjugate  of  the  beam  incident  on  the  crystal.  Figure  3  shows  the  resulting  phase- 
conjugate  return.  We  verified  that  the  return  was  along  a  direction  consistent  with  a 
retroreflection  of  the  incident  beam  and  that  it  built  up  with  the  oscillation.  This  is  consistent 
with  the  interpretation  of  a  phase-conjugate  return.  However,  we  have  not  yet  verified  the  phase- 
conjugate  nature  of  the  wavefront  by,  for  example,  demonstrating  distortion  correction.  Note 
that  the  return  is  elongated  in  the  direction  perpendicular  to  the  plane  of  incidence  defined  by  the 
pump  and  probe  beams.  A  similar  effect  was  previously  observed  in  photorefractive  ring 
conjugators.  This  research  was  sponsored  by  the  U.S.  Air  Force  Office  of  Scientific  Research 
under  contract  F49620-92-C-0023. 
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Figure  Captions 

Fig.  1  Nondegenerate  two-wave  mixing  gain  for  a  weak  probe  beam  interacting  with  a  strong 
pump  beam  in  a  crystal  of  Cr:YA103. 

Fig.  2  Far-field  patterns  corresponding  to  the  spatial  modes  of  a  self-oscillation  in  Cr:YA103. 
The  different  modes  result  for  different  translational  positions  of  the  crystal. 

Fig.  3  The  elliptical  spot  on  top  of  the  speckle  pattern  of  scattered  light  is  the  far-field  pattern  of 
the  phase-conjugate  return  generated  by  the  four-wave  mixing  interaction  with  the 
incident  beam  and  the  counterpropagating  beams  of  the  self-oscillation. 
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Abstract 


We  present  experiments  demonstrating  reciprocity  lor  conjugation  of  highly  scattered  light;  if  a  fraction  v  of  the  total  scattered 
igil  is  incident  on  a  phase-conjugate  mirror  having  reflectivity  p,  then  a  fraction  p  rj1  returns  back  through  the  scatterer  as  a 
phase-conjugate  wave,  with  the  remaining  fraction  p(r)-rf)  being  lost  due  to  scattering.  Low  efficiency  results  from  inability 
to  gather  a  large  fraction  of  scattering,  but  the  efficiency  is  higher  than  is  possible  without  conjugation 


Consider  the  case  of  light  incident  on  a  lossless, 
highly  scattering  media.  As  shown  in  Fig.  I  a,  if  a  phase 
conjugator  could  gather  and  reflect  all  of  the  scattered 
light,  then  by  the  principle  of  time  reversal,  all  of  the 
reflected  light  would  return  as  the  conjugate  [1-4]. 
However,  in  practice  it  is  difficult  for  an  optical  system 
to  gather  a  large  fraction  of  highly  scattered  light.  What 
is  the  efficiency  in  this  case  of  partial  phase  conjuga¬ 
tion?  According  to  the  reciprocity  theorem,  if  a  fraction 
V  of  the  total  scattered  light  is  incident  on  a  phase- 
conjugate  mirror  having  reflectivity  p,  then  a  fraction 
P  rf  returns  back  through  the  scatterer  as  a  phase- 


1  A 

P) 

|y 

Conjugator 


conjugate  wave,  with  the  remaining  fraction  p  ( 77  -  rf) 
being  lost  due  to  scattering  [5J.  This  case  is  shown  in 
Fig.  lb.  We  previously  demonstrated  the  principle  of 
reciprocity  as  it  applies  to  non-polarization  preserving 
phase  conjugation  of  light  from  multimode  fibers, 
where  when  a  fraction  1/2  of  the  incident  light  (one 
polarization  component)  is  reflected  from  a  phase-con¬ 
jugate  mirror  with  a  reflectivity  p,  a  fraction  p/4  returns 
as  the  conjugate  [6].  New  experimental  results  dem¬ 
onstrating  the  principle  of  reciprocity  for  the  case  where 
a  varying  fraction  of  the  fiber  modes  are  conjugated  are 
to  be  published  [7] .  For  multimode  fibers,  we  had  also 


Fit  .  (a  Light  incident  on  a  lossless,  highly  scattering  media  is  scattered  into  all  directions.  An  ideal  conjugator  gathers  and  reflects  all  of  the 
suilleied  light  with  unity  reflectivity  to  reproduce  the  incident  beam  with  perfect  efficiency.  In  this  case  all  of  the  reflected  light  would  return  as 
the  conjugate,  (b)  In  practice,  a  conjugator  gathers  a  finite  fraction  „  of  the  scattered  light.  According  to  the  reciprocity  theorem,  if  a  fraction 
r?  the  total  scattered  light  is  incident  on  a  phase-conjugate  mirror  having  reflectivity  p,  then  a  fraction  p  if  returns  back  through  the  scatterer 
ils  11  phase-conjugate  wave,  with  the  remaining  fraction  p  (77—  t?2)  lost  due  to  scattering. 
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Fig.  2.  Setup  used  to  measure  scattered  power  density  as  a  function 
of  angle. 


Fig.  3.  Scattered  power  density  as  a  function  of  angle  for  the  thick 
(3/16"  polypropylene)  scalterer  used  in  our  experiments.  The  meas¬ 
urements  indicated  by  the  circles  are  normalized  to  the  on-axis  power 
and  the  solid  line  is  a  lit  of  an  exponential  dependence. 


Fig.  4.  Scattered  power  density  as  a  function  of  angle  for  the  thin 
(20  pun  polypropylene)  scalterer  used  in  our  experiments.  The  scat¬ 
ter  for  the  thin  scatterer  is  nearly  isotropic,  except  for  the  cosine 
dependence  resulting  from  the  area  projection  of  a  flat  surface. 


previously  demonstrated  that  all  of  the  fiber  modes  can 
be  gathered  and  conjugated,  such  that  all  of  the  light 
reflected  by  a  phase-conjugate  mirror  will  return  as  the 
conjugate  [6].  Here,  however,  we  describe  the  first 
experiments,  to  our  knowledge,  demonstrating  the  prin¬ 
ciple  of  reciprocity  as  it  applies  to  the  case  of  highly 
scattered  light,  where  all  of  the  scattered  light  cannot 
be  gathered  and  conjugated. 

In  our  experiments  we  used  both  a  thick  (3/16" 
polypropylene)  and  a  thin  (20  fjim  polyethylene)  scat¬ 
tering  media  to  show  that  the  principle  of  reciprocity 
applies  to  at  least  two  scatterers  with  different  proper¬ 
ties,  and  therefore  to  lend  credibility  to  the  general 
application  of  reciprocity  to  other  scattering  media.  To 
characterize  the  two  scatterers  we  measured  the  scat¬ 
tered  power  density  as  a  function  of  angle  using  the 
setup  shown  in  Fig.  2.  The  beam  incident  on  the  scat¬ 
terer  was  100  mW  of  unfocused  light  from  an  argon 
laser  operating  at  515  nm.  Measurements  of  the  scat¬ 
tered  light  at  various  angles  were  taken  using  a  detector 
of  1  cm2  area,  placed  at  74  cm  from  the  scatterer.  These 
measurements  are  shown  as  circles  in  Fig.  3  and  Fig. 
4,  and  the  solid  lines  are  fits  to  the  data.  The  scatter  for 
the  thick  scatterer  has  an  exponential  dependence  on 
angle  with  a  fwhm  of  approximately  40°.  The  scatter 
for  the  thin  scatterer  is  nearly  isotropic,  except  for  the 
cosine  dependence  resulting  from  the  area  projection 
of  a  flat  surface.  For  both  scatters,  calculations  of  the 
integrated  scatter  indicate  that  losses  are  negligible. 

Fig.  5  shows  the  experimental  setup  used  to  measure 
the  efficiency  for  phase  conjugation  of  scattered  light. 


Detector 


Fig.  5.  Experimental  setup  used  to  measure  the  efficiency  for  phase 
conjugation  of  scattered  light.  Light  from  a  laser  is  focused  by  a  lens 
onto  the  surface  of  the  scatterer.  A  pair  of  condensers  gathers  and 
focuses  scattered  light  into  a  self-pumped  conjugator.  A  polarizer 
placed  between  the  condensers  selects  the  proper  polarization  for 
operation  of  the  conjugator.  The  conjugate  return  is  measured  by  a 
detector  as  a  function  of  the  amount  of  light  incident  on  the  conju¬ 
gator,  which  is  adjusted  with  a  variable  aperture. 
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Fig.  6.  Photograph  of  the  phase-conjugate  return  taken  from  a  posi¬ 
tion  corresponding  to  the  detector. 
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Fig.  7.  Phase-conjugate  reflectivity  as  a  function  of  the  square  of  the 
fraction  of  the  total  scattered  light  that  is  incident  on  the  conjugator, 
in  the  case  of  the  thick  scatterer.  The  circles  represent  the  measure¬ 
ments  and  the  solid  line  are  fits  of  Rvc  =  p  rp,  as  predicted  by  the 
reciprocity  theorem.  The  excellent  agreement  of  the  functional  form 
for  most  of  the  data  indicates  the  validity  of  reciprocity  theorem  for 
predicting  the  results  of  partial  phase  conjugation  of  scattered  light. 
The  saturation  of  the  phase-conjugate  reflectivity  is  thought  to  be 
due  to  aberrations  in  our  condensers  at  low  1-number  that  do  not 
permit  the  focusing  of  all  of  the  light  into  the  conjugator  and/or 
spatial  variations  in  the  phase-conjugate  reflectivity  when  the  cone 
angle  of  the  light  incident  on  the  conjugator  becomes  large. 


Light  from  an  argon  laser  at  515  nm  is  focused  by  a 
lens  of  1  m  focal  length  onto  the  surface  of  the  scatterer. 
A  pair  of  condenser  lenses  (32.5  mm  diameter,  23.5 
mm  focal  length)  is  used  to  gather  and  focus  the  scat¬ 
tered  light  into  a  crystal  of  barium  titanate  that  acts  as 
a  self-pumped  phase  conjugator  of  the  “cat”  confi- 
guartion  [4] .  The  distance  between  the  first  condenser 
and  the  scatterer  is  adjusted  to  collimate  the  light  exit¬ 
ing  from  first  condenser  as  well  as  is  possible.  The  laser 
power  is  adjusted  to  have  approximately  14  mW  inci¬ 
dent  on  the  crystal.  A  polarizer  placed  between  the 
condensers  insures  that  the  light  incident  on  the  crystal 
is  extraordinary.  The  conjugate  return  is  reflected  from 
a  beam  splitter  and  measured  by  a  detector.  A  photo¬ 
graph  of  the  conjugate  return  at  a  position  correspond¬ 
ing  to  the  detector  (approximately  2  m  from  the 
scatterer)  is  shown  in  the  Fig.  6.  A  variable  aperture 
placed  between  the  condensers  allows  us  to  vary  the 
amount  of  light  incident  on  the  phase  conjugator.  When 
the  measurements  of  the  phase-conjugate  return  were 
completed,  the  conjugator  was  removed  and  replaced 
with  a  detector  to  measure  the  amount  of  light  passing 
through  the  aperture,  which  was  previously  incident  on 
the  conjugator,  as  a  function  of  the  aperture  setting.  A 
calibration  for  100%  reflectivity  was  made  by  tempo¬ 
rarily  placing  a  mirror  at  the  position  of  the  scatterer 
and  aligning  it  to  retrorellect  the  light  from  the  laser 
back  to  the  detector. 

The  circles  in  Figs.  7  and  8  show  the  measurements 
of  the  phase-conjugate  reflectivity  as  a  function  of  the 
square  of  the  fraction  of  the  total  scattered  light  that  is 


(Power  at  Crystal  /  Incident  Power)2  Tj2 


Fig.  S.  Phase-conjugate  reflectivity  as  a  function  of  the  square  of  the 
fraction  of  the  total  scattered  light  that  is  incident  on  the  conjugator, 
in  the  case  of  the  thin  scatterer. 
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incident  on  the  conjugator  rf ,  and  the  solid  lines  are 
fits  of  Rpc  =  p  rj2,  as  predicted  by  the  reciprocity  theo¬ 
rem.  Representative  error  bars  are  given  for  two  of  the 
points  in  both  Fig.  7  and  Fig.  8.  The  excellent  agree¬ 
ment  of  the  functional  form  for  most  of  the  data  indi¬ 
cates  the  validity  of  reciprocity  theorem  for  predicting 
the  results  of  partial  phase  conjugation  of  scattered 
light.  The  constant  p  from  the  fits  is  the  phase-conjugate 
reflectivity  of  the  self-pumped  crystal  of  barium  titan- 
ate.  The  values  of  33%  and  26%  are  in  good  agreement 
with  an  independent  measurement  of  30%  using 
unscattered  light. 

We  are  not  sure  why  the  phase-conjugate  reflectivity 
saturates,  but  it  may  be  due  to  residual  aberrations  in 
our  condensers  at  low  f-number  that  do  not  permit  the 
focusing  of  all  ot  the  light  into  the  conjugator  and/or 
spatial  variations  in  the  phase-conjugate  reflectivity 
when  the  cone  angle  of  the  light  incident  on  the  con¬ 
jugator  becomes  large.  The  effective  f-number  at  which 
the  saturation  behavior  takes  place  is  quite  low,  and  the 
corresponding  cone  angle  is  rather  large.  Saturation 
begins  at  about  f/ 1 .3  or  42°  in  Fig.  7,  and  17  1 .0  or  53° 
in  Fig.  8.  The  fact  that  the  data  agree  with  reciprocity 
to  such  low  f-numbers  and  large  cone  angles  is  satis¬ 
fying  in  and  ol  itself.  It  is  probably  not  too  surprising 
that  the  phase  conjugate  reflectivity  saturates  al  lower 
f-numbers  and  larger  cone  angles.  After  all,  neither  the 
condensers,  nor  the  conjugator  are  ideal.  For  example, 
the  condensers  are  manufactured  for  light  gathering, 
not  for  imaging,  and  the  residual  aberrations  are  severe. 
The  supplier  of  the  condensers  believes  that  they  can 
be  modeled  by  an  odd  asphere  with  second  order  cor¬ 
rection  for  the  first  surface,  and  a  120  mm  radius  for 
the  second  surface.  This  would  yield  a  focal  length  that 
varies  by  about  3  mm  (or  13%  of  the  specified  focal 
length)  over  the  zones  of  the  lens  aperture.  An  ideal 
conjugator  will  correct  for  aberrations  and  this  is  veri¬ 


fied  by  the  lact  that  most  of  the  data  agree  with  reci¬ 
procity.  However,  in  practice  conjugators  have  both 
finite  spatial  and  angular  apertures,  and  their  operation 
depends  on  the  focusing  conditions.  The  reflectivity  as 
a  function  of  the  angles  of  incidence,  in  one  plane,  on 
a  cat  conjugator  for  a  single  narrow  beam  [8]  and  for 
a  pair  of  narrow  beams  [9]  has  been  measured,  but  it 
is  not  clear  that  these  apply  to  our  situation  of  a  single 
extremely  wide  angle  beam  in  two  dimensions  with  our 
particular  aberrations. 

In  conclusion,  for  any  reciprocal  optical  system,  if  a 
fraction  r)  of  the  incident  light  is  incident  on  a  phase- 
conjugate  mirror  having  reflectivity  p,  then  a  fraction 
p  rj2  returns  as  a  conjugate,  with  the  remaining  fraction 
P  ( V~V2)  distributed  over  other  modes  of  the  system. 
Our  experimental  results  demonstrate  this  principle  of 
reciprocity  in  one  of  the  most  complicated  cases  (the 
case  of  scattering). 

This  research  was  sponsored  by  the  Office  of  Naval 
Research  under  contract  N000 1 4-93-C-00 1 0  and  by  the 
U.S.  Air  Force  Office  ot  Scientific  Research  under  con¬ 
tract  F49620-92-C-0023. 
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A  gain  of  22  times  was  obtained  by  nondegenerate  two-wave  mixing  in  Cr:YAI03.  To  our  knowledge,  this  is  the 
largest  cw  two-wave  mixing  gain  obtained  in  a  bulk  solid-state  nonphotorefractive  material.  The  measured  gain 
appears  to  be  limited  by  beam  breakup  that  is  due  to  spatial  nonuniformities  in  the  nonlinear  refractive-index 
change  that  are  the  result  of  inhomogeneities  in  the  crystal.  Predictions  based  on  our  data  indicate  that  gain 
in  excess  of  200  times  should  be  possible  in  Cr;YA103  if  a  homogeneous  crystal  can  be  obtained. 


Recently  we  demonstrated  a  gain  of  six  times  by  non¬ 
degenerate  two-wave  mixing  (NDTWM)  in  chromium- 
doped  yttrium  aluminate  (CriYAlOs).1  Here  we 
report  new  experiments  in  which  we  achieved  a  gain 
of  22  times  by  increasing  both  the  interaction  length 
and  the  doping.  To  our  knowledge,  this  is  the  largest 
cw  NDTWM  gain  obtained  in  a  solid-state  nonpho¬ 
torefractive  material.  The  large  gains  possible  in 
this  material  should  permit  the  development  of  cw 
oscillators  and  self-pumped  conjugators  with  proper¬ 
ties  different  from  those  employing  photorefractives. 

Experiments  demonstrating  high-gain  NDTWM  in 
new  solids  by  use  of  low-power  cw  lasers  constitute 
an  essential  step  in  establishing  the  practicality  of  al¬ 
ternative  materials  for  applications  of  nonlinear  op¬ 
tics  in  areas  such  as  signal  processing.  Availability 
of  high  gain  can  bring  about  rapid  progress  in  ex¬ 
perimental  research,  as  occurred  with  the  discovery 
of  photorefractives  such  as  barium  titanate  (BaTi03). 
To  our  knowledge,  however,  no  one  has  demonstrated 
gain  comparable  with  that  of  photorefractives  by  us¬ 
ing  a  solid-state  bulk  (as  opposed  to  guided-wave) 
Kerr  medium  with  a  cw  source.  We  report  proper¬ 
ties  of  a  candidate  material,  Cr:YA103,  that  may  sig¬ 
nificantly  advance  experimental  efforts  on  nonlinear 
optics  in  Kerr-like  media,  much  like  the  way  in  which 
basic  studies  of  BaTi03  initiated  widespread  study  of 
nonlinear  optics  in  photorefractives.  Motivation  for 
this  study  was  provided  by  the  expectation  of  new 
properties  that  should  result  from  the  different  non¬ 
linear  mechanism  (Kerr-like  versus  photorefractive), 
such  as  improved  dynamic  range  in  mutually  pumped 
phase  conjugation,2  that  is  important  in  wave-front- 
matched  heterodyne  receivers.3 

In  some  metal-ion-doped  insulators,  such  as  ruby, 
metastable  states  result  in  low  saturation  intensi¬ 
ties  and  correspondingly  large  optical  nonlinearities. 
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The  origin  of  the  nonlinearity  was  previously  deter¬ 
mined  to  be  the  result  of  light-induced  population 
changes  and  the  difference  in  the  susceptibility  be¬ 
tween  the  ground  and  metastable  states.4  Interest 
in  these  materials  stems  from  the  fact  that  they 
are  convenient  solid-state  materials  with  large  op¬ 
tical  nonlinearities  that  can  be  utilized  with  rela¬ 
tively  low-power  cw  lasers.  Photorefractives  also 
meet  these  conditions  and  have  demonstrated  two- 
wave  mixing  gain  of  >  103,  but  the  differences  in  the 
nonlinear  mechanisms  make  investigation  of  the 
metal-ion-doped  insulators  of  interest  because  of 
their  potential  advantages.  As  one  example,  the 
beam  ratio  for  which  mutually  pumped  phase  conju¬ 
gation  will  take  place  in  photorefractives  is  limited, 
but  this  is  not  the  case  in  Kerr  materials.2  Our 
search  of  the  literature  addressing  nonlinearities  in 
metal-ion-doped  insulators  indicated  that  Cr;YA103 
stands  out  as  the  material  having  the  largest  poten¬ 
tial  index  change  for  a  given  amount  of  absorption, 
i.e.,  the  largest  value  of  n2Is/ot,  where  Is  is  the  satu¬ 
ration  intensity,  n2  is  the  optical  Kerr  coefficient,  and 
a  is  the  absorption  coefficient.5  We  therefore  chose 
this  material  for  experimental  attempts  to  obtain 
high  two-wave  mixing  gain. 

As  grown,  Cr:YA103  has  a  color-center  absorption 
that  competes  with  the  absorption  responsible  for 
the  optical  nonlinearity.  To  remove  the  color  cen¬ 
ters,  we  heated  as-grown  material  in  a  forming  gas 
atmosphere.6  Figure  1  shows  absorption  spectra  for 
the  12-mm-thick  end  of  a  boule  of  Cr:YA103  as  grown 
(dotted  curve),  after  a  first  anneal  (dashed  curve)  in 
95%  N2/5%  H2  at  1050  °C  for  8  h  with  a  flow  rate 
of  250  cm3/min  and  a  ramp  rate  of  1  °C/min,  and  af¬ 
ter  a  second  anneal  (solid  curve)  at  1300  °C  for  24  h 
with  the  same  atmosphere  and  the  same  flow  and 
ramp  rates  as  used  in  the  first  anneal.  After  the 
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Fig.  1.  Cr:YA103  absorption  spectra  showing  the  bleach¬ 
ing  of  color  centers  by  heat  treatment. 

second  anneal  the  boule  was  fabricated  into  several 
crystallographic-cut  cubes,  one  of  which  was  used  in 
these  experiments.  Microprobe  measurements  indi¬ 
cated  that  the  Cr  concentration  was  0.12  wt.%. 

Figure  2  is  a  simplified  schematic  of  the  experi¬ 
ment  used  to  measure  NDTWM  gain.7  Light  from 
an  Ar  laser  at  514.5  nm  is  split  by  a  beam  splitter 
into  a  weak  probe  beam  /2(0)  and  a  strong  pump 
beam  /i(0).  We  shift  the  frequency  of  the  pump  with 
respect  to  the  probe  by  reflecting  it  from  a  mirror 
mounted  to  a  piezoelectric  transducer  and  driven  by  a 
triangle-wave  voltage  source.  This  results  in  a  mov¬ 
ing  interference  pattern  in  the  crystal  and  a  popula¬ 
tion  grating  that  lags  in  phase  because  of  the  finite 
response  time.  As  in  photorefractives,  this  phase 
shift  is  responsible  for  the  energy  exchange.  The 
sign  of  the  phase  shift  depends  on  the  sign  of  the 
frequency  shift  and  determines  the  direction  of  en¬ 
ergy  flow.  Amplification  of  the  probe  results  when 
the  mirror  moves  so  that  the  pump  is  upshifted, 
and  attenuation  results  when  the  mirror  moves  in 
the  opposite  direction  so  that  the  pump  is  down¬ 
shifted.  The  triangle-wave  frequency  and  thus  the 
pump  frequency  was  varied  to  maximize  gain.  This 
maximum  occurs  for  a  frequency  shift  of  ~5  Hz,  con¬ 
sistent  with  the  theoretical  behavior  described  by7,8 
1/(27 tt)  =  4.8  Hz,  where  r  is  the  metastable-state 
lifetime  of  33  ms.  Steady  state  is  achieved  since  the 
period  of  our  triangle  wave  (~7  s)  greatly  exceeds  r. 

To  obtain  the  intensities  required  for  high  gain, 
we  focus  the  beams  to  a  1/e2  waist  of  210-/zm 
radius  at  the  crystal.  We  ensure  a  high  degree 
of  overlap  in  the  crystal  length  L  —  9.1  mm  by 
having  the  beams  intersect  at  the  small  angle 
of  0.02  rad.  The  beams  were  incident  nearly 
parallel  to  the  (001)  direction  and  s  polarized 
perpendicular  to  the  (110)  direction.  This  po¬ 
larization  and  crystal  orientation  yield  the  high¬ 
est  peak  in  the  absorption  spectrum  at  560  nm 
and  the  lowest  minimum  in  the  residual  ab¬ 
sorption  between  peaks  at  475  nm.  Thus  we 
believe  that  this  configuration  yields  the  greatest 
interaction  with  the  absorption  responsible  for  the 
optical  nonlinearity  and  the  least  interaction  with 
residual  color  centers.  Depletion  of  the  pump  is 
avoided  by  use  of  a  large  input  beam  ratio,  i.e., 

/i(0)//2(0)  =  1000. 

Figure  3  shows  the  gain  I2(L)/I2( 0)  as  a  function  of 
the  pump  intensity.  The  solid  curve  is  a  fit  of 


h{L) 

h(0) 


=  exp  (-aL) 


X  exp 


1  -  exp(-nL) 
aL 


27rn2Ii(0)L/ A  \ 
[1  +  /1(0)//sp  \)’ 
(1) 


with  aL  =  0.69,  A  =  514.5  nm,  Is  =  kv/crr  = 
1.27kW/cm2,  where  h  is  Planck's  constant,  v  is 
the  optical  frequency,  a  is  the  absorption  cross  sec¬ 
tion,  and  t  is  the  metastable-state  lifetime.  With 
n2  as  the  only  adjustable  parameter,  the  fit  yields 
n2  =  3.2  X  10-7  cm2/W.  The  nonlinearity  that  is 
due  to  population  effects  dominates  in  both  beam 
coupling  and  self-focusing  because  the  n2  s  for  the 
corresponding  thermal  effects  are  at  least  two  orders 
of  magnitude  lower.  We  can  derive  Eq.  (1)  from 
Ref.  8  by  including  the  effect  of  saturation,  which  re¬ 
duces  the  nonlinear  index  by  a  factor  [1  +  Ii(0)/Is]~2, 
for  7i(0)  »  J2(0). 

The  inset  in  Fig.  3  shows  a  plot  of  the  gain  for  the 
range  from  zero  to  400  W/cm2.  With  the  previously 
stated  parameters  and  rt2  from  the  fit  for  the  lower- 
intensity  range,  the  curve  in  the  inset  represents  a 
theoretical  projection  based  on  Eq.  (1).  A  gain  of 
22  times  was  obtained  at  approximately  400  W/cm2, 
but  the  measured  gain  does  not  increase  with  inten¬ 
sity  as  much  as  theory  predicts.  (Theory  predicts  a 
gain  of  200  at  400  W/cm2.)  We  believe  this  is  due 
to  beam  breakup  resulting  from  spatial  nonuniformi¬ 
ties  of  the  nonlinear  refractive-index  change  in  the 
interaction  region.  As  evidence,  Fig.  4  shows  the 
transmitted  pump  with  the  probe  off  (a)  at  low  in¬ 
tensity  and  normal  exposure  (0.6  W/cm2,  1/100  s) 
indicating  the  transmitted  beam  size,  (b)  at  low 
intensity  and  at  a  longer  exposure  (1.4  W/cm2,  1  s) 
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Fig.  2.  Simplified  schematic  of  the  experiment  used  to 
measure  NDTWM.  PZT,  piezoelectric  transducer. 


Fig.  3.  Gain  versus  pump  intensity.  The  solid  curve  is 
a  fit  of  Eq.  (1).  The  inset  shows  that,  at  high  intensity, 
the  measured  gain  does  not  increase  with  intensity  as 
much  as  theory  predicts.  We  believe  this  is  due  to  beam 
breakup  resulting  from  spatial  nonuniformities  of  the 
nonlinear  refractive-index  change. 
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(a)  (b)  (C) 

Fig.  4.  Transmitted  pump  with  the  probe  off  (a)  at  low 
intensity  and  normal  exposure  indicating  the  transmit¬ 
ted  beam  size,  (b)  at  low  intensity  and  long  exposure 
showing  an  X-shaped  pattern  of  linear  scattering  from 
crystal  striae,  and  (c)  at  high  intensity  and  an  exposure 
equivalent  to  that  in  ( b),  showing  the  beam  breakup.  For 
all  photographs  the  plane  of  incidence  is  the  horizontal. 

to  reveal  the  X-shaped  pattern  of  weak  linear  scat¬ 
tering  that  is  believed  to  be  the  result  of  scatter¬ 
ing  by  crystal  striae,  and  (c)  at  high  intensity  and 
an  exposure  equivalent  to  that  in  (b)  (140  W/cm2, 
1/100  s)  showing  the  nonlinear  scattering  that  is 
due  to  nonuniformities  in  the  nonlinear  refractive- 
index  change. 

Spatial  nonuniformities  of  the  nonlinear  refractive- 
index  change  in  the  interaction  region  can  be  the  re¬ 
sult  of  shadows  or  distortions  from  surface  or  bulk 
defects  or  from  spatial  variations  in  the  Cr-ion  con¬ 
centration.  Photographs  of  the  crystal  illuminated 
by  a  low-intensity  collimated  beam  show  radial  and 
concentric  patterns  thought  to  be  Cr  striae.  Images 
of  a  low-intensity  pump  beam  inside  the  crystal  also 
indicate  that  the  striae  have  a  significant  effect  on  the 
beam  profile.  Beam  breakup  can  result  both  from 
variations  in  Cr  concentration,  yielding  a  variation 
in  n2,  and  from  intensity  variations,  yielding  varia¬ 
tions  in  the  induced  nonlinear  index  change. 

Qualitative  behavior  indicated  that  the  scattering 
at  high  intensity  was  nonlinear,  but  to  confirm  this 
we  measured  the  intensity  of  a  small  region  of  the 
scattering  as  a  function  of  the  pump  intensity.  An 
aperture  was  closed  down  to  a  diameter  such  that, 
at  low  intensity,  it  would  pass  one  half  of  the  pump 
that  transmitted  through  the  crystal,  and  then  it 
was  placed  in  the  speckle  of  the  scatter  indicated 
by  the  arrow  in  Fig.  4(c).  Measurements  indicated 
that  the  nonlinear  portion  of  the  scatter  had  a  cubic 
dependence  on  the  pump  intensity. 

We  made  the  following  observations  concerning  the 
scattering:  (1)  the  strength  of  the  nonlinear  scat¬ 
tering  depended  on  feedback  into  the  crystal  and  (2) 
a  threshold  existed  above  which  the  speckle  pattern 
of  the  scattering  was  dynamic.  An  experiment  ex¬ 
amining  the  interference  between  a  portion  of  the 
transmitted  pump  and  a  speckle  did  not  detect  any 
frequency  shift.  This  indicated  that  the  dynamic 
scattering  at  high  intensities  was  not  due  to  beat¬ 
ing  between  frequency- shifted  stimulated  scattering 
resulting  from  NDTWM  gain  and  some  linear  scat¬ 
tering  from  crystal  inhomogeneities.  When  the  crys¬ 
tal  was  oriented  to  retroreflect  the  pump  or  when  an 
external  mirror  was  used,  the  strength  of  the  nonlin¬ 
ear  scattering  increased  and  the  threshold  of  dynamic 
scattering  decreased.  To  obtain  the  highest  gains  in 
the  NDTWM  measurements  we  rotated  the  crystal 
away  from  normal  to  reduce  the  nonlinear  scatter¬ 


ing.  An  angle  of  incidence  for  the  pump  of  0.03  rad 
yielded  the  largest  gains.  Further  reductions  in  the 
feedback  resulting  from  the  rear  surface  reflection 
might  improve  the  gain. 

In  conclusion,  we  demonstrated  a  gain  of  22  times 
by  NDTWM  in  Cr:YA103.  The  gain  appears  to  be 
limited  by  beam  breakup  resulting  from  spatial 
nonuniformities  in  the  nonlinear  refractive-index 
change  that  are  the  result  of  inhomogeneities  in 
the  crystal.  Theory  predicts  that  gain  in  excess 
of  200  times  should  be  possible  in  a  homogeneous 
crystal.  Thus,  like  photorefractives,  Cr:YA103  is  a 
convenient  solid  having  the  potential  for  high  gain 
by  use  of  cw  lasers.  Interest  in  Cr:YA103  stems 
from  the  difference  in  nonlinearity  (Kerr-like  versus 
photorefractive)  and  the  resulting  potential  for  new 
properties  rather  than  from  any  improvement  in  sen¬ 
sitivity  or  speed.  Indeed,  the  sensitivity  of  Cr:YA103 
(~10_n  m3/J)  is  2  orders  of  magnitude  below  that  of 
BaTi03.  Nor  is  Cr:YA103  particularly  fast;  the  grat¬ 
ing  decay  time  in  Cr:YA103  is  determined  by  the 
metastable  lifetime  of  33  ms,  and  this  is  compara¬ 
ble  with  the  fastest  grating  decay  times  obtained 
in  BaTi03  under  moderate  illumination.  Although 
speed  increases  with  intensity  in  photorefractives, 
it  is  not  clear  that  they  can  withstand  the  intensi¬ 
ties  used  with  Cr:YA103,  making  a  direct  comparison 
difficult.  Gain  is  determined  by  the  index  change, 
and  the  maximum  possible  index  change  in  Cr:YA103 
(Anmax  ~  n2Is  =  4  X  10"4)  is  comparable  with  that 
of  BaTi03.  In  addition  to  its  having  potential  for 
new  properties  resulting  from  the  different  nonlin¬ 
ear  mechanism,  we  note  that  Cr:YA103  is  only  one 
example  of  a  large  class  of  metal-ion-doped  insu¬ 
lators.  Admittedly,  it  has  the  best  performance  to 
date,  but  an  effort  to  develop  metal-ion-doped  insu¬ 
lators  as  nonlinear  materials  that  is  as  extensive  as 
that  undertaken  for  photorefractives  may  yield  new 
materials  with  improved  figures  of  merit. 
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Abstract 

We  study  the  problem  of  two  beams  phase  conjugating  each  other  via  mutually  generated 
reflection  gratings  in  an  electrostrictive  medium.  We  present  the  first  solutions  to  the 
transcendental  equation  that  illustrate  large  dynamic  range. 
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SUMMARY 

A  mutually  pumped  phase  conjugator  (MPPC)  generates  the  phase-conjugate  replicas  of 
two  incoherent  incident  beams.  Each  input  beam  is  converted  into  the  phase-conjugate  replica 
of  the  other  by  Bragg  diffraction  off  a  shared  grating.  In  photorefractive  media,  several 
configurations  differing  in  their  number  of  internal  reflections  from  the  crystal  surfaces  were 
demonstrated;1  for  efficient  operation,  the  two  incoherent  beams  must  have  comparable  input 
intensities.  However,  applications  like  phase-conjugate  heterodyne  detection2  require  the  device 
to  work  for  large  imbalance  of  input  beam  intensities,  when  a  weak  optical  signal  from  a  remote 
transmitter  is  combined  with  the  strong  beam  from  a  local  oscillator.  A  modest  extension  of  the 
dynamic  range  can  be  obtained  by  increasing  the  photorefractive  gain;1  however,  it  is  desirable 
to  extend  the  dynamic  range  of  the  device  to  several  orders  of  magnitude. 

MPPC  has  also  been  studied  in  electrostrictive  Kerr  media  by  utilizing  the  Brillouin  gain 
in  a  transmission  geometry.3  An  advantage  of  MPPC  in  electrostrictive  Kerr  media  over 
photorefractive  media  is  the  large  dynamic  range  of  input  beam  ratio  over  which  the  process  will 
occur.  This  is  because  if  one  beam  intensity  is  large  enough  to  satisfy  the  threshold  condition  for 
the  total  beam  intensity,  then  a  small  intensity  of  the  second  beam  will  initiate  the  process  of 
MPPC.  However,  none  of  this  work  has  been  corroborated  experimentally.  There  has  been 
recent  theoretical  and  experimental  work  on  MPPC  in  electrostrictive  Kerr  media  using  the 
reflection  geometry.4  That  large  dynamic  range  was  also  possible  in  reflection  geometry  was 
discussed  in  Ref.  4,  but  here  we  present  the  first  solutions  to  the  transcendental  equation  that 
illustrate  this  possibility. 


Electrostrictive 
Kerr  Medium 

(a)  (b) 


Fig.  1  (a)  Schematic  diagram  showing  MPPC  using  reflection  gratings  in  an  electrostrictive 
Kerr  medium,  (b)  Wave  vector  diagram  for  MPPC  using  reflection  gratings. 


The  geometry  for  the  nonlinear  interaction  responsible  for  MPPC  is  shown  in  Fig.l(a) 
and  Fig.l(b).  Input  beams  1  and  4  are  a  pair  of  mutually  incoherent  beams  that  enter  through  the 

front  face  of  the  nonlinear  medium  and  intersect  at  an  angle  20.  Acoustic  phonons  that  are 
initially  present  in  a  transparent  Kerr  medium  like  CS2  due  to  thermal  and  quantum  noise  cause 
Stokes  scattering  of  the  input  beams  in  all  directions.  Also,  a  photon  from  each  input  beam  can 
be  spontaneously  converted  into  a  frequency  downshifted  Stokes  photon  and  an  acoustic  phonon. 
The  coherently  generated  Stokes  waves  will  interfere  with  the  input  beam  to  produce  interference 
patterns  which  travel  at  the  acoustic  velocity,  and  which  drive  acoustic  waves  by 
electrostnction.5  These  acoustic  waves  induce  index  gratings  in  the  medium,  and  the  index 
grating  that  diffracts  each  incoherent  input  beam  into  the  phase-conjugated  output  of  the  other 
input  beam  beam  will  be  reinforced  by  both  the  beams.  Hence  this  mutual  index  grating  has  the 
maximum  gain  compared  to  all  the  other  possible  gratings,  and  we  label  the  corresponding 
Stokes  scattered  wave  of  input  beam  1  as  beam  3,  beam  3  being  the  phase-conjugate  of  input 
beam  4.  Similarly,  the  relevant  Stokes  wave  of  input  beam  4  is  beam  2,  and  this  is  the  phase- 
conjugate  of  pump  beam  1.  Coupled- wave  equations  for  the  complex  beam  amplitudes  Ai  (i=lto 
4)  describe  the  energy  exchanged  between  the  various  beams  at  steady-state: 

f  (A14  +  A+2A4)A3 


-  f  (aV^  +  a^  )A4 


“f  (A*1A3  +  A2<  )At 


(A1/f3  +  A*2A4)A2  ,  (1) 

where  g  is  the  Brillouin  gain  coefficient.  Assuming  non-zero  Stokes  seeds  at  the  z  =  L  backplane 
that  arise  from  the  mutual  scattering  of  each  input  beam  in  a  direction  counterpropagating  to  the 
other  input  beam:  A2L  =  f  A4(L),  A3L  =  f  A1(L),  f  being  the  mutual  scattering  coefficient  that  is 
equal  in  amplitude  and  phase  for  both  the  pumps,  we  obtain  the  following  transcendental 

equation  for  the  phase-conjugate  reflectivity  R  =  IA3(0)  /A40l2: 
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R  =  q  Ifl2  exp[G(l— R  /  q)]  .  (2) 

Here  q  is  the  intensity  ratio  of  the  two  incoherent  beams  at  the  z=  0  input  plane:  q  =  I10  / 140>  and 
G  =  S(Iio+^4o)L  is  die  total  Brillouin  gain.  The  phase-conjugate  reflectivity  S  in  the  other  arm  is 
equal  to  IA2(0)  /A*qI2  and  is  related  to  R  by:  S  =  R/q2.  The  device  transmissivity  T  measures  the 
fraction  of  each  input  beam  that  is  converted  into  the  phase-conjugate  output  of  the  other  beam, 
and  is  equal  for  both  inputs.  Note  that  if  Ifl2  =  0,  i.e.,  I2L  =  I3L  =  0,  then  the  phase-conjugate 
reflectivities  R  and  S  are  equal  to  zero.  Hence  finite  seed  values  of  the  Stokes  waves  are 
required  for  MPPC  to  work  in  reflection  geometry.  Fig.  2  is  a  plot  of  the  solutions  of  Eq.  (2) 
showing  the  phase-conjugate  reflectivities  R,  S  and  transmissivity  T  as  a  function  of  the  input 
intensity  ratio  q  on  a  log-log  scale,  for  a  fixed  value  of  Brillouin  gain  G4=18  satisfied  by  the 
strong  input  beam  4.  This  is  below  the  threshold  value  required  for  self-SBS  (  ~  23)  of  input 
beam  4.  We  take  Ifl2  =  10*12  for  noise  induced  Stokes  scattering.  Increasing  q  is  equivalent  to 


increasing  Iio.  keeping  I40  and  hence  G4  fixed.  Note  that  the  weak  input  beam  I^o  (q  =  10***)  is 
phase-conjugated  with  a  reflectivity  greater  than  unity.  (S  =  60).  Since  S  is  inversely 

proportional  to  Iio>  S  decreases  with  increasing  q  for  q  ^  10*k  As  q  is  further  increased,  the  total 
Brillouin  gain  increases,  leading  to  increased  conversion  of  input  beam  I40  into  the  Stokes  beam 
12(0),  as  seen  by  the  increase  in  T.  This  compensates  for  the  decrease  of  S  with  increasing  I10,  so 
that  S  actually  increases  for  a  limited  range  of  q  values.  When  the  total  Brillouin  gain  is  enough 
to  convert  most  of  the  input  beams  into  their  Stokes  waves,  then  any  further  increase  in  q  results 
in  decreasing  S  once  again.  R  is  a  monotonically  increasing  function  of  q  because  its 
denominator  I40  is  a  fixed  quantity,  and  increase  in  1 10  results  in  more  conversion  of  this  input 
beam  into  its  Stokes  wave  13(f)). 


10*  10*  10*  10*  10*  iff1  io°  101 


Input  Beam  Ratio  q  =  I  !i 

1  10  40 

Fig.  2  Phase-conjugate  reflectivities  R,  S  and  transmissivity  T  as  a  function  of  intensity  ratio  of 
the  two  input  beams  q,  for  Ifl2  =  10*12  and  G4  =  18. 

In  conclusion,  we  studied  Brillouin  induced  MPPC  in  the  reflection  geometry  and 
obtained  solutions  illustrating  that  a  weak  beam  can  be  phase-conjugated  with  reflectivities  larger 
than  unity  .  The  technique  works  for  a  large  dynamic  range  of  input  beam  ratio  as  the  intensity  of 
the  weak  input  beam  is  increased,  while  the  intensity  of  the  strong  incoherent  input  beam  is  kept 
constant  below  the  threshold  value  for  self-SBS. 
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Abstract 


A  gain  of  22  times  was  obtained  by  nondegenerate  two-wave  mixing  in  Cr:YA103. 
Theoretical  predictions  based  on  our  experimental  data  indicate  that  gain  in  excess  of  200  times 
should  be  possible  in  Cr:YA103  if  a  homogeneous  crystal  can  be  obtained. 
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Summary 

Recently  we  demonstrated  a  gain  of  6  times  that  was  obtained  by  nondegenerate  two- 
wave  mixing  in  CnYAlC^.1  Here  we  report  new  experiments  in  which  a  gain  of  22  times  was 
achieved  by  increasing  both  the  interaction  length  and  the  doping  of  the  material.  To  our 
knowledge,  this  is  the  largest  cw  two-wave  mixing  gain  obtained  in  a  bulk,  solid  state,  non- 
photorefractive  material.  Experiments  demonstrating  high  gain,  nondegenerate  two-wave  mixing 
in  new  solids  using  low  power  cw  lasers  constitute  an  essential  step  in  establishing  the  viability 
of  alternative  materials  for  applications  of  nonlinear  optics  in  areas  such  as  signal  processing. 
Availability  of  high  gain  can  bring  about  rapid  progress  in  experimental  research,  as  occurred 
with  the  discovery  of  photorefractive  materials  such  as  barium  titanate.  To  our  knowledge 
however,  no  one  has  demonstrated  gain  comparable  to  that  of  photorefractive  materials  by  using 
a  solid  state  bulk  (as  opposed  to  guided  wave)  Kerr  medium  with  a  cw  source.  We  report 
properties  of  a  candidate  material,  chromium-doped  yttrium  aluminate  (Cr:YA103),  that  may 
significantly  advance  experimental  efforts  on  nonlinear  optics  in  Kerr-like  media,  much  the  way 
basic  studies  of  barium  titanate  initiated  widespread  study  of  nonlinear  optics  in  photorefractive 
materials.  Motivation  for  the  present  work  was  provided  by  the  expectation  of  new  properties 
that  should  result  from  the  different  nonlinear  mechanism  (Kerr-like  vs.  photorefractive),  such  as 
improved  dynamic  range  in  mutually  pumped  phase  conjugation2  that  is  important  in  wavefront 
matched  heterodyne  receivers.3 

Figure  1  is  a  simplified  schematic  of  the  experiment  used  to  measure  gain  in 
nondegenerate  two-wave  mixing.4  Light  from  an  argon  laser  is  split  into  a  weak  probe  beam 
I2(0)  and  a  strong  pump  1,(0),  with  a  fixed  ratio  I2(0)/I,(0)  =  1/1000.  The  frequency  of  the  pump 
is  shifted  by  approximately  5  Hz  with  respect  to  the  probe  by  reflecting  it  from  a  mirror  mounted 
to  a  piezoelectric  transducer  and  driven  by  a  triangle-wave  voltage  source.  The  grating  formed 
by  the  interaction  of  the  pump  and  probe  in  the  crystal  of  Cr:YA103  results  in  amplification  for 
the  probe  when  the  mirror  is  moving  in  one  direction,  and  attenuation  when  the  mirror  moves  in 
the  opposite  direction. 
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Fig.  1  Schematic  of  experiment  used  to  measure  nondegenerate  two-wave  mixing  gain. 


Figure  2  (a)  shows  a  plot  of  the  gain  as  a  function  of  the  intensity  of  the  pump  for  the 
range  from  0  to  80  W/cm2.  The  circles  are  the  experimental  measurements,  and  the  line  is  a  fit 
of  the  theoretical  form. 
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with  aL  =  0.69,  L  =  9.1  mm,  X  =  514.5  nm  and  Is  =  hv/aT  =  1.27  kW/cm2,  where  h  is  Planck's 
constant,  v  is  the  optical  frequency,  a  is  the  absorption  cross-section  and  %  is  the  metastable  state 
lifetime.  Using  n2  as  the  only  adjustable  parameter,  the  fit  yields  n2  =  3.2  x  10'7  cm2/W. 


(a) 


(b) 


Fig.  2  Gain  vs.  pump  intensity:  (a)  from  0  to  70  W/cm2,  and  (b)  from  0  to  400  W/cm2. 
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Figure  2  (b)  shows  the  gain  for  the  range  from  0  to  400  W/cm2.  Using  the  previously 
stated  parameters  and  ti2  as  determined  from  the  fit  for  the  lower  intensity  range  of  Fig.  2(a),  the 
line  in  Fig.  2(b)  represents  projections  based  on  the  theoretical  form  expressed  by  Eq.  (1).  A 
gain  of  22  times  was  obtained  at  approximately  400  W/cm2,  but  the  measured  gain  does  not 
increase  with  intensity  as  much  as  theory  predicts.  (Theory  predicts  a  gain  of  200  at  400 
W/cm2.)  We  believe  this  is  due  to  "beam  break-up"  resulting  from  spatial  nonuniformities  of  the 
nonlinear  refractive  index  change  in  the  interaction  region.  As  evidence,  Fig.  3  shows 
photographs  of  the  transmitted  pump  beam  with  the  probe  beam  off,  (a)  at  low  intensity  and 
normal  exposure  indicating  the  transmitted  beam  size,  (b)  at  low  intensity  and  long  exposure 
showing  an  "X"  shaped  pattern  of  linear  scattering  that  is  believed  to  be  the  result  of  scattering 
from  surface  and  bulk  defects,  and  (c)  at  high  intensity  and  exposure  equivalent  to  case  (b) 
showing  both  the  nonlinear  scattering  due  to  nonuniformities  in  the  nonlinear  refractive  index 
change,  and  the  linear  scattering  due  to  surface  and  bulk  defects.  Spatial  nonuniformities  of  the 
nonlinear  refractive  index  change  in  the  interaction  region  can  be  the  result  of  shadows  or 
distortions  from  surface  or  bulk  defects,  or  from  spatial  variations  in  the  Cr-ion  concentration. 


(a)  (b)  (c) 

Fig.  3  Photographs  of  the  transmitted  pump  beam,  (a)  Gaussian  beam  size,  (b)  "X"  shaped 
pattern  resulting  from  linear  scattering,  and  (c)  nonlinear  scattering. 

In  conclusion,  we  have  demonstrated  a  gain  of  22  times  by  nondegenerate  two-wave 
mixing  in  Cr:YAlC>3.  The  experimentally  measured  gain  appears  to  be  limited  by  "beam  break¬ 
up"  due  to  spatial  nonuniformities  in  the  nonlinear  refractive  index  change,  that  are  the  result  of 
inhomogeneities  in  the  crystal.  Theoretical  predictions  based  on  our  experimental  data  indicate 
that  gain  in  excess  of  200  times  should  be  possible  in  Cr:YAlC>3  if  a  homogeneous  crystal  can  be 
obtained.  Due  to  its  potential  for  such  high  gain,  Cr:YA103  is  a  very  promising  material  for 
applications  of  nonlinear  optics  in  areas  such  as  signal  processing,  making  it  clear  that  additional 
experimental  efforts  are  warranted  to  explore  nonlinear  response  in  solid-state  Kerr-like  media 
generally,  and  to  investigate  further  the  potential  advantages  of  materials  like  Cr:YA103  over 
photorefractive  crystals. 
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Abstract 

Beam  coupling  gain  of  6  times  was  obtained  by  nondegenerate  two-wave  mixing  in 
Cr:YA103.  To  our  knowledge,  this  is  the  largest  cw  gain  obtained  in  a  solid  state,  non- 
photorefractive  material.  This  material  will  allow  the  study  of  cw  oscillators  and  self-pumped 
conju gators,  with  properties  different  from  those  using  photorefractives. 
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Summary 

Beam  coupling  gain  of  6  times  was  obtained  by  nondegenerate  two-wave  mixing  in 
Cr:YA103.  To  our  knowledge,  this  is  the  largest  cw  gain  obtained  in  a  solid  state,  non- 
photorefractive  material.  This  material  will  allow  the  study  of  cw  oscillators  and  self-pumped 
conju gators,  with  properties  different  from  those  of  photorefractives. 

In  some  metal-ion  doped  insulators,  such  as  ruby,  metastable  states  result  in  low 
saturation  intensities  and  correspondingly  large  optical  nonlinearities.  The  origin  of  the 
nonlinearity  in  these  materials  was  previously  determined  to  be  the  result  of  light  induced 
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population  changes  and  the  difference  in  the  susceptibility  between  the  ground  and  metastable 
states.1  Interest  in  these  materials  stems  from  the  fact  that  they  are  convenient,  solid  state 
materials  with  large  optical  nonlinearities  that  can  be  utilized  with  relatively  low  power  cw 
lasers.  Photorefractive  materials  also  meet  these  conditions  and  have  demonstrated  very  large 
beam  coupling  gain,  but  the  difference  in  the  nonlinear  mechanisms  makes  investigation  of  the 
metal-ion  doped  insulators  of  interest.  As  one  example,  the  beam  ratio  for  which  mutually- 
pumped  phase  conjugation  will  take  place  in  photoreffactives  is  limited,  but  this  is  not  the  case  in 
Kerr  materials.2  Our  search  of  the  literature  on  nonlinearities  in  metal-ion  doped  insulators 
indicated  that  Cr:YA103  stands  out  as  the  material  having  the  largest  potential  index  change  for  a 
given  amount  of  absorption  (i.e.,  largest  value  of  n2ls/oc,  where  Is  is  the  saturation  intensity  and 
n2  and  a  are  the  optical  Kerr  and  absorption  coefficients,  respectively).3  We  therefore  chose  this 
material  for  our  experimental  attempts  to  obtain  high  gain. 

Figure  1  is  a  simplified  schematic  of  the  experiment  used  to  measure  gain  in 
nondegenerate  two-wave  mixing.4  Light  from  an  argon  laser  at  515  nm  is  split  by  a  beam  splitter 
into  a  weak  probe  beam  and  a  strong  pump  beam,  with  the  probe  to  pump  beam  ratio  of  1:100. 
The  frequency  of  the  pump  beam  is  shifted  with  respect  to  the  probe  beam  by  reflecting  it  from  a 
mirror  mounted  to  a  piezoelectric  transducer  and  driven  by  a  triangle  wave  voltage  source.  The 
grating  formed  by  the  interaction  of  the  pump  and  probe  beams  in  the  nonlinear  medium  results 
in  amplification  for  the  probe  when  the  mirror  is  moving  in  one  direction,  and  attenuation  of  the 
probe  when  the  mirror  moves  in  the  opposite  direction.  After  transmission  through  the  nonlinear 
medium,  the  probe  is  detected  and  the  signal  is  recorded  on  an  oscilloscope.  Figure  2  shows  the 
signal  recorded  on  the  oscilloscope  at  a  pump  intensity  of  0.6  kW/cm2.  The  triangle  wave  is  the 
voltage  applied  to  the  piezoelectric  transducer,  the  solid  line  at  the  base  corresponds  to  zero 
signal,  the  dashed  line  corresponds  to  the  signal  level  with  the  pump  off,  and  the  remaining  line 
is  the  signal  with  the  pump  on.  When  the  triangle  wave  is  ramping  down,  there  is  gain  for  the 
probe,  and  at  its  maximum  amplitude  the  gain  is  approximately  6  times.  Figure  3  shows  the  gain 
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as  a  function  of  the  pump  intensity.  The  crosses  are  the  data  and  the  solid  line  is  a  fit  to  theory 
yielding  n2  ~  1 .4  x  107  cm2/W. 

The  6-times  gain  is  nearly  enough  to  demonstrate  a  ring  conjugator.  We  are  currently 
attempting  to  obtain  higher  gain  by  using  materials  with  higher  Cr  doping.  This  research  was 
sponsored  by  the  U.S.  Air  Force  Office  of  Scientific  Research  under  contract  F49620-92-C-0023. 
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Figure  Captions 

Fig.  1  Simplified  schematic  of  the  experiment  used  to  measure  gain  in  nondegenerate  two-wave 
mixing. 

Fig.  2  Signal  recorded  on  the  oscilloscope  at  a  pump  intensity  of  0.6  kW/cm2.  The  triangle 
wave  is  the  voltage  applied  to  the  piezoelectric  transducer,  the  solid  line  at  the  base 
corresponds  to  zero  signal,  the  dashed  line  corresponds  to  the  signal  level  with  the  pump 
off,  and  the  remaining  line  is  the  signal  with  the  pump  on. 

Fig.  3  Gain  as  a  function  of  the  pump  intensity.  The  crosses  are  data  and  the  solid  line  is  a  fit  to 
theory. 
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Abstract 

Two  incoherent  laser  beams  mutually  pump  each  other  to  generate  Stokes-shifted 
phase-conjugate  replicas  of  each  other  via  reflection  gratings  in  an  electrostrictive 
medium.  For  Stokes  scattering  initiated  by  noise,  a  phase-conjugate  reflectivity  of  1%  in 
each  arm  requires  a  combined  Brillouin  gain  of  23.  The  technique  works  for  a  large 
dynamic  range  of  input  beam  ratio  as  the  intensity  of  a  weak  input  beam  is  increased, 
while  the  intensity  of  the  second  input  beam  is  fixed  below  the  threshold  value  for  self- 
stimulated  Brillouin  scattering.  The  weak  beam  is  phase-conjugated  with  reflectivities 
larger  than  unity,  and  the  upper  limit  to  the  dynamic  range  is  set  by  competition  from 
self-stimulated  Brillouin  scattering. 
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1.  INTRODUCTION 

A  mutually  pumped  phase  conjugator  (MPPC)  generates  the  phase-conjugate 
replicas  of  two  incoherent  incident  beams.  Each  input  beam  is  converted  into  the  phase- 
conjugate  replica  of  the  other  by  Bragg  diffraction  off  a  shared  grating.  In 
photorefractive  media,  several  configurations  differing  in  their  number  of  internal 
reflections  from  the  crystal  surfaces  were  demonstrated  using  transmission  gratings.1'5 
MPPC  has  been  used  for  beam  cleanup6  of  a  pulsed  laser  beam,7  self-aligning  injection 
locking  of  a  laser  diode  array,8-9  two-way  optical  communication,10  interferometry,11 
dynamically  programmable  self-aligning  optical  interconnects,12  and  phase-conjugate 
heterodyne  detection.13  For  an  efficient  MPPC,  the  two  incoherent  beams  must  have 
comparable  input  intensities.14  However,  applications  like  phase -conjugate  heterodyne 
detection13  require  the  device  to  work  for  large  imbalance  of  input  beam  intensities,  when 
a  weak  optical  signal  from  a  remote  transmitter  is  combined  with  the  strong  beam  from  a 
local  oscillator.  A  modest  extension  of  the  dynamic  range  can  be  obtained  by  increasing 
the  photorefractive  gain;14  however,  it  is  desirable  to  extend  the  dynamic  range  of  the 
device  to  several  orders  of  magnitude.  We  also  note  that  MPPC  using  reflection  gratings 
in  photorefractive  media  has  been  demonstrated  for  two  configurations:  a  single 
interaction  region  with  no  internal  reflection,15  and  the  geometry  of  two  interconnected 
rings.16 

Saxena  and  Yeh17  have  shown  theoretically  that  MPPC  is  also  possible  in 
electrostrictive  Kerr  media  by  Brillouin  gain  in  a  transmission  geometry.  Their  results 
showed  that  the  Brillouin  intensity  gain  required  for  this  process  was  4,  much  lower  than 
that  required  for  self-stimulated  Brillouin  scattering  (~25)  or  Brillouin  enhanced  four- 
wave  mixing  (~10).18  This  value  is  identical  to  the  threshold  gain  required  for  the  double 
phase-conjugate  mirror  in  photorefractive  media.1*19  An  advantage  of  MPPC  in 
electrostrictive  Kerr  media  over  photorefractive  media  is  the  large  dynamic  range  of  input 
beam  ratio  over  which  the  process  will  occur.17  This  is  because  if  one  beam  is  strong 
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enough  to  satisfy  the  threshold  condition,  then  a  mere  presence  of  the  second  beam  will 
initiate  the  process  of  MPPC.  Ridley20  extended  the  steady-state  analysis  of  Ref.  17  to  a 
time-dependent  one,  and  arrived  at  an  expression  for  the  growth  rate  of  the  forward 
scattered  light  by  ignoring  pump  depletion.  To  the  best  of  our  knowledge,  none  of  this 
work  has  been  corroborated  experimentally.  There  has  been  recent  theoretical  and 
experimental  work  by  Stemklar21  on  MPPC  in  electrostrictive  Kerr  media  using  the 
reflection  geometry.  His  results  showed  that  the  combined  Brillouin  gain  required  for 
this  process  is  of  the  same  order  as  that  for  self- stimulated  Brillouin  scattering  (SBS); 
however,  the  issue  of  large  dynamic  range  for  the  input  beam  ratio  was  not  examined.  In 
this  paper,  we  study  Brillouin  induced  MPPC  in  the  reflection  geometry  and  compare  its 
performance  with  that  in  the  transmission  geometry.  We  show  that  the  device  works  for  a 
large  dynamic  range  of  input  beam  ratio  by  increasing  the  intensity  of  a  weak  input  beam, 
while  fixing  the  intensity  of  the  second  input  beam  below  the  threshold  value  for  self- 
SBS.  The  upper  limit  to  the  dynamic  range  is  set  by  competition  from  self-SBS  of  the 
first  beam.  We  also  note  that  Brillouin  induced  self-pumped  phase  conjugation  of  a 
single  beam  using  reflection  gratings  and  external  mirrors  in  linear  and  loop  cavities  has 
been  studied  in  the  context  of  lowering  threshold  pump  intensity  for  SBS.22 

The  paper  is  organized  as  follows.  In  Section  2,  we  present  the  general  theoretical 
formulation  of  the  problem  and  derive  the  coupled  wave  equations.  In  Section  3A,  we 
obtain  the  threshold  condition  ignoring  pump  depletion.  Analytic  solutions  that  take  into 
account  depletion  of  the  pump  beams  are  obtained  in  section  3B,  with  numerical  plots  of 
the  results  shown  in  section  3C.  Comparisons  with  MPPC  in  transmission  geometry  are 
made  in  section  4,  while  competing  effects  due  to  self-SBS  of  each  pump  are  discussed  in 
section  5.  Finally,  the  results  obtained  in  this  paper  are  summarized  in  Section  6. 

2.  THEORETICAL  FORMULATION 


3 


The  geometry  for  the  nonlinear  interaction  responsible  for  MPPC  is  shown  in  Fig. 
la.  Input  beams  1  and  4  are  a  pair  of  mutually  incoherent  beams  (called  the  "pump" 
beams)  that  enter  through  the  front  face  of  the  nonlinear  medium  and  intersect  at  an  angle 
20.  Acoustic  phonons  that  are  initially  present  in  a  transparent  Kerr  medium  like  CS2  due 
to  thermal  and  quantum  noise  cause  Stokes  scattering  of  the  pump  beams  in  all  directions. 
Also,  a  photon  from  each  pump  beam  can  be  spontaneously  converted  into  a  frequency 
downshifted  Stokes  photon  and  an  acoustic  phonon.  The  coherently  generated  Stokes 
waves  will  interfere  with  the  pump  beam  to  produce  interference  patterns  which  travel  at 
the  acoustic  velocity,  and  which  drive  acoustic  waves  by  electrostriction.23  These 
acoustic  waves  induce  index  gratings  in  the  medium,  and  the  index  grating  that  diffracts 
each  incoherent  input  beam  into  the  phase-conjugated  output  of  the  other  input  beam 
beam  will  be  reinforced  by  both  the  beams.  Hence  this  mutual  index  grating  has  the 
maximum  gain  compared  to  all  the  other  possible  gratings,  and  we  label  the 
corresponding  Stokes  scattered  wave  of  input  beam  1  as  beam  3,  beam  3  also  being  the 
phase-conjugate  of  pump  beam  4.  Similarly,  the  relevant  Stokes  wave  of  pump  beam  4  is 
beam  2,  and  this  is  the  phase-conjugate  of  pump  beam  1.  Steady-state  beam  coupling  is 
described  by  nondegenerate  four-wave  mixing  equations  in  the  nonlinear  medium  with 
appropriate  boundary  conditions. 

We  assume  that  the  various  optical  fields  are  plane  waves,  so  that  the  total  electric 
field  E  in  the  nonlinear  medium  can  be  represented  as 

e  4, 

E(r,t)  =  -f  X  Am(z)  exp[i(km«r  -  comt)]  +  c.c.  ,  (1) 

m=l 

where  em,  Am,  km,  and  com  are  the  unit  polarization  vector,  complex  steady-state 
amplitude,  wave  vector  and  frequency  of  the  m-th  wave.  The  z  axis  is  taken  normal  to 
the  surface  of  the  medium,  and  the  complex  amplitudes  are  assumed  to  be  functions  of  z 
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due  to  energy  transfer  by  beam  coupling.  The  intensity  pattern  generated  by  the 
interference  of  the  coherent  beams  1  and  3,  and  the  coherent  beams  2  and  4,  is  given  by 


* 

<E2(r,t)>  =  2  {X'Aml2+  [(e1.e3*)A1A*3exp[i(K13.r-a)i3t)] 

m=l 


+  (e4*e2*)  A4  A*  exp[i(K42*r  -  co42t)]  +  c.c.]  }  ,  (2) 


where  the  angular  brackets  indicate  averaging  over  times  long  compared  to  the  optical 
periods,  Kij  =  kj  —  kj,  and  coy  =  <0£  —  C0j.  Each  interference  term  induces  a  volume  index 
grating  in  the  medium  via  the  electrostrictive  Kerr  effect,  given  by23*24 


Any  = 


(e^e,*)  KjjV 


8nope0  (v  K;j  -  CDij 2  -  i  toy  Ty) 


AiAjexp[i(Kyr  -  a>ijt)]  +  c.c. 


(3) 


where  y  is  the  electrostrictive  coefficient,  n0  is  the  linear  refractive  index  of  medium  in 
absence  of  light,  p  is  the  mass  density,  Eq  is  the  dielectric  constant  of  vacuum,  v  is  the 
acoustic  speed  and  Ty  is  the  Brillouin  linewidth  that  depends  on  experimental  geometry. 
We  wish  to  consider  the  case  of  equal  frequencies  and  wavevectors  for  the  two  index 
gratings  induced  by  the  intensity  pattern  to  yield  a  common  or  shared  grating:  coi3  =  co42 
=  Aco,  Ki3  =  K42  =  K,  and  Ti3  =  r42  =  T.  Then  the  total  light  induced  change  in  the 
refractive  index  of  the  medium  is  given  by 


An  =  An, 3  +  An42 -  *¥  Kgl^)  A,Aj  +(e4^«)  A^j 

8n0pe0  (v2K2  -  Aco2  -  i  Aco  T) 


exp[i(K-r  -  Acot)] 


+  c.c. 


(4) 
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To  describe  the  propagation  and  coupling  of  the  waves  in  the  nonlinear  medium  at 
steady-state,  we  start  with  the  scalar  wave  equation 


v2e+2^e  =  o 


(5) 


where  n  =  (no  +  An)  is  the  total  refractive  index  of  the  medium.  Since  Aco  « co,  we  have 
ignored  the  small  frequency  differences  between  the  waves  and  taken  com  »  co  for  all  the 
waves.  We  have  also  neglected  the  small  absorption  of  the  waves  in  the  transparent  Ken- 
medium.  Substituting  for  n  and  E  in  the  wave  equation,  using  the  slowly  varying 
envelope  approximation  and  ignoring  terms  that  are  quadratic  in  the  index  change,  we 
arrive  at  the  following  coupled  wave  equations 


IE1  =  2  Kel*e3*)  +  (e4*e2*)  ^2A4]  A3 

=  \  [(et**e3)  ^1A3  +  (e4**e2)  A2^*4  ]  A4 

-fa  =  -§~  [(ei*#e3)  ^1A3  +  (e4**e2)  A^4  ]  Ax 

=  2  ^el*e3*)  A1^3  +  (e4*e2*)  ^2A4]  A2  •  (6) 

Here  g  is  the  complex  Brillouin  gain  coefficient  that  is  given  by 

ito  K-V 

g  -  - - — - — *- - — - -  (7) 

ccos9  4nopeo(v2K2- Am2-iAcor) 

Note  that  g  has  units  of  inverse  of  the  product  of  intensity-length.  The  complex  gain 
coefficient  is  a  negative  real  number  when  the  frequency  difference  between  two  coherent 
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waves  is  equal  to  the  acoustic  frequency:  coi  -  (03  =  (04  -  0)2  =  Aco  =  vK,  indicating  a  90° 
phase  shift  between  the  index  grating  and  the  intensity  pattern.23,24  In  Stemklar's 
experiment,21  the  two  input  beams  were  orthogonally  polarized  with  respect  to  one 
another  to  enable  tagging  of  the  Stokes  photons  of  each  pump  by  virtue  of  their 
polarization.  Let  input  beam  1  be  linearly  polarized  in  the  plane  of  incidence  (p- 
polarization),  and  let  input  beam  4  be  linearly  polarized  perpendicular  to  the  plane  of 
incidence  (s-polarization).  Since  each  coherently  generated  Stokes  wave  has  the  same 
polarization  as  its  pump  wave,  =  cos20,  e4»e2  =  1.  The  angle  20  between  the 
incoherent  beams  is  typically  small  (~10),21  so  that  the  dot  product  between  the  unit  p- 
vectors  may  also  be  taken  equal  to  unity.  At  resonance  (Aco  =  vK),  the  coupled  wave 
equations  (6)  reduce  to 

~6z  =  “  ¥  *  A>^3  +  A*4  *  A3 

^  =  (A>j+A2<)A4 

=  (A*A3+A2<,Ai 

IT  -  -f  (a,4+^a4)a2 

where 

gR  s  mKR t 

4cnoeovprRcos0 

and  the  subscript  R  indicates  the  shared  reflection  grating  under  consideration,  and  Kr  = 
2k  cos0.  The  electrostrictive  response  of  the  medium  has  been  stated  in  general  terms, 
and  can  be  easily  changed  from  reflection  to  transmission  geometry17  by  changing  the 
interference  terms.  From  the  coupled  wave  eqs.  (8),  we  note  that  that  there  is  no  phase 


mismatch  associated  with  the  four-wave  mixing  term  along  the  z-direction  for  equal 
frequencies  of  the  two  incoherent  pump  beams.  We  now  solve  the  coupled  wave 
equations  (8)  with  and  without  pump  depletion. 


3A.  Threshold  Condition  using  Undepleted  Pump  Approximation 

If  the  two  input  pump  beams  are  much  stronger  than  the  self-generated  Stokes  beams  (Aj, 
A4  »  A2,  A3),  then  we  can  take  dAj/dz  =  dA4/dz  =  0,  so  that  A  j ,  A4  are  constant  for  all 
z,  i.e.,  Aj(z)  =  Aj(0)  s  A10,  and  A4(z)  =  A4(0)  =  A40.  Hence  the  coupled- wave  equations 
(8)  reduce  to 

dA9  jic 

dz  =  2"  (  A10  A4Q  A3  + 140  A2) 

dAi 

=  2~  ^  ^10  A3  +  A40  A2)  »  (10) 


where  Ij  =  IAi|2.  The  two  equations  can  be  easily  decoupled  and  solved  to  give 


A2(z)  = 


(to  I40  exp(- gRInz/2) 

Iq  a2(^1 


A10  A40  [l  -  exp(-  gRlpz/2)] 

In  A3(0) 


A3(z)  = 


I4o  +  I]0exp(-gRInz/2) 

I0  a3  ((') 


AiqA40[1  -  exp(-  gRlnz/2)] 
1q 


A2(0)  ,  (11) 


where  Io  =  I10+I40  is  the  sum  of  the  input  pump  intensities.  If  A2(L)  =  A2L  and  A3(L)  = 
A3L  represent  the  small,  nonzero  seed  values  of  the  Stokes  waves  at  the  z  =  L  plane, 

where  L  is  the  interaction  length,  then  the  expressions  for  the  Stokes  waves  at  the  z  =  0 
plane  simplify  to 


A2(0)  = 


I10  +  I4o  e°-5G 


V2L 


A*o  a4o  ( eQ5G  -  1) 


L3L 
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a3(0)  = 


A3L  + 


(12) 


A1QA4o(e°  5G  ~  1) 

Io 

where  G  =  gRIo  L  is  the  total  intensity  Brillouin  gain  for  the  two  incoherent  pumps.  These 
simple  solutions  contain  interesting  information  about  MPPC  using  reflection  gratings. 
There  will  not  be  any  Stokes  waves  at  the  z  =  0  plane  if  both  the  seeds  are  zero  at  the 
back  plane:  A2l  =  A3L  =  0.  However,  only  one  nonzero  seed  is  sufficient  to  generate 
both  the  Stokes  waves  at  the  output.  As  discussed  by  Tikhonchuk  and  Zozulya25  for 
similar  configurations,  the  MPPC  process  under  consideration  in  this  paper  is  not  an 
oscillator  but  an  amplifier  of  small  input  seed  values  of  the  scattered  radiation.  The 
reflection  geometry  MPPC  has  a  hard  excitation  threshold  in  that  nonzero  seeds  are 
required  for  operation.  Notice  also  that  each  Stokes  wave  has  contributions  from  two 
terms.  The  first  term  is  due  to  the  exponential  buildup  of  the  seed  in  the  presence  of  two 
pumps,  one  coherent  and  the  other  incoherent  with  respect  to  the  Stokes  wave.  The 
incoherent  pump  also  contributes  to  the  Brillouin  gain  due  to  cross  readout  of  the  mutual 
reflection  grating.  The  second  term  is  the  4WM  term  that  leads  to  phase  conjugation.  If 
one  of  the  pump  beams  is  zero  along  with  its  coherently  generated  Stokes  wave,  then  the 
MPPC  solutions  of  Eqs.  (12)  reduce  to  the  familiar  exponential  buildup  of  the  remaining 
Stokes  wave.  The  finite  Stokes  seeds  at  the  z  =  L  backplane  arise  from  the  mutual 
scattering  of  each  pump  in  a  direction  counterpropagating  to  the  other  pump: 


A2L  —  (  A4(L)  »  A3L  —  f  Aj(L) 


(13) 


where  f  is  the  mutual  scattering  coefficient,  assumed  equal  in  amplitude  and  phase  for 
both  the  pumps.  A  realistic  picture  should  describe  seeding  by  distributed  noise  terms, 
yet  we  use  the  simplified  approach  of  self-SBS  in  assuming  deterministic  seed  levels  at  a 
point  boundary.  Since  the  pumps  are  assumed  to  be  constant  throughout  the  interaction 
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region,  we  have  A2L  =  f  A40,  A3L  =  f  A10.  Substituting  boundary  conditions  (13)  into  the 
solutions  (12)  and  simplifying,  we  obtain 

I2(0)  =  Ifl2  I40  eG  f  I3(0)  =  Ifl2  I10  eG  .  (14) 


In  SBS,  a  Stokes  wave  originates  from  spontaneous  scattering  of  the  pump  beam  by  index 
gratings  induced  by  the  acoustic  noise  that  is  always  present  in  a  medium  at  thermal 
equilibrium.  Since  the  spontaneous  scattering  coefficient  is  very  small,  and  the  Brillouin 
gain  is  large  only  in  a  narrow  solid  angle  for  a  focused  pump  beam,  the  intensity  of  the 
Stokes  seed  is  taken  to  be  10-11  to  1013  of  the  pump  intensity.26  In  our  calculations,  we 
take  Ifl2  =  1012.  If  we  define  threshold  as  the  condition  when  1%  of  each  pump  is 
converted  to  its  Stokes  wave  at  the  output,  i.e.,  l2(0)/l4o  =  I3(0)/Ii0  =  0.01,  then  the 
Brillouin  gain  at  threshold  is  given  by 

°th  =  gR(Iio+l4o)  L  =  10  lne10  =  23.2  ,  (15) 


so  that  the  sum  of  the  input  intensities  of  both  the  pumps  have  to  yield  a  Brillouin  gain  of 
23  at  threshold.  This  number  is  comparable  to  Brillouin  gain  required  for  SBS  initiated 
by  noise.  Hence  the  intensity  of  each  pump  can  be  below  the  threshold  value  required  for 
self-SBS,  and  yet  MPPC  will  build  up  from  noise  to  1%  phase-conjugate  reflectivities  if 
the  sum  of  the  Brillouin  gain  for  the  two  pumps  satisfies  Eq.  (15). 

3B.  Analytic  solutions  with  Pump  Depletion 

The  boundary  conditions  appropriate  for  this  problem  are  the  field  strengths  A10,  A40  of 
the  two  pumps  at  the  z  =  0  input  plane,  and  the  finite  seed  values  of  the  Stokes  waves 
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A2l>  A3L  at  the  z  =  L  plane.  As  the  beams  propagate  in  the  nonlinear  medium,  the 
following  quantities  are  conserved:27-28 


Ii(z)-l3(z)  =  Ci 

^(z)  -  I2(z)  =  c2 

Ai(z)A2(z)  -  A3(z)A4(z)  =  di 
A2(z>^3(z)  -  X’Ci(z)A4(z)  =  d2 


(16) 


where  Ci,  c2,  di,  and  d2  are  integration  constants  that  are  yet  to  be  determined.  With  the 
help  of  these  conservation  laws,  the  four  coupled  equations  of  Eqs.  (8)  can  be  decoupled 
to  yield  the  following  two  first-order  differential  equations  for  the  new  variables  A2j  = 
A2/A*  and  A34  =  A^A^ : 


dA 


21 


dz 


y  [di*A21  +  N  A2i  —  di  ] 


dA 


dz 


M  =  -  f-d,*A2 


2  di  A^  +  N  A34  +  di  ] 


(17) 


where  N  —  ci+c2.  Equations  (17)  can  be  readily  solved  to  yield  the  following  solutions 


A2i(z)  = 


QA21(0)  +  tanh  [2di  -  N  A21(0)] 

Q  +  tanh  (^^)  [N  +  2  di*A21(0)] 


A34(z)  = 


QA34(0)  -  tanh  (^^)  [2di  +  N  A^O)] 
Q  -  tanh  (^)  [2  di*A34(0)-N] 


(18) 


where  the  parameter  Q  is  defined  as 


Q  =  ( N2  +  4  Id!  I2  )1/2 


(19) 


In  obtaining  the  solutions  for  A21(z)  and  A34(z)  given  by  Eqs.  (18),  no  approximations 
have  been  made  regarding  di  and  N.  If  we  define  the  phase-conjugate  reflection 
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coefficients  as:  p  s  A3(0)  /A*40  =  A34(0)  and  a  =  A2(0)  /A%  =  A21(0),  then  from  the  above 
equations,  we  obtain 


p  =  QA34(L)  +  Tl  [N  A34(L)  +  2dfi 
Q  +  TL  [2  di  *A34(L)  -  N] 

a  =  QA2i(L)  +  TL[NA21(L)-2d!] 

Q  -  TL  [N  +  2  di*A21(L)]  ’  (  } 

where  Tl  =  tanh  (gRQL/4).  We  now  assume  the  same  form  of  boundary  conditions  as 
given  by  Eqs.  (13)  of  the  previous  section:  A2L  =  f  A4(L),  A3L  =  f  Aj(L),  i.e.,  the  mutual 
scattering  is  equal  in  amplitude  and  phase  for  both  the  pumps  at  the  end  face  of  the 
nonlinear  medium.  Then  the  first  and  second  conservation  laws  yield  the  pump 
intensities  at  the  end  face 


Ii(L) 


Iio  -  R I40 
1-  |fP 


U(L) 


I40  -  S  lip 
l-|fl2 


(21) 


where  R  and  S  are  the  power  phase-conjugate  reflectivities:  R  =  |p|2,  S  =  |a|2.  With  the 
boundary  conditions  given  by  Eqs.  (13),  the  third  conservation  law  shows  that  the 
integration  constant  di  =  0.  Consequendy 

R  - 

S  =  q2  ,  (22) 


where  q  is  defined  as  the  intensity  ratio  of  the  pump  beams  at  the  input:  q  =  I10  / 140.  This 
relationship  is  also  true  for  MPPCs  in  photorefractive  media.1'5-28  From  the  definition  of 
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phase-conjugate  reflectivities  R  and  S,  it  is  obvious  that  their  maximum  values  occur 
when  all  of  the  pump  beams  are  converted  into  the  Stokes  beams.  Hence  Rm  =  Ik/Uo  = 
q,  and  Sm  =  I40/I10  =  1/q,  where  the  subscript  m  denotes  maximum  value.  Another 
relevant  unknown  for  MPPC  is  the  device  transmissivity  in  each  direction:  Ty  =  Ii(0)/Ij0, 
and  it  reveals  the  fraction  of  each  input  beam  that  is  converted  into  the  phase-conjugate 
output  of  the  other  beam.  Using  Eq.  (22),  it  can  be  shown  that  T31  =  T24  a  T,  so  that  both 
input  beams  have  equal  transmissivity.21  This  relationship  is  also  true  for  MPPCs  in 
photorefractive  media. 1  The  maximum  value  of  the  device  transmissivity  is  unity,  when 
all  the  energy  of  an  input  beam  is  converted  into  the  phase-conjugate  of  the  other  beam. 
Since  di  =  0,  Q  =  N  and  Tl  =  tanh  (gRNL/4).  Hence  Eqs.  (20)  simplify  to 

R  =  exp(gRNL)  IA34(L)|2  ,  s  =  exp(gRNL)  IA21(L)I2  .  (23) 

Now  N  =  I10  + 140  - 12(0)  - 13 (0)  is  still  an  unknown  quantity  and  may  be  expressed  as  N 
=  Io(l-R/q).  Using  Eqs.  (21),  we  may  also  show  that  IA34(L)I2  =  q  Ifl2.  Substituting  N 
and  IA34(L)I2  in  the  above  equation  for  R,  we  obtain  the  following  transcendental 
equation  for  R: 


R  =  q  Ifl2  exp[G(l-R  /  q)]  .  (24) 

Note  that  if  Ifl2  =  0,  i.e.,  I2L  =  I3L  =  0  in  eqs.  (13),  then  the  phase-conjugate  reflectivities  R 
and  S  are  equal  to  zero.  Hence  finite  seed  values  of  the  Stokes  waves  are  required  at  the 
back  surface  of  the  nonlinear  medium  for  MPPC  to  work  in  reflection  geometry.  This  is 
in  contrast  to  MPPC  using  transmission  gratings,17  where  steady-state  phase-conjugate 
reflectivities  can  be  obtained  with  zero  Stokes  seeds  at  either  edge  of  the  nonlinear 
medium.  Equation  (24)  can  be  solved  numerically  for  R  using  known  values  for  the 
combined  Brillouin  intensity  gain  G,  intensity  ratio  q  for  input  beams,  and  the  mutual 
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scattering  coefficients  of  the  pumps  equal  to  Ifl2.  Knowing  R,  we  can  calculate  S  using 
Eq.  (22):  S  =  R/q2. 

3C.  Numerical  Results 

Figs.  2a, b  are  plots  of  the  phase-conjugate  reflectivities  R,  S  and  transmissivity  T  as  a 
function  of  Brillouin  gain  G4  =  gR  I40  L  for  Ifl2  =  1012,  and  different  values  of  beam  ratio 
q:  (a)  q  =  1,  (b)  q  =  0.5.  Note  that  the  input  beam  ratio  q  =  1 10  / 140  affects  the  total 
Brillouin  gain:  G  =  gR(Iio+l4o)  L  =  G4(l+q)  as  we  keep  I40  and  hence  G4  constant.  For 
equal  input  intensities  of  beam  1  and  4,  the  two  phase-conjugate  reflectivities  and 
transmissivity  are  all  equal  and  approach  the  limiting  value  of  unity  for  large  coupling 
gain.  The  threshold  Brillouin  gain  for  beam  4  is  around  11.5,  and  since  q  =  1,  this 
corresponds  to  a  total  Brillouin  gain  of  23  at  threshold,  as  shown  in  section  3A.  If  pump 
beam  4  is  twice  as  intense  as  pump  beam  1  at  the  input  (i.e.  q  =  0.5,  as  in  Fig.  2b),  then 
the  weaker  pump  has  a  larger  phase-conjugate  reflectivity  (=  S)  that  approaches  the 
limiting  value  of  2  (  =  1/q)  for  large  Brillouin  gain.  The  stronger  pump  has  a  smaller 
phase-conjugate  reflectivity  (=  R)  that  approaches  the  value  0.5  (=  q)  for  large  Brillouin 
gain.  Transmissivity  T  (shown  as  a  dashed  curve)  approaches  the  maximum  value  of 
unity  for  large  gain.  Since  pump  beam  1  is  now  weak,  G4  at  threshold  has  to  be  larger  in 
Fig.  2b  than  the  threshold  value  in  Fig.  2a  so  that  condition  given  by  Eq.  (15)  is  still 
satisfied  for  the  total  gain. 

Figs.  3a,b  are  a  plot  of  the  phase-conjugate  reflectivities  R,  S  and  transmissivity  T 
as  a  function  of  the  input  intensity  ratio  q  of  the  two  pumps  on  a  log-log  scale,  for 
different  values  of  Brillouin  gain  G4  satisfied  by  pump  beam  4.  The  values  of  G4  are 
chosen  below  the  threshold  value  required  for  self-SBS  (  ~  23)  of  pump  beam  4.  The 
upper  limit  for  q  value  is  set  by  self-SBS  of  beam  1;  since  the  coupled  wave  eqs.  (8)  do 
not  include  self-SBS  terms  for  each  pump,  our  formulation  is  not  valid  when  competing 
effects  such  as  self-SBS  set  in  when  each  pump  has  a  Brillouin  gain  >  23.  However,  in 
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section  5,  we  show  that  the  gain  for  phase-conjugation  by  MPPC  is  twice  as  large  as 
phase-conjugation  by  self-SBS,  so  that  the  Brillouin  gain  would  have  to  be  very  large 
before  self-SBS  begins  to  compete  with  MPPC.  Increasing  q  is  equivalent  to  increasing 
Iio,  keeping  I40  and  hence  G4  fixed.  For  G4  =  12  (Fig.  3a)  and  small  q  values,  (10*8  to  10' 
2)  the  total  Brillouin  gain  is  roughly  half  that  required  for  threshold.  Consequently,  very 
little  energy  of  the  pump  beams  is  converted  into  the  Stokes  waves,  and  since  S  is 
inversely  proportional  to  Iio,  S  decreases  with  increasing  q.  As  q  is  further  increased,  the 
combined  Brillouin  gain  approaches  threshold,  leading  to  increased  conversion  of  pump 
beam  I40  into  the  Stokes  beam  12(0),  as  seen  by  the  increase  in  T.  This  compensates  for 
the  decrease  of  S  with  increasing  Iio,  so  that  S  actually  increases  for  a  limited  range  of  q 
values.  When  the  total  Brillouin  gain  is  enough  to  convert  most  of  the  pumps  into  their 
Stokes  waves,  then  any  further  increase  in  q  results  in  decreasing  S  once  again.  R  is  a 
monotonically  increasing  function  of  q  because  its  denominator  I40  is  a  fixed  quantity, 
and  increase  in  Iio  results  in  more  conversion  of  this  pump  beam  into  its  Stokes  wave 
13(0).  Fig.  3b  is  a  similar  plot  for  G4  =  18.  Notice  that  the  variation  between  the 
maximum  and  minimum  values  of  S  is  nearly  3  orders  of  magnitude  smaller  in  Fig.  3b 
than  in  Fig.  3a  due  to  the  larger  Brillouin  gain  G4  to  begin  with  in  Fig.  3b.  Even  when 
the  Brillouin  gain  of  each  beam  is  smaller  than  the  threshold  gain  for  self-SBS,  figures  3 
reveal  that  the  weak  beam  can  be  phase-conjugated  with  a  reflectivity  greater  than  1, 
opening  the  possibilities  for  a  variety  of  applications  in  phase-conjugate  resonators, 
phase-conjugate  interferometers  and  phase-conjugate  gyros.  When  G4  =  12,  S  >  1  for  q  < 
10'7  (fig.  3a),  and  when  G4  =  18,  S  >  1  for  q  <  104. 

4.  Comparison  with  MPPC  in  transmission  geometry 

To  facilitate  this  comparison,  figures  4a  and  4b  show  MPPC  in  electrostrictive  media 
using  transmission  geometry.  The  corresponding  coupled-wave  equations  from  Ref.  17 
are: 
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^A-l  §T  /•  *  .  **  a  iAkz\  » 

dz  “  2  (  ^1^4  +  A2A3  e  )  A4 


^  =  f(A*A4e^  +  A24)A3 
^  =  -f(AAe-“*  +  A*2A3)A2 
^  =  -f  (A'A.  +  A^'^A, 


where 

coKtY2 

gT  = - 11 -  (26) 

4cnoeovprTcos0 

For  transmission  geometry  denoted  by  subscript  T,  20  is  the  angle  between  coherent 
beam  pairs  1,  4,  and  2,  3,  so  that  Kt  =  2k  sin0.  To  compare  the  Brillouin  gain 
coefficients  g,  and  phonon  lifetimes  X;  for  the  various  geometries,  we  must  express  the 
geometrical  dependence  contained  in  the  acoustic  wavevector  magnitude  Kj  and  Brillouin 
linewidth  Tj  explicitly.  We  use  the  relationship24’29:  Tj  =  if1  =  riKi2/p,  T|  being  the 
viscosity  of  the  medium.  This  equation  is  true  only  for  acoustic  frequencies  larger  than 
the  vibrational  relaxation  frequency  of  the  electrostrictive  medium,  when  acoustic 
absorption  can  be  described  by  classical  mechanisms  only.30  Contradirectional  SBS 
corresponds  to  0  =  0  in  reflection  geometry  (the  angle  between  a  pair  of  coherent  beams 
is  rc-20),  and  gsBS  =y2'/Ser\v,  Tsbs  =  p/4rik2,  £  being  the  uniform  dielectric  constant  of  the 
medium.  Hence  gR/gsBS  =  l/cos2©  =  Xr/xsbs.  so  that  a  2k  SBS  grating  has  a  smaller  gain 
than  the  reflection  grating,  and  acoustic  phonons  corresponding  to  contradirectional  SBS 
have  the  shortest  lifetime.  Similarly,  gT  =  2gsBs/sin20,  and  xt/xSbs  =  l/sin2©.24*29  Hence 
the  Brillouin  gain  is  minimum  for  contradirectional  SBS,  and  increases  with  smaller 
beam  angles  in  the  forward  direction.  Yet  the  acoustic  phonons  are  difficult  to  build  up 
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in  forward  MPPC  as  they  tend  to  walk  off  the  interaction  region.  Also,  the  lifetime  of  the 
acoustic  phonons  is  much  longer  in  the  forward  direction,  so  that  MPPC  in  transmission 
geometry  requires  not  only  intense  pump  beams  but  also  large  pulsewidths  to  achieve 
steady-state.  This  may  explain  why  MPPC  in  electrostrictive  media  has  not  been 
observed  in  the  transmission  geometry.  Notice  the  phase  mismatch  Ak  =  (2  n  Am  /  c) 
associated  with  the  transmission  geometry  that  is  absent  in  the  reflection  geometry 
(compare  with  Eqs.  (8)  of  section  3A).  Also  note  that  the  phase-conjugate  reflectivities  in 
Ref.  16  were  obtained  with  the  boundary  conditions  of  zero  seeds  at  their  startup  planes: 
Iio  =  I3L  =  0.  In  reflection  geometry,  non-zero  seeds:  I2L  *  0, 13L  *  0  have  to  be  assumed 
at  the  end  face  of  the  nonlinear  medium  to  obtain  finite  phase-conjugate  reflectivities. 
Hence  MPPC  in  reflection  geometry  has  a  hard  excitation  threshold.  Another  related 
difference  in  the  two  formulations  is  the  manner  in  which  the  threshold  condition  is 
derived,  along  with  the  value  of  the  Brillouin  gain  required  for  threshold.  For  equal 
intensity  incoherent  beams,  the  threshold  condition  in  Ref.  16  is  derived  from  the 
transcendental  equation:  gT(l40+l2L)L  =  4.  No  physical  solutions  will  be  obtained  below 
this  threshold  value  of  intensity  Brilloiun  gain,  and  the  phase-conjugate  reflectivities  are 
92%  at  threshold  gain.  In  this  paper,  the  threshold  condition  is  defined  as  the  Brillouin 
gain  required  to  give  1%  phase-conjugate  reflectivities.  For  Stokes  waves  initiated  by 
noise  (Ifl2  =  10'12),  this  requires  a  combined  Brillouin  gain  of  gR(Iio+l4o)L  =  23. 

5.  Competition  due  to  self-SB  S 

The  coupled  wave  Eqs.  (8)  are  no  longer  valid  if  one  or  both  the  pump  beams  is  intense 
enough  to  satisfy  the  threshold  for  self-SBS.  Two  additional  waves  must  be  included  in 
the  formalism  to  describe  the  Stokes  scattering  of  each  pump  beam  due  to  backward  SBS. 
Let  E5  (Eg)  be  the  coherently  generated  Stokes  wave  of  pump  beam  1  (pump  beam  4)  due 
to  2k  gratings  (see  Fig.  5).  Wave  E5  (Eg)  is  collinear  with  beam  2  (beam  3),  but  it  has  a 
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frequency  equal  to  co  —  AcOsbs*  where  Acdsbs  is  the  SBS  acoustic  frequency  equal  to  2kv. 

Hence  the  modified  coupled  wave  Eqs.  (8)  with  terms  due  to  self-SBS  are 

dAi 


■d^=-T  (A1^3  +  A*2A4)A3-^  AjIAsI2 


^  =  -f  (A*A3+A2<)A4 

^  =  ~^2  (^lA3  +  A2^4  )A1 


dA4 

dz 

dAj 

dz 

dA^ 

dz 


=  -^  (A1Ar3  +  A*2A4)A2  -pp  A4IA6|2 


gSBS 


-  AjIA^2 

-  A6IA4I2 


(27) 


where  gsBS  has  been  defined  in  Section  3.1.  Note  that  since  both  Stokes  waves  3  and  5 
are  coherently  generated  from  pump  beam  1,  they  could,  in  principle,  write  a  transmission 
grating  between  themselves.  Such  secondary  gratings  written  among  the  various  noise 
terms  would  be  very  weak  compared  to  the  primary  gratings  written  by  the  pumps  and 
noise,  except  near  the  input  face  of  the  nonlinear  medium  where  the  Stokes  waves 
undergo  considerable  amplification.  In  what  follows,  we  will  ignore  such  gratings. 
Similar  remarks  apply  for  the  Stokes  waves  2  and  6  that  are  coherently  generated  from 
pump  beam  4. 


If  we  assume  that  both  the  pumps  are  much  stronger  than  their  Stokes  waves,  we 
can  ignore  pump  depletion  and  consider  the  input  pump  amplitudes  as  constant  quantities 
throughout  the  interaction  region.  In  that  case,  the  above  equations  reduce  to 

dA2  £p  * 
dz  =  ~  ~2  (  Aio  A40  A3  + 140  A2) 


18 


^  =  -f  (IioA3  +  A10A'J0A2) 


_  _  gSBS 
dz  "  2 

^6  _  _  gSBS 
dz  -  2 


A5  IlO 


A6  140 


(28) 


In  this  limit,  the  coupled  wave  equations  for  beams  2  and  3  are  still  independent  of  beams 
5  and  6.  With  A10  and  A40  as  constant  quantities,  the  solutions  given  by  Eqs.  (12)  for 
A2(0)  and  A3(0)  are  still  valid,  despite  the  presence  of  additional  Stokes  wave  5  and  6. 
Hence  self-SBS  does  not  modify  the  solutions  for  MPPC  for  negligible  pump  depletion. 

The  solutions  for  waves  2  and  3  will  be  affected  by  self-SBS  of  the  pumps  only  if 
we  take  into  account  pump  depletion.  To  do  this,  one  must  solve  the  six  nonlinear 
coupled  wave  equations  (26)  numerically,  with  two-point  boundary  values.  However,  we 
will  now  show  qualitatively  that  the  Brillouin  gain  coefficient  for  the  mutual  reflection 
grating  can  be  twice  as  large  as  the  2k  SBS  gratings,  so  that  the  latter  may  be 
discriminated  against  owing  to  a  lower  gain.  We  use  the  formulation  developed  by  Yeh31 
for  MPPC  in  photoreffactive  media  to  arrive  at  this  qualitative  result.  We  have  seen  in 
sections  3.2  and  3.3  that  the  two  phase-conjugate  reflectivities  are  comparable  for  similar 
intensity  of  the  two  input  beams.  Let  the  amplitudes  of  the  various  Stokes  waves  be 
expressed  as  the  phase-conjugates  of  the  relevant  pumps: 

A2  =  PMPPC  A*1  ,  A3  =  pMPPCA4  >  A5  =  PsbsA*  ,  A6  =  pSBSA4  (29) 

where  pR  and  Psbs  are  arbitrary  constants.  Substitution  of  Eqs.  (29)  into  (27)  results  in 
the  following  set  of  coupled-wave  equations: 

=  -  ( gR  'Pmppc  I2  IA4I2  +  IpsBsl2  IAjI2  )  Aj 
=  -  gR  Pmppc  IA4I2 
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dA3 


=  “gR  Pmppc  A4IAJ2 


dA4 


dT  =  “  (  8r  1  Pmppc  I2  lA^2  +  yP  IpsBsl2  IA4I2 )  A4 


^=-IfSPSBSA*1  IAJ2 

_6  _  _  1§M  psBs  ^  |a4|2 


(30) 


Hence  if  the  p's  and  g's  have  comparable  magnitudes,  then  the  above  equations  show  that 
phase  conjugation  due  to  MPPC  builds  at  twice  the  gain  as  SBS  of  each  pump  beam 
because  of  the  reinforcement  of  the  mutual  reflection  grating  in  MPPC  from  both  pairs  of 
coherent  beams.  This  factor  of  2  in  the  exponential  gain  is  similar  to  the  gain  for  the 
phase-conjugate  wave  due  to  SBS  being  twice  as  large  as  the  gain  for  any  other  scattered 
wave,  so  that  the  scattered  waves  are  discriminated  against  owing  to  their  lower  gain.26 


6.  Conclusions 

We  studied  Brillouin  induced  MPPC  in  the  reflection  geometry,  and  showed  that  a  phase- 
conjugate  reflectivity  of  1%  in  each  arm  requires  a  combined  Brillouin  gain  of  23  for 
Stokes  scattering  initiated  by  noise.  MPPC  works  for  a  large  dynamic  range  of  input 
beam  ratio  by  increasing  the  intensity  of  a  weak  input  beam,  and  fixing  the  intensity  of 
the  second  input  beam  below  the  threshold  value  for  self-SB  S.  Our  numerical  plots  show 
that  the  weak  beam  is  phase-conjugated  with  reflectivities  larger  than  unity,  and  the  upper 
limit  to  the  dynamic  range  is  set  by  competition  from  self-SBS  of  the  first  beam.  We 
have  also  compared  the  performance  of  MPPC  with  reflection  gratings  to  that  using 
transmission  gratings,  and  commented  on  why  the  latter  may  be  difficult  to  observe 
experimentally.  We  have  shown  that  the  performance  of  MPPC  is  not  affected  by  self- 
SBS  for  negligible  pump  depletion  near  threshold.  Even  under  conditions  where  pump 
depletion  takes  place,  we  have  presented  an  argument  indicating  that  the  gain  for  phase 
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conjugation  by  MPPC  is  twice  that  for  self-SBS  due  to  the  mutual  reflection  grating 
formed  in  MPPC. 
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FIGURE  CAPTIONS 


Fig.  la 

Fig.  lb 
Figs.  2 

Figs.  3 

Fig.  4a 

Fig.  4b 
Fig.  5 


Schematic  diagram  showing  MPPC  using  reflection  gratings  in  an 
electrostrictive  Kerr  medium. 

Wave  vector  diagram  for  MPPC  using  reflection  gratings. 

Phase-conjugate  reflectivities  R,  S  and  transmissivity  T  as  a  function  of 
Brillouin  gain  G4  for  Ifl2  =  10'12  and  (a)  q  =  1  (b)  q  =  0.5. 

Phase-conjugate  reflectivities  R,  S  and  transmissivity  T  as  a  function  of 
intensity  ratio  of  the  two  input  beams  q  for  Ifl2  =  1012  and  (a)  G4  =  12  (b) 
G4  =  18. 

Schematic  diagram  showing  MPPC  using  transmission  gratings  in  an 
electrostrictive  Kerr  medium. 

Wave  vector  diagram  for  MPPC  using  transmission  gratings. 

Schematic  diagram  showing  MPPC  reflection  gratings  and  2k  SBS 
gratings. 
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2  Nonlinear  Optical  Phase  Conjugation  by  Four-Wave 
Mixing  in  Kerr  Media 

R.  Saxena  and  P.  Yeh1 


2.1  Introduction 


Optical  phase  conjugation  is  a  technique  for  reversing  both  the  direction  of 
propagation  and  the  overall  phase  factor  of  an  incoming  light  wave  [2.1-5]. 
This  retroreflecting  property  of  phase-conjugate  waves  is  useful  in  applications 
requiring  automatic  pointing  and  tracking  [2.1-5],  phase  aberration  correc¬ 
tions  [2.1-5],  phase-conjugate  resonators  [2.6-8],  interferometers  [2.9-15], 
gyros  [2.16-20],  and  many  other  devices. 

Various  nonlinear  optical  effects  can  be  utilized  to  conjugate  the  phase  of 
an  optical  beam  in  real  time.  Nonlinear  optical  phase  conjugation  was  first  ob¬ 
served  in  the  backward  scattered  light  produced  by  stimulated  Brillouin  scatter¬ 
ing  (SBS)  [2.21,  22].  No  pump  beams  are  required  for  phase  conjugation  by 
SBS,  but  the  intensity  of  the  incident  wave  must  exceed  a  threshold  value  (~  1 
MV/cm2)  for  the  process  to  occur.  Nonlinear  photon  echoes  [2.23]  in  an  atomic 
medium  provide  another  technique  for  generating  phase-conjugate  waves;  how¬ 
ever,  atomic  motion  has  a  serious  adverse  effect  on  phase  conjugation.  Three- 
wave  mixing  [2.24]  in  asymmetric  crystals  can  also  yield  wave-front  reversed 
replicas,  but  the  phase  matching  condition  restricts  the  angular  field  of  view  of 
the  input  wave  to  very  small  values.  This  restriction  is  eliminated  in  phase  con¬ 
jugation  by  four-wave  mixing  (FWM)  [2.25-27],  a  process  that  will  be  studied 
in  detail  in  this  chapter.  FWM  uses  the  optically  induced  polarization  of  a 
nonlinear  medium  that  is  cubic  in  the  electric  field  strength,  and  can  occur  in 
all  isotropic  materials  like  gases,  liquids,  and  glasses,  as  well  as  crystals.  In  this 
chapter  we  will  restrict  our  attention  to  nonsaturable,  transparent  Kerr  media. 
We  will  consider  the  frequency  of  all  the  interacting  optical  waves  to  be  equal 
and  far  removed  from  any  material  resonance,  and  will  refer  to  this  as  degener¬ 
ate  four-wave  mixing  (DFWM). 

We  start  by  defining  phase-conjugate  waves  and  give  a  simple  proof  of  how 
aberrations  due  to  propagation  in  a  distorting  medium  can  be  corrected  by 
phase  conjugation.  We  then  describe  the  nonlinear  optical  polarization  induced 
in  the  medium  in  the  presence  of  electromagnetic  waves,  and  discuss  FWM  via 
a  third-order,  nonlinear  susceptibility.  The  discussion  is  then  directed  to  the 
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derivation  of  the  coupled  wave  equations  describing  the  spatial  evolution  of  the 
interacting  waves.  Analytic  solutions  for  the  phase-conjugate  reflectivity  are 
obtained  with  negligible  pump  depletion  for  both  equal  and  different  pump 
intensities.  The  case  of  pump  depletion  is  solved  exactly  for  transmission  and 
reflection  geometries.  We  also  consider  mutually  pumped  phase  conjugation 
(MPPC)  in  Kerr  media  with  electrostrictive  nonlinearity.  We  conclude  the 
chapter  with  a  summary  of  our  results. 


2.2  Phase-Conjugate  Electromagnetic  Waves 

Consider  the  propagation  of  an  electromagnetic  wave  along  the  forward  z- 
direction.  Its  electric  field  can  be  written  as 

E  =  A  cos  (cot  —  kz  —  <f>),  (2.1) 

where  co  is  the  frequency  and  k  is  the  wave  number  of  the  monochromatic 
optical  wave.  The  amplitude  A  and  the  phase  (j>  are  real  functions  of  posi¬ 
tion  (x,y,z).  If  the  amplitude  is  a  slowly  varying  function  of  z  compared  with 
cos  (cot  —  kz  —  <f>),  then  the  propagation  of  the  wave  can  be  easily  understood  in 
terms  of  the  motion  of  the  wave  fronts:  three-dimensional  surfaces  defined  by 
kz  +  <j>(r)  =  constant.  The  phase-conjugate  of  an  electromagnetic  wave  defined 
by  (2.1)  is  given  by 

Ec  =  A  cos(cut  +  kz  +  <j>).  (2.2) 

By  examining  (2.1)  and  (2.2)  we  find  that  these  two  waves  have  exactly  the  same 
wave  fronts  at  any  point  in  space.  However,  the  motion  of  these  two  sets  of  wave 
fronts  is  in  opposite  directions.  In  fact,  if  we  reverse  the  sign  of  the  time  variable 
t  in  E  we  obtain  Ec,  and  thus  the  phase-conjugate  wave  is  often  referred  to  as 
the  time-reversed  wave. 

Using  analytic  representation  we  may  write  (2.1)  as 

E=  A(r)Qii(0,-kz~^r)\  (2.3) 

where  only  the  real  part  of  the  right-hand  side  is  the  physical  value  of  the  electric 
field.  If  we  further  define  a  complex  amplitude 

Ai(r)  =  A(r)e-W'\  (2.4) 

then  the  electric  field  of  (2.3)  may  be  written  in  analytic  representation  as 

E  =  Al(rW{a'~kz).  (2.5) 

Using  these  definitions  the  phase-conjugate  wave  given  by  (2.2)  can  be  written  as 

£c  =  ^(r)ei("'+fc*>.  (2.6) 

We  note  that  the  phase-conjugate  wave  propagates  along  the  backward  z-direc- 
tion  with  a  complex  amplitude  A1{.  To  get  Ec  from  E  we  take  the  complex 
conjugate  of  the  spatial  part  only,  leaving  intact  the  factor  exp(icor). 
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Distorting 

Medium 


z  =  o  z  =  L 
(a)  Incident  Wave 


z  =  o  z  =  L 
(b)  Phase  Conjugate  Wave 
Fig.  2.1.  Propagation  of  phase-conjugate  waves 


To  appreciate  the  property  of  phase  conjugation  we  consider  the  propaga¬ 
tion  of  a  plane  wave  through  a  distorting  medium  such  as  the  atmosphere 
(Fig.  2.1a).  Let  the  distorting  medium  be  confined  to  the  region  between  z  =  0 
and  z  =  L.  As  a  result  of  the  nonuniformity  of  the  refractive  index  n(x,  y,  z)  of 
the  distorting  medium,  the  wave  fronts  of  the  transmitted  wave  are  no  longer 
planar.  Let  the  electric  field  of  the  beam  propagating  in  the  forward  z-direction 
be  represented  by  (2.5).  The  wave  equation  in  MKS  units  is  given  [2.1]  by 


V2E 


d2E  dE  d2PNL 
^°£~dt2  dt2  ’ 


(2.7) 


where  /i0  is  the  permeability  of  the  nonmagnetic  medium  and  a  is  its  conductiv¬ 
ity,  which  is  included  to  account  for  possible  absorption  losses.  e(x,  y,  z)  is  the 
dielectric  constant  of  the  distorting  medium  and  PNL  is  the  nonlinear  polariza¬ 
tion  induced  in  the  medium.  In  deriving  the  wave  equation  (2.7),  we  have  as¬ 
sumed  that  the  fractional  change  of  s  in  one  optical  wavelength  is  small.  If  we 
consider  the  propagation  of  the  wave  given  by  (2.5)  in  a  medium  with  negligible 
absorption  (a  —  0)  and  negligible  nonlinear  effects  (PNL  =  0),  then  its  complex 
amplitude  A  x  must  satisfy  the  following  equation: 
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1  + 


CD  ? 

—  nz(x,y,z)  -  k 
c 


-2  ik  —  A. 
dz  1 


0, 


(2.8) 


where  the  subscript  1  denotes  spatial  derivatives  transverse  to  the  z-axis. 
In  deriving  (2.8)  we  have  used  the  slowly  varying  envelope  approximation 
(\kjdAj/dz\  »  | d2Aj/dz2\),  and  the  relation  e  =  e0n2.  This  is  a  first-order  differen¬ 
tial  equation  in  z;  thus,  if  Al  is  known  at  z  =  z0  for  all  x  and  y,  then  Ax  is 
determined  for  all  points  in  space,  provided  that  n{x,y,  z)  is  known. 

We  now  consider  the  complex  conjugation  of  (2.8)  as  a  purely  mathematical 
operation,  and  assume  that  the  refractive  index  of  the  distorting  medium  is  real 
li.e.,  the  medium  is  lossless): 


V\A\  + 


or 


n2(x,y,z)  -  k 


1  A* 


*  +  2ik  —  A*  =  0. 
dz  1 


(2.9) 


This  is  a  scalar  wave  equation  and  is  obeyed  by  a  wave  propagating  in  the 
backward  z-direction  with  the  form 


i  (<ot+kz) 


provided  that  we  put 


(2.10) 


A2  =  pAf  for  all  x,  y,  z,  (2.11) 

where  p  is  an  arbitrary  constant.  We  have  thus  shown  that  a  wave  A2 
exp[i(cot  +  kz)],  whose  spatial  portion  is  everywhere  the  complex  conjugate  of 
/Cexp[i(cof  -  fcz)],  satisfies  the  same  scalar  wave  equation.  Since  the  scalar 
wave  equation  is  first  order  in  z,  A2(x,y,z)  can  be  determined  for  all  x,  y  and  z 
provided  that  A2  is  known  at  some  z  =  z0  for  all  x  and  y.  In  addition,  if 

A2(x,y,z0)  =  pAf(x,y,z0)  for  all  x,  y,  (2.12) 

then 


A2(*,y,z)  =  pAf(x,y,z)  for  all  x,  y,  z.  (2.13) 

Physically,  this  means  that  if  in  some  plane  z  =  z0  (say  z0>  L  in  Fig.  2.1)  we 
generate  a  wave  Ec  whose  amplitude  is  the  complex  conjugate  of  that  off,  then 
Ec  will  propagate  backward  and  remain  the  phase  conjugate  of  E  everywhere. 
Its  wave  fronts  will  thus  coincide  with  those  of  E  at  all  points  in  space.  Figure 
—  lb  illustrates  that  the  propagation  of  Ec  is  free  of  aberrations  after  traversing 
the  distorting  medium. 


2.3  Phase  Conjugation  by  Externally  Pumped  DFWM 

2.3.1  Coupled  Wave  Equations 

We  consider  transparent  Kerr  media  that  are  nonsaturable  and  nondispersive. 
In  phase  conjugation  by  FWM  [2.25-27]  (Fig.  2.2),  three  optical  waves  are 
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PUMP  PHASE  CONJUGATE 


Fig.  2.2.  The  basic  geometry  of  phase  conjugation  via  FWM 


incident  simultaneously  on  the  nonlinear  medium:  two  intense,  counterprop- 
agating  pump  beams  £t  and  £2;  and  a  third,  weak  probe  beam  £4  at  a  small 
angle  (  =  20)  to  the  pump  beams.  Owing  to  nonlinear  interactions  in  the  medium, 
a  fourth  beam  £3  is  generated  which  is  the  phase  conjugate  of  the  probe  beam. 
The  frequency  of  each  interacting  optical  wave  is  equal  to  co,  and  is  far  removed 
from  any  material  resonance.  If  we  assume  that  all  the  fields  are  linearly  polar¬ 
ized  in  the  same  direction,  then  the  total  electric  field  £  is  given  by 

E(r,  0  =  Z  Ej(r>  0  =  x  Z  AAZ)  exP  [*(«>*  “  ki ' ')]  +  c-c->  (2- ] I4) 

j= i  z  J= i 

where  Ai  is  the  complex  amplitude  of  the  j th  wave  in  the  steady  state,  and  c.c. 
denotes  the  complex  conjugate.  The  z-axis  is  defined  to  be  normal  to  the  surface 
of  the  medium,  with  the  complex  amplitudes  assumed  to  be  functions  only  of  z 
owing  to  energy  transfer  by  beam  coupling.  The  polarization  induced  in  the 
lossless,  isotropic  Kerr  medium  is  in  the  same  direction  as  the  electric  fields,  and 
is  given  by 

P  =  x(1)£  +  x{2)E2  +  x<3)£3  +  "'»  (2.15) 

where  x{n)  is  the  nth-order  electric  susceptibility  of  the  medium. 

In  the  following  discussion  we  will  consider  the  polarization  term  propor¬ 
tional  to  x(3)  only,  since  it  can  generate  the  phase-conjugate  wave  by  DFWM. 
Of  all  the  possible  third-order  terms  that  arise  from  £3  we  are  interested  only  in 
terms  which  satisfy  the  frequency  and  momentum  relation  for  phase  conjuga¬ 
tion  by  FWM:  co1  +  a>2  =  co3  +  co4  =  2 co,  and  +  k2  =  k3  +  k4.  The  quan¬ 
tum-mechanical  picture  of  the  process  involves  annihilation  of  a  photon  from 
each  pump  wave  with  the  simultaneous  creation  of  a  photon  for  both  the  probe 
and  the  phase-conjugate  wave.  The  wave  equation  (2.7)  is  satisfied  separately  by 
each  frequency  component  of  the  total  electric  field.  Using  the  slowly  vary¬ 
ing  envelope  approximation  we  obtain  the  following  coupled  wave  equations 
[2.28-31] 
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c IA1 
dz 

d  A2 
dz 

dA3 

dz 

dA4 


=  -*At  -  ifc|\41(3|/l1|2  +  6\A2\2  +  6\A3\2  +  6|/44|2)  +  6A$A3A41, 

—  %A  2  +  \k[A2(6\A1\2  +  2\A2\2  +  6\A3\2  +  6|  ^l^]2)  +  6/4*zl3/l4], 

(2.16) 


=  <xA3  +  i/c[/l3(6|/l1|2  +  6\A2\2  +  3|/13|2  +  6\A4\2)  +  6AyA2AX\, 
dz  =  —ocA4  -  k[/l4(6|/l1|2  +  6\A2\2  +  6\A3\2  +  3M4|2)  +  6A,A2A*3l 


where 


8  cos  9 

In  arriving  at  (2.16)  and  (2.17)  we  have  assumed  that  the  z-axis  coincides  with 
the  bisectrix  of  the  beams.  The  first  term  on  the  right-hand  side  of  each  equation 
of  (2.16)  represents  linear  absorption  in  the  Kerr  medium.  In  isotropic  media 
the  phase-matched,  third-order  polarization  terms  in  the  above  equations  have 
different  origins:  phase  modulation,  both  self  and  crossed  (denoted  by  PM), 
transmission  gratings  (T),  reflection  gratings  (R),  and  two-photon  excitation 
(2P).  The  relative  magnitude  of  each  process  depends  strongly  on  the  non¬ 
linear  medium,  and  on  the  polarization  and  the  relative  angles  of  the  interacting 
waves.  Writing  the  contribution  due  to  each  process  separately,  we  obtain  the 
following  coupled  wave  equations: 


U/l  .  _ 

=  —  a/^i  ~  2i/cPM/41(|,41|  +  \A2\  4-  |/13|  +  |zl4|2) 

-  2ifcTzl4(/l1/l*  +  A%A3)  -  2i/cRzl3(/l1/l*  +  A$A4) 

2ikrA2(A  1  A* )  ix 2p[_A* (A l A i)  +  2A*(AlA2  -1-  A3A4) 

+  2A$(AlA3)  +  2At(AlA4)l 

d  A2 

~  aA2  +  2iKTPMv42(| >1 1|  +  \A2\2  +  |zl3|2  +  |/44|2) 

+  2iKTA3(A1[A4  +  A2A$)  +  2iKRA4(A*A2  +  A2A%) 

+  2i^zl1(zl?zl2)  +  \K2p\lA*{A,A2  +  A3A4)  +  A*2{A2A2) 

+  2A*(A2A3)  +  2A*(A2A4)2,  (2.18) 

d/I3 

^2  =  a^3  +  2i/cPM/43(|y41|2  +  \A2\2  +  \A3\2  +  |/44|2) 

+  2\KTA2(AyA1  +  A$A3)  +  2iKRAr  {A\A3  +  A2A%) 

+  2iK'RA4{A3AX)  +  iK2P[2/lt(/l1zl3)  +  2A${A2A3)  +  A*3{A3A3) 

+  2At(AlA2  +  A3A4)l 
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-j- -  =  —ocA4  —  2itcPM/l4(|/l1|2  +  \A2\2  +  \A3\2  +  \A4\2) 
d  z 

-  Hk^AMX^a  +  A2A$)  -  2ikrA2(A1A*  +  A*A4) 

-  2\k'rA3(A*A4)  -  iKn&AXiA^)  +  2A*(A2A4) 


+  2A$(AlA2  +  A3A4)  +  A*(A4A4)l 

where 

_  S 0  .  _  IF  ,  _  Q>Xr3)  IF 

KpM  8  cos  0  \j  e  ’  Kj  8  cos  0  \J  e  Kr  8  cos  6  \j  e 

_  c°Xr3)  IF  ,  _  w'/jp  IF 

Kr  8  cos  0  v  e  K2P  8  cos  0  v  £  ' 


(2.18) 


(2.19) 


Note  that  if  all  the  third-order  electric  susceptibilities  are  equal  (x(Pm  =  Zt3>  = 
Xr3>  =  Xr(3>  =  X(2P  =  X<3>)  then  we  recover  the  coupled  wave  equations  of  (2.16). 

The  first  nonlinear  term  on  the  right-hand  side  of  each  equation  in  (2.18) 
represents  self  and  crossed  phase  modulation  terms,  so  called  because  these 
terms  modify  the  phase  of  each  wave  without  contributing  to  the  energy  ex¬ 
change.  The  next  two  mechanisms,  with  coupling  constants  jct  and  kr,  kr,  can 
be  viewed  as  a  holographic  process  as  shown  in  Fig.  2.3.  Transmission  gratings 
are  formed  by  the  interference  of  the  forward  pump  with  the  probe  beam,  and 
the  backward  pump  with  the  phase-conjugate  beam,  both  gratings  having  large 
spacing  and  the  same  wavevector.  Reflection  gratings  are  formed  by  the  interfer¬ 
ence  of  the  forward  pump  with  the  phase-conjugate  beam,  and  the  backward 


GRATINGS 

E1*E4,E2E3* 


GRATING 

Ei*E2 


2k  REFLECTION 
GRATING 

e3*e4 


Fig.  2.3.  Gratings  formed  in  FWM 
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pump  with  the  probe  beam,  both  gratings  having  the  same  wavevector  and 
smaller  spacing  than  the  transmission  gratings.  Two  additional  reflection  grat¬ 
ings  (known  as  the  ‘ 2k ’  gratings  because  of  the  magnitudes  of  their  grating 
waveveciors)  are  formed  by  both  pairs  of  counterpropagating  beams  in  the 
FVVM  geometry:  one  pair  is  formed  by  the  oppositely  travelling  pump  waves, 
and  the  other  is  formed  by  the  probe  and  phase-conjugate  waves.  The  ‘ 2k ' 
gratings  have  a  very  fine  grating  spacing  and  the  response  of  the  medium  at  such 
high  spatial  frequencies  is  usually  small.  The  last  mechanism  has  no  holographic 
analog,  and  corresponds  to  a  nonlinear  index  that  has  no  spatial  modulation 
but  temporal  oscillations  at  frequency  2 co.  It  arises  because  of  self-  or  induced 
two-photon  absorption.  Readout  (or  interaction)  of  the  nonlinear  index  with  a 
third  wave  generates  the  desired  polarization. 

The  formidable  task  of  solving  the  complete  nonlinear  coupled  equations 
(2.16)  in  the  presence  of  linear  absorption  and  pump  depletion  in  the  medium, 
with  the  effects  of  nonlinear  refraction  included,  was  undertaken  numerically  by 
Guha  and  Conner  [2.31].  Their  method  does  not  require  the  determination  of 
the  first  integrals  of  the  coupled  nonlinear  differential  equations,  and  the  effects 
of  linear  absorption  can  be  incorporated  in  the  theory  in  a  straightforward 
manner.  Marburger  and  Lam  [2.28]  also  examined  this  problem  for  collinear 
geometry  of  all  the  interacting  waves,  with  the  pumps  polarized  orthogonally  to 
the  probe  and  phase-conjugate  beams.  This  results  in  the  elimination  of  terms 
with  coupling  constants  jct  and  kr,  as  well  as  two  terms  with  coupling  constant 
k'2P,  on  the  right-hand  side  of  the  coupled  wave  equations  (2.18)  because  the 
scalar  product  of  the  two  fields  is  zero  for  this  choice  of  polarization.  For  negli¬ 
gible  absorption,  the  equations  have  also  been  solved  by  Maruani  [2.29], 
Ujihara  [2.30],  and  Lytel  [2.33].  Kaplan  and  Law  [2.32]  showed  the  existence 
of  multistable  roots  and  isolas  in  the  input-output  response  of  FWM  in  Kerr 
media. 


2.3.2  Analytic  Solutions  with  Negligible  Absorption  and  Undepleted  Pumps 

We  now  limit  our  attention  to  the  case  of  negligible  absorption  in  the  medium 
(a  =  0),  and  equal  third-order  electric  susceptibilities  for  all  the  processes.  We 
also  assume  that  the  pump  fields  are  much  stronger  than  the  probe  and  conju¬ 
gate  fields,  so  that  |i4x|,  \A2\  »  \  A3\,  \A4\.  The  coupled  wave  equations  (2.16)  then 
simplify  to 

dA,  , 

-^=  -3kA1(|/l1|2  +  2|A2|2), 

dA2 

=  3ikA2(2\A  t\2  +  |-42l2), 

dA* 
dz 


6iK['43(|/41|2  +  \A2\2)  +  AlA2A*2, 
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+  \A2\2)  +  AxA2A%\  (2.20) 

The  terms  involving  the  square  of  the  modulus  of  the  complex  amplitudes  on 
the  right-hand  side  modify  the  phases  of  the  waves  without  affecting  the  energy 
exchange  between  the  waves.  In  the  limit  of  strong  pumps  the  pump  equa¬ 
tions  are  not  coupled  to  the  weak  probe  and  its  phase  conjugate,  and  can  be 
solved  separately.  One  can  show  from  the  pump  equations  that  dlT^/dz  = 
d|/l2|2/dz  =  0,  so  that  as  the  pumps  propagate  through  the  medium  their  in¬ 
tensities  remain  constant  at  \A  iol2  and  \A2L\2  respectively,  using  the  notation 
^(O)  =  Al0,  and  A2(L)  =  A2L.  When  the  electric  fields  of  the  optical  waves  are 
described  by  equations  (2.14),  then  the  intensity  of  the  ;'-th  wave  is  equal  to 
e.0nc\Aj\2l2  in  MKS  units  [2.5],  Ignoring  the  numerical  factor,  we  use  \Aj\2 
to  represent  the  intensity  of  the  j-th  wave.  The  pump  phases,  however,  vary 
linearly  with  distance  owing  to  the  change  in  their  propagation  vectors: 

Ay{z)  =  T10exp[-3i/cz(|T10|2  +  2|,42L|2)], 

(2.21) 

A2(z)  =  A2Le\p[3iK{z  -  L){2\Al0\2  +  \A2L |2)]. 

The  product  of  the  two  pump  amplitudes  is  now  given  by 

A^z)A2{z) 

=  ^  10^21, exP[3iKz(|/410|2  -  |T2L|2)]exp[  — 3i/<L(2|/l10|2  +  \A2L\2)']. 

(2.22) 

Note  that  the  z  dependence  of  the  product  vanishes  for  equal  pump  intensities 
(I  A  jo  |2  =  \A2L\2).  For  unequal  pump  intensities  this  product  will  introduce  a 
phase  mismatch  in  the  coupling  terms  AlA2A*  and  AlA2A*  that  occur  in  (2.20) 
for  waves  3  and  4  respectively  [2.26], 


2.3.2. 1  Equal  Pump  Intensities 


We  now  solve  for  the  probe  and  the  phase-conjugate  waves  assuming  equal 
pump  intensities  (M10|2  =  \A2L \2  =  /10).  Defining  new  variables 

B3  =  T3exp[  — 6i/cz(|/l10|2  +  |/t2L|2)], 

(2.23) 

=  /t4exp[6kz(|/l10|2  +  |T2J2)], 
the  equations  for  A3  and  A4  in  (2.20)  become 


M3 

dz 


=  i/cB|, 


di?4 

dz 


=  -i  >cBf, 


(2.24) 


where 
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k-  =  6/C/410/42Lexp[  — 3iK'L(2|/l10|2  +  |42L|2)].  (2.25) 

Boundary  conditions  specify  the  complex  amplitudes  A3{L)  =  A3L  and  /44(0)  = 
q40  of  the  two  waves  at  their  input  planes  z  =  L  and  z  -  0,  respectively.  The 
solutions  of  (2.24)  have  the  form 


p  M  _  cos(|£'iz)  B  ,  •  *  sin[|/c|(z  -  L)] 
3(Z)  cos(|/c|L)  3L  |/c|  cos(|k|L) 


Bt o, 


54(z)  =  -i 


k  sin(|/c|z) 
|k|  cos(|k|L) 


BL  + 


cos[|/c|(z  —  L)] 
cos(|/c|L) 


B. 


’40- 


(2.26) 


For  phase  conjugation  by  DFWM  the  probe  wave  B40  is  nonzero  at  its  input 
plane  z  =  0,  while  there  is  no  conjugate  wave  at  the  z  =  L  plane,  i.e.,  A3L  =  0  = 
B3L.  In  this  case  the  solutions  in  (2.26)  simplify  to 


sin[|>c|(L-  z)] 
cos(|jc|L) 


Bto, 


S4(z) 


cos[|£-|(L-z)]g 

cos(|/c|L) 


(2.27) 


Note  that  wave  B3,  generated  by  nonlinear  interactions  in  the  medium,  is 
linearly  proportional  to  the  complex  conjugate  of  the  input  probe  wave  B*0. 
Consider  an  input  probe  beam  that  is  not  a  plane  wave  but  possesses  a  complex 
wave  front;  it  can  be  considered  as  a  continuous  superposition  of  plane  waves 
all  having  the  same  frequency  but  different  wave  vectors.  Owing  to  the  linearity 
of  (2.27),  each  of  the  plane  wave  components  of  the  complex  input  wave  front 
is  phase-conjugated  independently,  resulting  in  a  phase-conjugate  wave  front 
travelling  in  the  backward  direction  [2.5,  26].  Hence  faithful  distortion  correc¬ 
tion  is  possible  for  an  input  probe  beam  with  large  angular  field  of  view  if  we  use 
a  phase  conjugator  that  relies  on  FWM  in  a  Kerr  medium.  For  small  coupling 
strength  {\ic\L  «  1),  the  phase-conjugate  wave  at  the  output  plane  simplifies  to 
-43(0)  %  —6iKLAloA2LA*oexp\_  —  2iKL(2\Alo\2  +  \A2L\2)'],  which  shows  how 
the  nonlinear  mixing  of  the  three  input  waves  in  the  medium  led  to  the  genera¬ 
tion  of  the  fourth  wave.  The  phase-conjugate  reflectivity  R  and  the  coherent 
transmissivity  T  are  given  by  [2.1,  5,  26] 


„  M  3(°)l2 

M40l2 
t_\abl)\2 
\a40  I2 


|b3(0)|2 
UW2  ■ 
\BBD\2 

\b,0\2 


tan2(|/c|L), 

1 

cos2(|ic|L)’ 


(2.28) 


The  expression  for  T  shows  that  the  transmitted  probe  beam  is  amplified  for  all 
values  of  pump  intensities  and  coupling  strength,  while  R  >  1  onlyif|/c|L>  n/4. 
When  \  k\L  =  n/2  then  R  =  T  =  oo,  which  corresponds  to  oscillation  -  nonzero 
outputs  of  the  phase-conjugate  wave  and  the  transmitted  probe  wave  exist  at 
the  output  planes  even  in  the  absence  of  an  input  probe  wave  zl40.  This  depen- 
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NORMALIZED  COUPLING  STRENGTH  |k|L 

Fig.  2.4.  Phase-conjugate  reflectivity  R  and  transmissivity  T  as  a  function  of  normalized  coupling 
strength  |\jL. 


dence  of  R  and  T  on  \ic\L  is  shown  in  Fig.  2.4.  In  an  experiment  on  phase 
conjugation  by  FWM  in  a  x(3)  nonlinear  medium  with  interaction  length  L, 
such  a  behavior  would  be  obtained  by  measuring  R  and  T  as  a  function  of  the 
equal  pump  intensity  /10.  For  small  nonlinearity  (i.e.  |k|L«  1),  R  ~  |ic|2L2  = 
36k2L2/20  and  T  «  1  +  \ic\2L2  =  I  +  3 6k2L2I210.  Hence  a  plot  of  phase-conju¬ 
gate  reflectivity  as  a  function  of  pump  intensity  on  a  log-log  scale  would  have 
a  slope  of  2  and  an  intercept  equal  to  2  log(6/<L).  An  experimental  measure¬ 
ment  of  the  phase-conjugate  reflectivity  R  for  known  values  of  the  equal  pump 
intensity  /10  and  the  interaction  length  L  would  yield  the  value  of  the  coupling 
constant  k,  which  in  turn  is  a  known  function  [see  (2.17)]  of  the  nonlinear 
susceptibility  x(3)  of  the  FWM  medium. 

Figure  2.5  is  a  plot  of  the  beam  intensities  /3(z)  and  /4(z),  normalized  to  the 
input  probe  intensity  /40,  as  a  function  of  distance  inside  the  interaction  region 
for \i<\L  =  n/3.  The  conjugate  wave  is  zero  at  the  input  plane  z  =  L,  and  exceeds 
the  probe  input  intensity  at  the  output  plane  z  =  0,  leading  to  a  reflectivity 
greater  than  unity  (R  =  3).  At  the  same  time  the  transmitted  probe  beam  inten¬ 
sity  exceeds  its  input  intensity,  leading  to  a  coherent  transmission  amplification 
(T  =4).  The  amplifications  of  the  probe  and  phase-conjugate  beams  are  at  the 
expense  of  the  pump  beams. 

The  case  of  nearly  degenerate  FWM  (NDFWM),  in  which  the  weak  inci¬ 
dent  probe  field  is  at  a  slightly  different  frequency  to  the  frequency  of  the  two 
pump  beams,  has  been  considered  by  Pepper  and  Abrams  [2.34,  26].  For  weak 
nonlinear  coupling  their  results  show  that  the  phase-conjugate  reflectivity  has  a 
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NORMALIZED  DISTANCE  z  /  L 

Fig.  2.5.  Normalized  beam  intensities  inside  the  interaction  region  for  \k\L  =  nj 3 


sinc2(x)  dependence  on  the  phase  mismatch,  which  is  a  typical  result  for  coher¬ 
ent,  phase-mismatched,  nondepleted  interactions  [2.5]. 


2.3.2.2  Unequal  Pump  Intensities 


The  coupled  wave  equations  (2.24)  for  the  probe  and  the  phase-conjugate  waves 
will  now  have  a  phase  mismatch  factor  due  to  the  unequal  intensities  of  the 
pump  beams: 


~  =  i/ce'^B*, 
dz 


d  B4 
dz 


—  ijce'^'B*, 


(2.29) 


where  the  phase  mismatch  wavevector  /?  is  given  from  (2.22)  by  the  following 
quantity: 

P  =  M\Al0\2  -\A2L\2)l  (2.30) 

Note  that  /?  =  0  for  equal  pump  intensities.  Assuming  the  boundary  conditions 
appropriate  for  phase  conjugation  by  DFWM  ( B3L  =  0,  B40  ^  0),  the  solutions 
of  (2.29)  have  the  following  form: 

B.M--  Mfifee-n.  sinW--z>]  , 

ycos(yL)  —  i^sin(yL) 
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B4(z)  =  £40e^2 


(3 

y  cos  [y(L  -  z)]  +  i-sin[y(L  -  z)] 
B 

ycos(yL)  +  i-sin(yL) 


(2.31) 


where  y  =  ^/| /c|2  +  /?2/4.  If  =  0,  then  y  =  |k|  and  the  solutions  given  by  (2.31) 
reduce  to  those  in  (2.27)  for  equal  pump  intensity.  The  phase-conjugate  reflec¬ 
tivity  R  is  given  by  the  simple  expression  [2.31] 


R  = 


sin2(yL) 


sin2(yL) 


(2.32) 


which  reduces  to  (2.28)  when  [3  =  0.  For  R  to  be  equal  to  infinity  (self-oscilla¬ 
tion),  y2/|/c|2  =  sin2(yL)  which  can  only  be  true  if  {3  =  0,  i.e.,  only  for  equal  pump 
intensities  [2.28,  31].  If  the  coupling  constant  is  much  smaller  than  the  phase 
mismatch  (i.e.,  \k\  « (3/2),  then  the  phase-conjugate  reflectivity  has  a  sinc2(x) 
dependence  [2.26,  34]:  R  x  [|K|Lsin(0L/2)/(£L/2)]2. 

Figure  2.6  is  a  plot  of  the  phase-conjugate  reflectivity  R  as  a  function  of  the 
normalized  phase  mismatch  /?L  introduced  by  the  unequal  pump  intensities,  for 
several  values  of  the  coupling  strength  \k\L.  For  \k\L  =  nil,  self-oscillation  oc¬ 
curs  only  at  (3  =  0;  finite  values  of  imbalance  in  the  pump  intensities  serve  to 
reduce  the  phase-conjugate  reflectivity  from  this  optimum  value.  The  qualitative 
features  of  this  figure  are  similar  to  the  frequency  bandpass-filter  characteristics 
of  the  phase-conjugator  when  the  probe  frequency  has  a  small  detuning  from 
the  pump  frequency  (NDFWM)  [2.26,  34]. 


NORMALIZED  PUMP  IMBALANCE  PL 

Fig.  2.6.  Phase-conjugate  reflectivity  R  as  a  function  of  normalized  pump  imbalance  /?L  for  various 
coupling  strengths. 
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2.3.3  Analytic  Solutions  with  Pump  Depletion  but  Negligible  Absorption 
2.33 A  The  Transmission  Grating 


We  now  consider  a  lossless  (a  =  0)  Kerr  medium  with  negligible  response  at 
high  spatial  frequencies,  so  that  kr  =  k'r  =  0.  In  addition,  we  neglect  two- 
photon  interactions  in  the  medium,  so  that  the  coupled  wave  equations  (2.18) 
simplify  to 

d  /4 

— -  =  2i/CpM/4 ! 70  -  HkjAMi**  +  A*As\ 

dz 


A  A 

-r -  =  2iKmA2I0  +  2\ktA3{A1A*  +  A2At\ 

d  z 
a  a 

—  =  2ifCpMA3/o  +  2\ktA2(A1A*  +  A*A2), 
dz 


(2.33) 


dA4 


dz 


=  -2\kpmA4I0  -  2iKTAl(AfA4  +  A2A*), 


where  I0  is  proportional  to  the  total  intensity  of  all  the  beams: 


/ o=Z  I  A/.  (2.34) 

j= i 

The  differential  equations  (2.33)  verify  that  d/0/dz  =  0,  so  that  the  total  intensity 
is  conserved.  As  in  the  previous  section  the  phase  modulation  terms  in  each 
equation  modify  the  phase  of  the  wave,  and  can  be  eliminated  by  introducing 
new  variables  Bj  in  the  following  manner: 

B,  =  A .  e2iKWo:,  B2  =  A2e~2iKpM'o:, 

(2.35) 

Bz  =  A2e-z,K™,oZ,  B4  =  A4e2lKm'°:. 


In  terms  of  the  new  variables  the  coupled  wave  equations  are: 


dBl 

dz 


2i/cTB4(B1  B*  +  B*B2), 


^7  =  2i ktB3(B*B4  +  B2B$), 
i^  =  2iKTB2(BlBt  +  B*2B3), 


(2.36) 


~=~2  iktB1(B*B4  +  B2B*3). 

These  equations  are  similar  in  form  to  the  coupled  wave  equations  for  DFWM 
in  photorefractive  media  in  transmission  geometry  [2.35,  36],  but  with  a  cou- 
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pling  constant  kt  that  is  purely  real.  As  the  beams  propagate  in  the  nonlinear 
medium  we  find  that  the  following  quantities  are  conserved: 

Ii(z)  +  I4(z)^dlt 


I2{z)  +  I3{z)  =  d2,  (2.37) 

Bl(z)B2(z)  +  B3{z)B4{z)  =  c , 


where  the  constants  and  d2  are  determined  by  the  boundary  values  of  the 
problem,  and  the  integration  constant  c  is  yet  to  be  determined.  The  first  two 
equations  show  that  the  total  forward  and  backward  energies  are  conserved 
separately  at  each  plane,  while  the  last  conservation  law  is  an  outcome  of  the 
reciprocity  theorem  [2.37],  With  the  help  of  these  conservation  laws,  the  four 
coupled  equations  of  (2.36)  can  be  decoupled  to  yield  the  following  two  first- 
order  differential  equations  for  the  new  variables  Bl2  =  By/B £  and  B34  = 
B3/Bf. 


—  2i/cT[c  +  (d1  —  d2)Bi2  “ 

d  B 

j,  =  2iKT[c  +  (d2  —  dl)B34  —  c*£f4]. 


(2.38) 


These  equations  can  be  readily  solved  to  yield  the  following  solutions: 


B  (z)  =  QBlli0)  +  ’  tan(KTQz)l4Bl2(0)  -  2c] 
12  Q-itan(^Tez)[d  +  2c*B12(0)] 

B  (z)  =  2^34(0)  +  i  tan(/cTQz)  [zf534(0)  +  2c] 
34  Q  —  i  tan(KTQz)  [d  —  2c*534(0)] 

where  the  parameters  A  and  Q  are  defined  as 
A  =d2-dl, 

Q  =  (A2  +  4\c\2)1'2. 


(2.39) 


(2.40) 


For  phase  conjugation  by  DFWM  the  boundary  conditions  are  the  nonzero 
values  of  the  two  pump  fields  at  z  =  0  and  z  =  L(Al0  and  A2L,  respectively)  and 
the  probe  field  at  z  =  0  (/140),  while  the  phase-conjugate  field  is  zero  and  its 
input  plane  z  =  L  (A3L  =  0).  With  jB34(L)  =  0,  we  obtain  the  phase-conjugate 
reflection  coefficient  p(  =  B34( 0)): 


2ic  tan(KTQL) 

Q  +  iA  tan  {ktQL)' 


(2.41) 


If  we  assume  that  the  probe  and  phase-conjugate  waves  are  much  weaker  than 
the  two  strong,  equal  intensity  pump  waves  (/3,  /4  «  Ix  =  /2),  then  zf  =  0  and 
Q 2  =  4|c|2.  From  the  third  conservation  law  in  (2.37)  we  find  that  c  =  BXB2,  so 
that  Q  =  2|5162|.  Hence  the  expression  for  p  simplifies  to: 
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p  a  - i  tan<2,cTl  B1  Bjl O-  (2-42) 

1^1  #2l 

which  is  identical  to  the  complex  phase-conjugate  reflection  coefficient  that  can 
be  obtained  from  (2.27),  for  negligible  pump  depletion,  but  with  2 kt  in  the  above 
equation  replaced  by  6k  for  equal  contributions  from  all  the  different  mecha¬ 
nisms  besides  the  transmission  grating.  The  phase-conjugate  intensity  reflec¬ 
tivity  is  defined  as  R  =  \p\2,  and  an  examination  of  (2.41)  reveals  that  R  is  a 
function  of  |c|2  only,  and  not  of  any  linear  combinations  of  c  and  c*.  We  may 
solve  for  |c|2  by  recalling  that  Bl2{L)  =  c/I2L •  Also,  B12(0)  =  /i0/(c*  -  p*U 0), 
where  we  have  used  the  third  conservation  law  of  (2.37).  Substituting  these 
values  of  Bl2  and  p  into  the  first  equation  of  (2.39)  evaluated  at  z  =  L,  we  obtain 
the  following  transcendental  equation  for  |c|2: 

(kl2  -/io^2l)[Q2  +  ^ 2  tan2  (ktQL)]  +  4\c\2 1 2LI40 tan2 {ktQL)  =  0. 

(2.43) 

To  study  the  phase-conjugate  reflectivity  JR  numerically,  we  define  beam  ratios 
at  input  as  follows:  probe  ratio,  q  =  /40/(/ 10  +  Iil)\  pump  ratio,  r  =  I2L/I i0; 
and  I0  =  1.  In  terms  of  these  parameters  the  input  beam  intensities,  normalized 
to  the  total  input  intensity,  are  given  by:  Il0/1 0  =  1/[(1  +  q)(l  +  »*)],  = 

r/[(l  +  q){  1  +  r)],  and  /40//0  =  q/{  1  +  q).  The  argument  of  the  tangent  function 
(  =  KjQL)  in  (2.43)  can  be  expressed  as  GQn,  where  G  =  ktI0L  is  the  net  cou¬ 
pling  strength  in  the  Kerr  medium  and  Q„  =  Q/I0. 

Figure  2.7  is  a  plot  of  the  phase-conjugate  reflectivity  R  in  transmission 
geometry  versus  the  net  coupling  strength  G  for  equal  intensity  pump  beams 
(r  =  1),  and  with  a  probe  as  strong  as  the  sum  of  the  two  pumps  {q  =  1).  The 
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Fig.  2.7.  Phase-conjugate  reflectivity  R  as  a  function  of  net  coupling  strength  G  for  q  =  r  =  1 
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PUMP  RATIO  r 

Fig.  2.8.  Phase-conjugate  reflectivity  R  as  a  function  of  pump  ratio  r  for  G  =  2  and  q  =  1 


reflectivity  R  exhibits  oscillations  as  a  function  of  G,  and  can  even  be  multi¬ 
valued  for  a  suitable  choice  of  the  parameters  q  and  r,  leading  to  bistability 
[2.28,  32].  Smaller  probe  intensity  compared  to  equal  intensity  of  the  strong 
pumps  tends  to  lead  to  instabilities  in  the  behavior  of  R  versus  G.  Figure  2.8  is 
a  plot  of  R  versus  the  pump  ratio  r  for  G  =  2  and  q  —  1.  Note  that  the  reflec¬ 
tivity  is  maximized  when  the  pump  beams  are  of  equal  intensity 


2.3.3. 2.  The  Reflection  Grating 


We  now  assume  that  the  response  of  the  medium  to  two-photon  interactions, 
transmission  gratings,  and  ‘2 k'  gratings  is  negligible,  so  that  k2P  =  jct  =  k'r  —  0. 
The  coupled  equations  (2.18)  then  simplify  to 

&A, 

— — 21/Cpm/Ij/o  —  2iK:R^43(y41  >1*  +  A^Af), 

dA2  „ 

~  ^1KPM-^2^0  2lKR/l4(/4fi43  A2A*), 

<m,  .  .  (2'44) 

2i/cPM/4370  4-  2\KRAfA*A2  +  A2A*), 

d/t4 

=  —  2ikPM/44/0  —  2iKRA2{AiA%  +  A%A4). 


Once  again,  these  differential  equations  verify  that  the  total  intensity  /0  is  a 
conserved  quantity.  Defining  variables  Bj  in  the  same  manner  as  before,  we 
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have: 


-1  =  -HkMB.B*,  +  B*2Ba), 
dz 

^  =  2i  KRB4(B*B3  +  B2Bt), 
dz 

<^±  =  2iKRBl(B*B3  +  B2B*), 

dz 


(2.45) 


which  are  similar  in  form  to  the  coupled  wave  equations  for  DFWM  in  photo- 
refractive  media  for  reflection  geometry  [2.35,  38,  39]  but  with  a  coupling  con¬ 
stant  k*r  that  is  purely  real.  As  the  beams  propagate  in  the  nonlinear  medium  we 
find  that  the  following  quantities  are  conserved: 


Mz)  +  U{z)  =  d„  /2(z)  +  /3(z)  =  d2,  Mz)  -  /3(z)  =  </3, 

I2(z)  -  /4(z)  =  <f4,  B^B^z)  -  B3(z)B4(z)  =  c, 


(2-46) 


where  the  constants  d^d2,  d3  and  d4  are  determined  by  the  boundary  values  of 
the  problem,  and  the  integration  constant  c  is  yet  to  be  determined.  With  the 
help  of  these  conservation  laws  the  four  coupled  equations  (2.45)  can  be  de¬ 
coupled  to  yield  the  following  two  first-order  differential  equations  for  the  new 
variables  Bl2  =  B{/B*  and  U34  =  B3/B *: 


-r ~  =  2ixR[c  +  (d,  -  d2)Bl2  -  c*B2l2l 
dz 

dR 

=  2iKR[c  +  (d,  -  d2)B^  - 


(2.47) 


Comparing  the  equation  for  U34  in  reflection  geometry  with  the  corresponding 
equation  in  transmission  geometry  (see  (2.38))  reveals  that  besides  kr  replacing 
kt,  the  sign  of  A  is  reversed  in  reflection  geometry.  Now  the  phase-conjugate 
reflectivity  R  (  =  |p|2)  is  an  even  function  of  A  (see  (2.41)).  One  can  also  show  that 
the  transcendental  equation  in  reflection  geometry  for  determining  |c|2  is  identi¬ 
cal  to  (2.43).  Hence  the  reflection  and  transmission  geometries  yield  the  same 
phase-conjugate  reflectivities  for  identical  coupling  strengths. 


2.4  Mutually  Pumped  Phase  Conjugation  (MPPC)  in 
Electrostrictive  Kerr  Media 

Recently,  several  papers  [2.40-44]  have  reported  a  new  type  of  phase  conjuga¬ 
tion  process  in  photorefractive  crystals,  in  which  two  mutually  incoherent  laser 
beams  pump  each  other  to  generate  the  respective  phase-conjugate  reflections. 
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The  nonlinear  medium  acts  like  a  phase-conjugate  mirror  to  both  of  the  inci¬ 
dent  beams  and  generates  their  phase-conjugate  outputs  via  NDFWM.  Several 
configurations,  differing  in  the  number  of  internal  reflections  from  the  sides 
of  the  photorefractive  crystals,  were  used  to  demonstrate  this  process.  As  in 
[2.45]  we  refer  collectively  to  these  phase  conjugators  as  mutually  pumped 
phase  conjugators  (MPPCs)  because  of  their  common  property  of  mutual 
pumping.  Although  photorefractive  crystals  are  very  efficient  at  low  operating 
powers,  the  process  is  inherently  slow  due  to  the  need  to  generate  a  large  num¬ 
ber  of  photoinduced  charge  carriers  [2.46].  In  [2.47]  the  authors  have  examined 
the  possibility  of  the  same  process  occurring  in  a  transparent  Kerr  medium  with 
an  electrostrictive  nonlinearity,  since  it  offers  the  advantage  of  faster  response 
times. 


2.4.1  Theoretical  Formulation 

The  geometry  for  the  nonlinear  interaction  responsible  for  MPPC  is  shown  in 
Fig.  2.9.  Input  beams  2  and  4  are  a  pair  of  mutually  incoherent  beams  that  enter 
the  medium  from  opposite  sides.  Both  beams  may  be  obtained  from  a  laser  with 
a  short  coherence  length,  such  that  the  difference  between  the  propagation  dis¬ 
tances  from  the  laser  to  the  Kerr  medium  for  the  two  incident  beams  exceeds  the 
coherence  length  of  the  laser.  In  a  transparent  Kerr  medium  like  CS2,  acoustic 
phonons  that  are  initially  present  in  the  medium  owing  to  thermal  and  quantum 
noise  will  help  in  Stokes  scattering  the  incident  photons  via  acousto-optic  Bragg 
scattering  [2.22],  Also,  a  photon  from  each  input  beam  can  be  spontaneously 
converted  into  a  frequency-downshifted  Stokes  photon  and  an  acoustic  phonon. 
The  interference  between  each  incident  beam  and  its  coherently  generated  scat¬ 
tered  beam  will  produce  a  set  of  fringes  which  travel  at  the  acoustic  velocity,  and 
which  drive  an  acoustic  wave  by  electrostriction.  This  additional  sound  wave 
reinforces  the  initial  acoustic  phonons,  thereby  increasing  the  coupling  between 
each  input  beam  and  its  scattered  beam.  The  same  acoustic  wave  may  also 


Fig.  2.9.  Schematic  diagram  showing  the  geometry  of  NDFWM  in  electrostrictive  Kerr  media 
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diffract  the  other  incoherent  input  beam  to  generate  the  cross-readout  phase- 
conjugated  beam  at  each  port,  labelled  in  Fig.  2.9  as  beams  1  and  3.  Thus,  when 
we  inject  two  input  beams,  we  will  observe  the  buildup  of  a  pair  of  Stokes 
downshifted  phase-conjugate  beams.  Interestingly,  one  beam  by  itself  will  not 
generate  a  phase-conjugate  reflection,  and  both  the  incident  beams  are  required 
for  the  nonlinear  medium  to  phase  conjugate  them  in  MPPC.  Since  the  read- 
write  processes  occur  simultaneously,  the  dynamics  of  the  beam  coupling  are 
described  by  NDFWM  equations  in  the  nonlinear  medium  with  appropriate 
boundary  conditions. 

In  Fig.  2.9  beams  2  and  4  are  the  two  input  beams  of  frequency  co,  so  that 
aj2  z=  a)4  =  co.  Beams  1  and  3  are  the  Stokes-scattered  fields  from  the  input  fields 
4  and  2  respectively,  hence  their  frequency  is  downshifted  and  given  by  col  = 
o3  =  co  —  Q,  where  Q  is  the  acoustic  frequency.  This  acoustic  frequency  is  equal 
to  {2n0cov  sin  0)/c,  where  v  is  the  speed  of  sound  in  the  medium,  and  26  is  the 
angle  between  the  beams  1  and  4,  and  between  beams  2  and  3,  respectively,  in 
the  medium.  The  wave  equation  (2.7)  is  satisfied  separately  by  each  frequency 
component  of  the  total  electric  field.  Using  the  slowly  varying  envelope  approxi¬ 
mation,  and  neglecting  linear  absorption  in  the  medium,  we  obtain  the  follow¬ 
ing  coupled  wave  equations: 

d/4 

— -  =  —  2i/cPM/41/0  —  2iKjA4_(A1A*  +  A*A3e  lAk~), 
dz 


-j  —  2ix-pM/l2/o  +  2iK*A3(A*A^e  lAk~  +  A2A*), 
az 

—  =  2ifcPM/43/0  +  2\KTA2(AlA%e'Jk-  +  A2A3), 

az 

=  —  2iK'PM/44/0  —  T\k*A  ^{A*A^  +  A2A*  e>Ak~), 


(2.48) 


where  we  have  assumed  that  the  angle  between  the  coherent  beam  pairs  is  small 
(cos  9  %  1)  in  calculating  the  phase  mismatch.  Since  Q  «  co  we  have  ignored  the 
small  frequency  differences  between  the  waves  and  taken  coj  ^  co  for  all  the 
waves  in  the  definition  of  the  coupling  constants  (see  (2.19)).  The  wavevector 
mismatch,  Ak  =  (2 nQ/c),  is  the  phase  mismatch  associated  with  the  geometry 
shown  in  Fig.  2.9.  For  media  whose  Kerr  effect  is  dominated  by  electrostriction 
the  complex  third-order  susceptibility  can  be  expressed  as  [2.22,  47,  48] 


xV2,= 


~KYe 


4p0(Q2 


v2k2  +  i  nry 


(2.49) 


where  K  is  the  acoustic  wave  number,  ye  is  the  electrostrictive  coefficient,  p0  is 
the  mass  density,  and  r  is  the  inverse  of  the  phonon  lifetime.  The  coupling 
constant  kt  is  now  complex,  and  is  a  function  of  the  frequency  difference  be¬ 
tween  each  incident  wave  and  its  Stokes  scattered  wave.  At  resonance,  Q  = 
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±vK,  and  jcx  is  purely  imaginary  (  =  ig),  indicating  a  ±90°  phase  shift  between 
the  index  grating  and  the  intensity  pattern  [2.48].  If  the  Bragg  condition  is 
satisfied  for  the  incident  wave,  the  sound  wave,  and  the  scattered  wave,  then  a 
phase  shift  of  90°  is  automatically  realized  for  a  wave  scattered  with  a  frequency 
downshift  (Stokes  scattering). 

The  coupled  wave  equations  in  (2.48)  describe  MPPC  in  a  transparent  Kerr 
medium  where  depletion  of  the  incident  beams  has  been  taken  into  account.  The 
equation  for  each  beam  has  contributions  from  two  terms:  one  due  to  diffraction 
from  the  grating  written  by  the  beam  itself  (the  usual  two-beam  coupling  term 
that  arises  due  to  acousto-optic  Bragg  scattering  for  MPPC  in  Kerr  media),  and 
the  second  due  to  cross-readout  of  the  grating  written  by  the  other  pair  of 
coherent  beams  (the  FWM  term).  Both  of  the  self-generated  beams  1  and  3  are 
amplified  as  they  propagate  through  the  interaction  region,  since  both  are  at  the 
Stokes  frequency  of  their  respective  forward-going  pumps  (recall  that  beam  3  is 
propagating  in  the  negative  z-direction). 

We  assume  that  the  phase  mismatch  is  small  for  the  geometry  and  the  Kerr 
medium  under  consideration,  so  that  AkL«  1.  The  coupled  wave  equations  in 
(2.48)  then  reduce  to  equations  (2.33)  for  transmission  geometry,  with  the  differ¬ 
ence  that  jct  =  (i Ig)  now,  owing  to  electrostrictive  resonance.  Hence  we  may  use 
the  solutions  obtained  in  Sect.  2.3. 3.1  with  this  replacement  for  jct.  We  now 
define  the  phase-conjugate  reflection  coefficients  in  the  two  arms  as  the  follow¬ 
ing  amplitude  ratios:  p  =  B3(0)/B^(0)  =  £34( 0)  and  o  =  Bj(L)/Bf(L)  =  Bl2{L). 
For  MPPC  the  boundary  conditions  are  the  nonzero  intensities  of  the  two  input 
beams  at  the  opposite  faces  of  the  nonlinear  medium  (i.e.,  I2L,  /40  ^  0),  while 
the  two  self-generated  phase-conjugate  beams  are  zero  at  their  inputs  (/10  = 

1 31,  =  0).  It  follows  that  B12(0)  =  jB34(L)  =  0,  so  that  the  solutions  given  in  (2.39) 
simplify  to 

_  2c  tanh  pL  2c  tanh  pL 

P  Q  —  A  tanh  pL’  a  Q  +  A  tanh  pL’  (2.50) 

where  p  =  gQ,  A  =  I2i  —  /40,  and  I0  =  I2L  +  /40  for  MPPC.  The  phase-conju¬ 
gate  intensity  reflectivities  in  the  two  arms  are  defined  as  R  =  \p\2  and  S  =  |oj2. 
We  still  have  to  determine  the  integration  constant  c.  An  examination  of  (2.50) 
reveals  that  R  and  S  are  functions  of  |c|2  only,  and  not  of  any  linear  combina¬ 
tions  of  c  and  c*.  We  may  solve  for  |c|2  by  recalling  that  /43(0)A4(0)  =  c  from 
the  boundary  conditions  and  from  the  third  conservation  law  in  (2.37).  This 
result  can  also  be  written  as  p  =  c//40.  Equating  this  expression  for  p  with  that 
given  in  (2.50)  we  finally  arrive  at  the  following  transcendental  equation  for  jc|2: 

tanh(g()L)  =  (2.51) 

Once  the  root  of  the  nonlinear  transcendental  equation  is  determined,  the  inten¬ 
sity  reflectivities  can  be  easily  calculated  from  the  relationships  R  =  |c|2//|0  and 
S  =  |c|2//2l.  It  is  interesting  to  note  that  R/S  =  /2L//40,  so  that  the  ratio  of  the 
two  phase-conjugate  reflection  coefficients  is  equal  to  the  ratio  of  the  two  input 
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beam  intensities.  For  equal  intensity  input  beams  (I2L  =  /40),  A  =  0,  and  Q2  = 
4jr|2  from  (2.40).  Since  the  phase-conjugate  reflectivity  of  either  beam  is  propor¬ 
tional  to  |c|2,  and  hence  also  to  Q ,  the  value  of  Q  must  be  small  for  weak 
phase-conjugate  signals.  The  threshold  for  the  coupling  strength  is  then  ob¬ 
tained  by  using  tanh(x)  %  x  in  (2.51): 


(#f  0^)ih  =  ^th  —  1, 


(2.52) 


where  2 g  is  the  Kerr  amplitude  coupling  coefficient  in  equations  (2.48).  This  is 
lower  than  the  typical  threshold  values  of  25  quoted  for  phase  conjugation  by 
SBS  [2.22]  and  of  10  quoted  for  Brillouin-enhanced  FWM  (BEFWM)  [2.49-52]. 
Equation  (2.52)  is  valid  for  equal  intensity  input  beams.  For  unequal  beam 
intensities  the  threshold  coupling  strength  may  be  obtained  in  the  following 
manner.  Since  4|c|2  =  Q2  -  A1  from  (2.40),  Q2  >  A2  for  real  values  of  |c|2  and 
phase-conjugate  reflectivities.  The  threshold  condition  corresponds  to  Q  =  A,  in 
which  case  the  phase-conjugate  reflectivity  is  zero.  Substituting  this  in  (2.51)  and 
simplifying,  we  obtain 


1 40 


(2.53) 


Figure  2.10  shows  the  variation  of  the  threshold  coupling  strength  with  input 
beam  ratio  The  threshold  value  of  the  coupling  strength  is  a  mini¬ 

mum  for  equal  intensity  input  beams.  Unequal  intensities  of  the  two  input 
beams  require  higher  values  of  the  coupling  strength  for  threshold  operation. 
For  larger  values  of  the  coupling  strength  the  product  gI0L  >  1,  and  R  and  S 
reach  their  limiting  values  of  /2z.//40  and  /40//2L  respectively  in  the  absence  of 
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Fig.  2.10.  Threshold  coupling  strength  Gth  as  a  function  of  input  beam  intensity  ratio  U0/I2l 
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absorption  losses.  Knowing  Bl2(z)  and  B34(z),  the  intensities  of  the  various 
beams  as  a  function  of  the  distance  within  the  interaction  volume  may  be  calcu¬ 
lated  from  the  conservation  laws: 


h(z)  =  Ini*) 


dz  ~  d i  ^34(2) 
l-/12(z)/34(z)’ 


^2  djJjA?) 

1  —  ^  12(Z)^  34-(z) 


13(2)  =  /34W 


dj  ~  ^2^12(Z) 

1  ~  ^  12(Z)^34(Z)’ 


d  1  —  d21 12  (z) 
l-/12(z)/34(z)’ 


(2.54) 


where  /12(z)  =  \Bl2(z)\2  and  /34(z)  =  |£34(z)|2.  The  known  parameters  for  this 
problem  are  the  net  power  flux  A;  the  Kerr  field  coupling  coefficient  2 g;  and  the 
interaction  length  L.  We  now  discuss  some  numerical  results  obtained  in  our 
study  of  this  problem. 


2.4.2  Numerical  Results 

We  consider  CS2  as  a  typical  transparent  Kerr  medium  with  electrostrictive 
nonlinearity.  The  field  Brillouin  gain  coefficient  for  backward  SBS  in  CS2  is 
about  0.1  cm/MW  [2.22,  52],  so  that  for  a  half-angle  of  6  %  10°  in  the  forward 
direction,  the  amplitude  Kerr  coupling  coefficient  2 g  is  2  cm/MW  [2.48].  The 
Brillouin  frequency  of  CS2  near  room  temperature  is  Q  =  7.6  GHz  [2.52].  We 
choose  an  interaction  length  L  =  10  cm.  From  (2.52)  it  follows  that  the  thresh¬ 
old  value  of  the  total  intensity  /0  is  0.1  MW/cm2.  Figure  2.11  is  a  plot  of  the 
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Fig.  2.11.  Phase-conjugate  reflectivities  R  ( - )  and  S  ( - )  in  both  the  arms  versus  the  input 

beam  ratio 
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phase-conjugate  intensity  reflectivities  R  (solid  curves)  and  S  (dashed  curves)  in 
both  the  arms  versus  the  ratio  of  the  input  beam  intensities  q  =  /40//2 l-  F°r  a 
given  beam  ratio,  R  and  S  are  higher  for  a  larger  net  Brillouin  gain  G  since 
this  leads  to  more  coupling  of  beams.  If  the  input  intensity  in  one  arm  is 
much  smaller  than  in  the  other,  which  in  our  study  corresponds  to  /40  «  I2L  or 
/,L  «  /40,  then  the  device  can  phase  conjugate  the  weak  beam  with  a  very  high 
reflection  coefficient,  much  larger  than  unity  (Fig.  2.11).  The  peak  value  of  R 
corresponds  to  all  of  the  power  of  the  strong  input  beam  I2L  being  transferred 
to  the  diffracted  beam  3,  this  intensity  being  further  divided  by  the  small  input 
intensity  of  beam  4  to  give  the  large  value  of  the  ratio  R.  A  similar  description 
is  true  for  maximum  S,  with  the  roles  of  the  input  beams  reversed.  However,  we 
expect  that  phase  mismatch  and  linear  absorption,  which  were  neglected  in 
obtaining  the  analytic  results  of  previous  section,  will  decrease  the  reflectivity 
appreciably.  When  R  is  high,  S  is  correspondingly  low:  the  diffraction  of  a  weak 
beam  is  further  scaled  down  by  the  large  input  intensity  in  that  arm  to  give  the 
ratio  S.  Hence  S  is  appreciable  only  when  R  is  negligible,  much  like  a  seesaw! 
Note  that  the  harder  the  stronger  arm  is  pumped,  the  smaller  is  the  minimum 
intensity  of  the  signal  that  can  be  phase  conjugated  in  the  other  arm. 

Figure  2.12  shows  the  same  effect  on  a  logarithmic  scale,  with  a  wider 
dynamic  range  of  the  phase  conjugator  for  higher  coupling  strengths.  Note  that 
even  a  slight  change  in  the  total  intensity,  and  hence  in  G,  can  change  the 
minimum  value  of  the  weak  beam  intensity  that  can  be  phase  conjugated  by  a 
few  orders  of  magnitude.  This  is  because  the  coupling-strength  product  appears 
in  the  argument  of  an  exponential  function  in  the  transcendental  equation, 


Fig.  2.12.  Phase-conjugate  reflectivities  R  ( - )  and  S  ( - )  in  both  the  arms  versus  the  input 

beam  ratio 
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Fig.  2.13.  Phase-conjugate  reflectivity  R  in  the  left  arm  as  a  function  of  the  interaction  length 


hence  any  slight  variations  in  it  will  cause  a  large  change  in  the  results.  The 
sensitivity  of  the  MPPC  device  to  weak  input  beam  energy  in  one  arm,  and  the 
rapid  decrease  in  reflectivity  in  that  arm  with  increasing  probe  strength,  is  simi¬ 
lar  to  the  experimental  results  recently  reported  by  Lebow  and  Ackerman  [2.51] 
for  phase  conjugation  by  BEFWM. 

Figure  2.13  is  a  plot  of  R  as  a  function  of  the  interaction  length  for  different 
values  of  /40,  keeping  /2L  fixed.  As  expected  from  (2.52)  a  higher  total  input 
intensity  requires  a  smaller  interaction  length  to  start  the  process.  The  interest¬ 
ing  feature  revealed  in  this  graph  is  that  once  the  threshold  interaction  length  is 
crossed,  it  is  not  advantageous  to  use  much  longer  lengths  as  the  curves  of  R 
versus  L  remain  pretty  flat  at  longer  lengths,  this  being  more  so  at  higher  input 
intensities.  The  curves  for  different  /40  cross  because  for  a  fixed  /2L,  a  smaller 
value  of  /40  leads  to  a  higher  phase-conjugate  reflectivity,  with  larger  interaction 
lengths  required  for  the  process  to  occur.  Figure  2.14  is  a  similar  graph  for  the 
phase-conjugate  intensity  reflectivity  S  in  the  other  arm,  the  parameter  values 
being  the  same  as  in  the  set  of  Fig.  2.13.  In  contrast  to  Fig.  2.13,  S  is  higher  for 
larger  /40  because  of  the  increased  diffraction  of  beam  4  to  create  the  phase- 
conjugated  beam  1.  A  larger  value  of  /40  also  requires  a  smaller  interaction 
length  for  observing  the  process,  hence  there  is  no  crossing  of  the  curves  in  Fig. 
2.14  as  is  seen  in  Fig.  2.13. 

The  intensity  distribution  of  the  various  beams  within  the  interaction  re¬ 
gion  is  shown  in  Figs.  2.15-2.18.  The  Kerr  coupling  coefficient  and  the  input 
intensities  of  beams  2  and  4  are  kept  fixed  at  1  cm/MW,  0.5  MW/cm2,  and  0.6 
MW/cm2  respectively.  Figure  2.15  is  for  an  interaction  length  of  1  cm,  and  the 
coupling  strength  corresponding  to  this  length  is  only  slightly  above  the  thresh¬ 
old  value.  Most  of  the  energy  remains  in  the  incident  beams,  leading  to  low 
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Fig.  2.14.  Phase-conjugate  reflectivity  S  in  the  right  arm  as  a  function  of  the  interaction  length 


Fig.  2.15.  Field  distribution  for  MPPC  inside  the  interaction  region  with  L  =  1  cm 


phase-conjugate  reflectivities  in  both  arms.  With  an  increase  in  the  interaction 
length  to  2  cm  (Fig.  2.16),  the  phase-conjugate  beams  emerge  with  more  energies 
at  the  output  planes  than  in  Fig.  2.15,  and  since  the  input  beams  are  not  totally 
depleted,  R  and  S  do  not  reach  their  maximum  value.  For  L  =  4  cm  (Fig.  2.17) 
there  is  a  complete  transfer  of  energies  from  the  incident  beams  to  the  phase- 
conjugate  beams.  Further  increase  in  L  causes  the  region  where  energy  transfer 
occurs  to  shrink  (see  Fig.  2.18  for  L  =  10  cm),  because  now  the  process  does  not 
occur  gradually  throughout  the  interaction  length.  Instead,  the  intensities  of  the 
incident  beams  remain  fairly  constant  for  most  of  the  distance  traversed,  there 
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Fig.  2.16.  Field  distribution  for  MPPC  inside  the  interaction  region  with  L  =  2  cm 


Fig.  2.17.  Field  distribution  for  MPPC  inside  the  interaction  region  with  L  =  4  cm 


being  a  sharp  power  transfer  which  is  confined  to  a  narrow  region  around  the 
center  of  the  interaction  length.  This  behavior  can  be  understood  if  we  recall 
that  beams  1  and  3  are  zero  at  the  z  —  0  and  z  =  L  planes  respectively,  so  that 
there  are  no  gratings  to  begin  with  at  both  the  ends  of  the  interaction  length. 
The  diffracted  beams,  and  hence  the  gratings,  slowly  increase  and  are  apprecia¬ 
ble  near  the  middle,  where  maximum  interaction  occurs.  Complete  depletion  of 
the  incident  beams  in  this  region  once  again  results  in  the  absence  of  the  grat¬ 
ings,  and  hence  the  behavior  shown  in  Fig.  2.18.  An  interesting  outcome  is  that 
for  MPPC  in  which  two  beam-pairs  undergo  acousto-optic  Bragg  scattering  by 
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Fig.  2.18.  Field  distribution  for  MPPC  inside  the  interaction  region  with  L  =  10  cm 

a  common  sound  wave,  a  complete  energy  transfer  occurs  for  a  finite  thickness 
of  the  medium,  in  contrast  to  the  infinite  interaction  region  required  for  100% 
efficiency  for  standard  acousto-optic  Bragg  scattering  of  two  beams. 


2.5  Summary 


We  have  described  nonlinear  optical  phase  conjugation  by  four-wave  mixing  in 
Kerr  media.  Analytic  results  are  obtained  for  undepleted  pumps  and  negligible 
linear  absorption  in  the  medium  both  for  equal  and  different  pump  intensities. 
Pump  depletion  is  accounted  for  in  the  transmission  and  reflection  geometries 
separately.  We  have  also  considered  mutually  pumped  phase  conjugation  in 
electrostrictive  Kerr  media. 
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ABSTRACT 

Volume  phase  gratings  in  photorcfractive  crystals  have  nonuniform  amplitude  and 
phase  due  to  the  energy  exchanged  by  the  writing  beams  within  the  material.  Analytic 
expression  is  obtained  for  the  diffraction  efficiency  of  a  weak  reading  beam  that  does  not 
perturb  the  grating  and  has  a  different  polarization  from  the  write  beams.  For  a  read  beam 
of  arbitrary  intensity,  the  diffraction  efficiency  is  a  nonlinear  function  of  the  read  beam 
intensity  and  is  nonreciprocal  with  respect  to  readout  from  the  two  input  ports.  These 
properties  of  photorefractive  gratings  arc  studied  for  arbitrary  phase  shifts  of  the  index 
grating  from  the  intensity  pattern. 


I,  INTRODUCTION 


The  photorcfractive  crystal  has  proven  to  be  an  almost  ideal  holographic 
medium,  demonstrating  its  ability  to  write  efficient  volume  holograms  at  low  levels  of 
light  without  the  need  for  special  development  procedures.  The  amplitude  and  phase  of 
such  holograms  are  nonuniform1  within  the  crystal  due  to  two-beam  coupling  effects2 
between  the  pair  of  coherent  writing  beams.  Further,  the  read  beam  and  its  diffracted  beam 
induce  a  new  photorcfractive  grating  that  may  be  in  or  out  of  phase  with  the  original 
grating  being  read.3  The  diffraction  is  therefore  quite  different  from  that  of  a  simple  Bragg 
scattering  off  a  fixed  thick  grating,  and  the  coupled  wave  theory  of  Kogclnik  must  be 
modified  to  describe  the  diffraction  efficiency  r|  of  such  gratings.  A  recent  study 
examined  diffraction  efficiency  of  photorcfractive  gratings  when  polarization  asymmetry 
occurs  during  writing  and  readout.  The  new  photorcfractive  grating  formed  by  the  read 
beam  and  its  diffracted  component  is  ignored  for  a  reading  beam  that  is  weak  compared  to 
the  writing  beams.  For  this  case,  r|  is  independent  of  the  read  beam  intensity,  and  the 
argument  of  the  sine  function  in  rj  is  a  nonlinear  function  of  grating  thickness,  coupling 
constants  during  writing  and  readout,  and  input  intensity  ratio  of  writing  beams.  Inclusion 
of  the  effects  arising  from  a  read  beam  having  an  intensity  comparable  to  that  of  the 
writing  beams  has  shown  r\  to  be  a  nonlinear  function  of  the  read  beam  intensity  and 
nonreciprocal  with  respect  to  readout  from  the  two  input  ports.  Similar  nonrcciprocity 
with  respect  to  direction  of  readout  has  been  observed  for  fixed  volume  gratings  in 
photorefractive  media.7  A  theoretical  explanation8  of  this  effect  was  obtained  by  taking 
into  account  the  new  photoinduccd  grating  formed  by  the  read  beam  and  its  diffracted 
component,  besides  the  fixed  grating  present  initially.  The  nonlinear  and  nonreciprocal 
behavior  of  t\  has  been  examined  for  arbitrary  phase  shifts  of  the  photorefractive  grating 
from  the  intensity  pattern.9  Analytic  expressions  for  the  diffraction  efficiency  are  obtained 
for  limiting  values  of  the  phase  shift  between  index  change  and  fringe  pattern. 


IL  THEORY 


TA1  Diffraction  Efficiency  of  Fixed,  Uniform  Gratings 

Periodic  media  with  spatially  varying  conductivity  a  are  known  as  absorption 
gratings,  while  periodic  media  with  spatially  varying  dielectric  constant  e  are  known  as 
phase  gratings;  when  both  a  and  e  are  periodically  varying  functions  in  space,  then  we 
have  a  “mixed”  grating.  An  absorption  grating  is  necessarily  lossy  because  the  nonzero 
conductivity  of  the  medium  leads  to  the  decay  of  an  incident  light  beam.  An  example  of  a 
1-d  phase  grating  with  a  sinusoidal  index  variation  is  given  by 

n(x)  =  nQ  +  nj  cos  Kx  ,  (1) 

where  nQ  is  the  average  index  of  refraction  of  the  medium,  is  the  uniform  amplitude  of 
the  fixed  index  grating,  K  is  the  grating  wave  number  equal  to  2tt/A  (A  being  the  grating 
spacing),  and  x  is  the  spatial  coordinate  along  the  grating  wave  vector.  The  “thickness”  of 
a  grating  is  determined  by  the  dimensionless  parameter  Q  =  2jdL/n0A2,  where  X  is  the 

wavelength  of  light  and  L  is  the  thickness  of  the  grating  in  the  direction  of  light 
propagation.  A  grating  is  planar  or  thin  if  Q  <  1,  and  it  is  a  thick  volume  grating  if 
Q  >  10.  For  grating  spacing  comparable  to  wavelength  (A  =  X),  Q  is  less  than  or  greater 
than  unity  depending  on  whether  the  grating  thickness  L  is  smaller  or  greater  than  A~X. 
If  the  grating  planes  are  perpendicular  (parallel)  to  the  grating  surface,  then  we  have 
unslanted  transmission  (reflection)  gratings.  The  diffracted  beam  appears  on  the  same  side 
of  the  grating  as  the  transmitted  (reflected)  read  beam  for  a  transmission  (reflection) 
grating  (Fig.  1).  Table  I  by  Kogelnik10  summarizes  the  various  diffraction  orders  and  the 
maximum  diffraction  efficiency  that  can  be  obtained  theoretically  from  uniform  fixed 
gratings.  Note  that  the  diffraction  efficiency  of  thin  phase  gratings  is  limited  to  about 
34%  in  the  first  orders  in  the  Raman-Nath  regime.  For  volume  phase  gratings,  maximum 
diffraction  efficiency  is  100%  to  the  Bragg-matched  first  order,  all  higher  orders  being 
suppressed. 

The  fixed  volume  phase  grating  of  Eq.  (1)  is  formed  by  recording  the  interference 
pattern  of  two  coherent  writing  beams  IAW  and  Iuw 

Iw(x)  =  Iow  (1+m  cos  Kx)  ,  (2a) 

where  IQW  =  IAW  +  IfiW  is  the  sum  of  the  intensities  of  the  two  writing  beams, 

m  =  2-\J  IAWIBW  /  (IAW+IBW)  *s  modulation  depth  of  the  intensity  pattern,  and  the 
subscript  W  indicates  writing  beams.  If  p  =  IBW/IAW  IS  intensity  ratio  of  the  writing 

beams,  then  m  =  2V]}  /  (1  +  p).  The  recording  medium  (a  photographic  plate,  for 
example)  is  exposed  to  the  intensity  pattern  of  Eq.  (2a).  After  processing  with  a  developer 
and  a  fixer,  the  plate  is  bleached  to  convert  the  recorded  amplitude  modulation  into  a  phase 
modulation.  Assuming  a  linear  response  of  the  photographic  plate  to  the  incident 
exposure,  the  refractive  index  is  given  by 
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Fig.  1.  Diffraction  of  a  read  beam  from  various  types  of  fixed,  uniform  gratings. 
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Diffraction  Properties  of  Various  Gratings 


^  2nXL 

Type 

Geometry 

Diffraction 
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Bmax 
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Transmission 

Two:  ±1 

6.25% 

Q  <  i 

a(x)=a(x+A) 

Thin 

Phase 

Transmission 
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33.9% 

n(x)=n(x+A) 

lst-order 

Absorption 

Transmission 

One 

3.7% 

Q>  1 

a(x)=a(x+A) 

Reflection 

Bragg- 

7.2% 

Thick 

Phase 

Transmission 

Matched 

100% 

n(x)=n(x+A) 

Reflection 

Order 

100% 

n(x)  =  n0  +  mo  m  cos  Kx  =  n0  +  2mo  — cos  Kx 

1  +  p 


(2b) 


% 
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where  n10  is  a  constant  of  proportionality  that  depends  on  the  incident  exposure  and  the 
recording  medium. 

To  derive  Kogelnik’s  formula4  for  diffraction  efficiency  of  volume  phase  gratings 
with  uniform  amplitude,  consider  the  transmission  grating  shown  in  Fig.  2  with  an  index 
of  refraction  given  by  Eq.  (1).  The  electric  field  of  the  incident  read  beam  and  its  diffracted 
component  are  described  as  plane  waves 

EarO-.O  =  \  Ar(z)  cxp[i(kAR*r  -  cot)]  +  c.c.  ,  (3) 


where  AR  is  the  complex  amplitude  of  the  read  beam  at  steady-state,  co  is  the  frequency 
and  kAR  is  the  corresponding  wave  vector  in  the  medium  at  an  angle  0AR  with  respect  to 
the  z-axis.  To  describe  the  diffracted  beam,  A  is  replaced  by  B.  The  subscript  R  denotes 
the  read  beam  and  its  coherently  generated  diffracted  beam.  When  the  electric  field  of  the 
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Fig.  2.  The  schematic  representation  of  diffraction  from  a  uniform  phase  grating. 


optical  wave  is  described  by  Eq.  (3),  then  the  intensity  is  equal  to  e0n0cIAR(z)i2/2  in 
mks  units,11  where  e0  is  the  electric  permittivity  of  free  space  and  c  is  the  velocity  of 
light.  Ignoring  the  multiplicative  constant  factor,  Iar(z)=IAr(z)I2/2  is  taken  to  represent 
the  intensity  of  the  beam  in  the  rest  of  this  paper.  The  z  axis  is  taken  normal  to  the 
surface  of  the  medium,  and  the  complex  amplitudes  arc  assumed  to  be  functions  of  z  due 
to  absorption  and  diffraction.  The  amplitudes  arc  independent  of  the  transverse  coordinate 
x  because  the  grating  is  assumed  to  be  infinite  in  that  direction.  To  describe  the 
propagation  and  coupling  of  the  two  waves  at  steady  state  in  the  thick  transmission  phase 
grating,  we  start  with  the  scalar  wave  equation 


v2e  + 


2  2 
n  0) 


E  =  0 


(4) 


where  E  =  EAR  +  EBR  is  the  total  electric  field.  Substituting  for  n  and  E  in  the  wave 
equation,  using  the  slowly  varying  envelope  approximation  and  ignoring  phase* 
mismatched  terms  and  terms  that  are  quadratic  in  the  index  change,  we  arrive  at  the 
following  coupled  wave  equations 
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(5a) 


where 


Akz  =  k  (cos0BR  -  cos0AR)  ,  k  sin0BR+  K  =  k  sin0AR 


(5b) 


aR  is  the  intensity  absorption  coefficient  of  the  grating  at  the  reading  wavelength  that  has 
been  introduced  phenomenologically  to  the  coupled  wave  equations  and  Akz  is  the  phase- 

mismatch  along  the  z-dircction.  Integration  of  the  coupled  wave  equations  (5a)  over  the 
transverse  coordinate  x  lead  to  the  last  relation  expressed  in  Eq.  (5b).  For  phase-matched 
readout,  Akz  =  0  (Fig.  3).  Ignoring  the  trivial  solution  of  0BR  =  0AR  which  corresponds 

to  the  diffracted  beam  propagating  along  the  read  beam,  the  relevant  solution  is 
0ar  =  -  0BR  =  0R  for  Akz  =  0  (Bragg-matched  readout),  with  0R=  sin'^R^noA). 

Hence  the  Bragg-matched  diffracted  beam  propagates  at  an  angle  to  the  grating  planes  that 
is  equal  and  opposite  to  the  reading  beam  angle.  Using  boundary  conditions  AR(0)  =  ARo 

and  Br(0)  =  0,  the  solutions  for  the  beam  amplitudes  arc  given  by 


Ar(z)  =  ARo  cxp(-  )  C0SK2  , 

oc  z 

Br(z)  =  i  ARo  cxp(  -  2pos0-  )  sinKZ’  (6) 

where  k  =  jh^ArCOsOr,  0R  being  the  Bragg  angle  for  the  reading  beam.  Figure  4  shows 
the  normalized  beam  intensities  I Ar(z)I2/I ARqI2  (solid  curve)  and  IBr(z)!2/I ARqI2  (dashed 
curve)  as  a  function  of  distance  z  into  a  lossless  grating  (aR  =  0).  At  the  beginning  of  the 
thick  phase  grating,  there  is  no  diffracted  beam,  and  the  normalized  read  beam  intensity  is 
unity.  Periodic  exchange  of  energy  takes  place  between  the  read  beam  and  its  diffracted 
component.  The  ratio  of  the  diffracted  beam  intensity  at  the  output  of  the  grating  to  the 
read  beam  intensity  at  the  input  of  the  grating  determines  the  diffraction  efficiency: 
T)  =  IBr(L)I2  /  IARQr.  Substituting  for  k,  wc  arrive  at  Kogclnik^s  formula4 


T[  =  exp(- 


cos0R 


o  /  Ttn  i  L  v 

)  S*n  ™cA  ) 


ARCOS0R 


(7) 


The  argument  of  the  sine  function  depends  on  the  product  of  grating  amplitude  times  its 
thickness  and  on  the  wavelength  and  Bragg  angle  at  readout,  but  the  argument  is 
independent  of  the  read  beam  intensity.  For  a  lossless  grating  (aR  =  0)  and  Bragg-matched 
readout  (Akz  =  0),  diffraction  efficiency  is  equal  to  unity  whenever  the  argument  of  the 

sine  function  is  an  integral  multiple  of  rc/2,  implying  total  transfer  of  the  read  beam 
energy  to  the  diffracted  beam.  Substituting  Eq.  (2b)  into  Eq.  (7),  wc  obtain 


Fig.  4.  Intensities  of  the  read  beam  (solid  curve)  and  the  diffracted  beam  (dashed  curve)  as 
a  function  of  normalized  depth  kz  within  the  medium. 


tj  =  exp(- 


3*Luin2  r  ^  2^3. 

COS0R  LXRCOS0R  1+p 


]. 


(8) 


The  above  equation  illustrates  the  variation  of  the  diffraction  efficiency  of  uniform  phase 
grating  with  the  intensity  ratio  of  the  writing  beams. 


TB1.  Photorcfractive  Gratings  with  Nonuniform  Amplitude  and  Phase 

Two  coherent  light  beams  intersect  in  a  photorcfractive  crystal  when  writing  a 
grating  (Fig.  5);  due  to  photorcfractive  two-beam  coupling,  one  writing  beam  gains 
energy  at  the  expense  of  the  other  writing  beam.  The  nonreciprocal  energy  transfer 
between  the  two  beams  is  due  to  Bragg  scattering  off  the  grating  the  beams  themselves 
wrote.  This  phenomenon  is  known  as  “Self-Diffraction”  ,  and  the  direction  of  energy 
transfer  is  determined  by  the  orientation  of  the  crystal  and  the  nature  of  charge  carriers 
(electrons  or  holes).  Such  an  energy  exchange  is  not  found  in  most  nonlinear  materials 
because  they  respond  locally  to  optical  beams,  with  no  phase  shift  between  index  grating 
and  intensity  pattern.  We  will  now  show  that  the  amplitude  and  phase  of  the  index 
modulation  are  no  longer  uniform  due  to  the  exchange  of  energy  between  the  two  writing 
beams  in  the  photorcfractive  medium.  Similar  to  the  previous  section,  we  describe  the 
two  writing  beams  as  plane  waves  with  complex  amplitudes  Aw  and  Bw.  These  two 
beams  interfere  to  produce  a  sinusoidal  intensity  pattern  given  by 
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Fig.  5.  Schematics  of  photorefractive  two-beam  coupling. 


Iw(x,z)  =  Iow(z)  (l+  c'Kx  + 


m*(z) 


-iKx> 


(9) 


where  Iow(z)  is  the  local  sum  of  the  two  beam  intensities:  Iow(z)  ~  lAw(z)+W(z)  = 

[  IAw(z)l2  +  IBw(z)l2  ]/2.  K  is  the  magnitude  of  the  grating  wavevector:  IKI  =  lkAW  - 

kBWl  =  2tc/A,  A  being  the  fringe  separation  equal  to  Xw/(2sin9WE)  when  the  two  beams 

are  symmetrically  incident  at  an  external  angle  of  20WE  about  the  crystal  face  normal. 

is  the  wavelength  of  the  two  writing  beams  and  x  is  the  direction  along  the  grating  wdve 

vector.  m(z)  is  the  local  modulation  depth  of  the  incident  intensity  pattern:  m(z)  = 
* 

Aw(z)Bw(z)/I0(z);  it  depends  on  the  polarizations  and  intensity  ratio  of  the  two  writing 
beams.  The  high  intensity  regions  of  the  spatially  periodic  intensity  pattern  cause 
photoionization  of  the  impurity  levels  in  the  bandgap.  The  mobile  charges  drift  and 
diffuse  before  recombining  and  being  trapped  in  regions  of  low  intensity,  leading  to  a 
separation  of  ionic  charges  inside  the  crystal.  The  resulting  space-charge  Held  produces  a 
photorefractive  index  grating  via  the  linear  elcctro-optic  effect.  Assuming  electrons  to 
be  the  charge  carriers,  the  photorefractive  index  grating  can  be  written  as 


n  =  nb  + 


AWBW 


iw  f  i<> 

2  1  IAWIZ  +  lBu,lz 


elKx  +  c.c.] 


'W' 


(10a) 


where 


n.,w  “  Hi*  r  .. ...  E, 


Eo>Hd 


iw  “  "b  ‘crf,w  ^  \  E0+(Ed+Eq) 


tan<|>  = 


Ep2+Ed(Ed+E0) 


Here  nb  is  the  background  refractive  index  of  the  pholorcfractive  medium  in  the  absence  of 
light,  niw  is  the  amplitude  of  the  index  grating  at  the  fundamental  harmonic  of  the  grating 

vector  and  <J>  represents  the  spatial  phase  shift  of  the  index  gratings  with  respect  to  the 
light  intensity  pattern.  rcffW  is  an  effective  electro-optic  coefficient  that  depends  on  the 

polarization  of  the  writing  beams,  the  orientation  of  the  grating  wave  vector  that  is 
written,  and  the  non-zero  Pockcls  coefficients  of  the  pholorcfractive  medium.  E0  is  the  dc 
electric  field  applied  to  the  pholorcfractive  crystal,  Ed  is  the  characteristic  electric  field  due 
to  diffusion,  and  Eq  is  the  saturation  value  of  the  space  charge  field.  Substituting  for  n 
and  E  in  the  scalar  wave  Eq.  (4),  using  the  slowly  varying  envelope  approximation  and 
ignoring  phase-mismatched  terms  and  terms  that  arc  quadratic  in  the  index  change,  we 
arrive  at  the  following  coupled  wave  equations 
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where  aw  is  the  intensity  absorption  coefficient  at  the  writing  wavelcgth  and  0WI  is  the 
half  angle  between  the  beams  inside  the  crystal.  Expressing  the  complex  beam  amplitudes 

in  polar  form:  Aw  =  ^2 IAW  exp(i\|/AW),  we  obtain  the  following  equations  for  beam 
intensities  and  phases 
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The  intensity  equations  and  subsequently  the  phase  equations  can  be  solved  as  functions 
of  distance  z  in  the  crystal 

W2)  =  W°)  1  +  P  exP(  “  COS0W,  )  ’ 

!bw(z)  =  W°)  P  +  e-^  CXP(  “  cose^  )  •  (,3a) 

Vaw(z)  “  Vbw(z)  =  Vaw(°)  -  Vii\v(°)  “  Cz  +  ln[^p“]  .  (!3b) 

where  p  =  IbwC^awCO)  is  the  intensity  ratio  of  die  writing  beams  at  the  input.  We  find 
that  the  intensity  of  each  beam  at  a  distance  z  depends  on  the  input  intensities  of  both 
the  beams,  and  there  is  an  exponential  dependence  on  the  product  yz.  If  the  beam  being 
amplified  (=  probe  or  writing  beam  B)  is  initially  much  weaker  than  the  beam  being 
«  depleted  (=  pump  or  writing  beam  A)  such  that  P«l,  and  p«exp(-yL),  we  obtain  an 

exponential  buildup  for  the  probe;  while  the  pump  is  undcplcted  and  remains  constant  at 
its  input  value  in  a  lossless  photorcfractivc  medium:  IliW(L)  ~  IBW(0)  cyL  and  IAW(L)  = 
IAW(0).  Substituting  the  solutions  for  the  beam  intensities  and  phases  given  by  Eqs.  (13) 
into  the  expression  for  the  photorcfractivc  index  grating  given  by  Eq.  (10a),  we  arrive  at 
the  following  expression 

n(x,z)  =  nb  + - n'W  /2  cos[Kx  +  <J>  +  A\|/(0)  -  Cz  +  ln(^rt~)] 

Vp  d1*12  +  °—j=r  7  ^ 

Vp 


=  nb  +  nicff(z)  cos[Kx+(J)+A\)/cff(z)]  (14) 


where  A \j/  =  \|/AW  -  VBW-  This  is  the  general  expression  for  a  photorefractive  index 
grating  using  transmission  geometry  that  is  valid  for  arbitrary  phase  shift  4>  of  index 
grating  from  the  intensity  pattern.  In  contrast  to  the  fixed,  volume  phase  grating  with 
uniform  amplitude  and  phase  of  Eq.  (1),  the  photorefractive  index  grating  in  the  above 
equation  has  an  amplitude  and  phase  that  vary  with  the  depth  in  the  grating.  Bending  of 
the  photorefractive  hologram  through  the  crystal  thickness  is  caused  by  A\j/Cff(z)  -  there  is 
no  bending  for  $  =  tc/2,  and  bending  increases  as  <j)  shifts  from  n/2.  Note  that  Eq.  (14)  is 
independent  of  absorption  losses  (<*  aw)  in  the  medium  at  the  writing  wavelength. 
Without  loss  of  generality,  we  arbitrarily  set  the  relative  phase  between  the  two  writing 
beams  equal  to  zero  at  the  input  (A\(/(0)  =  0).  Figure  6a  is  a  plot  of  the  nonuniform 
grating  amplitude  2nicff(z)/nj  w  and  Fig.  6b  is  a  plot  of  the  nonuniform  phase  of  the 


grating  Aycff(z)  as  a  function  of  depth  z  into  the  grating.  The  phase  shift  <t>  of  the  index 
grating  from  the  intensity  pattern  is  0  (dashed  curve),  n/6  (dotted-dashed  curve)  and  jc/2 
(solid  curve),  while  (3  =  0.01  and  2rcniw  /  (^.\yC0S^Wl)  =  20  em'f .  Note  that  A\yeff(z) 
changes  sign  in  the  bulk  of  the  grating  for  <j)  =  n/6,  so  that  the  surfaces  of  equal  index 
change  are  curved  and  tilted.^ 

For  <[>  =  0,  there  is  no  phase  shift  of  the  index  grating  from  the  intensity  pattern, 
so  that  the  intensity  coupling  constant  y  between  the  writing  beams  is  equal  to  zero,  and 
the  phase  coupling  constant  C  is  equal  to  its  maximum  value  of  Cm  =  nn^wCosOw,.  In 
this  case,  Eq.  (14)  simplifies  to 


n(x,z) 


VI 

nb  +  niw  j+p 


COS(Kx-Cmz) 


(15a) 


which  is  a  grating  with  uniform  amplitude  and  a  spatial  frequency  that  is  modified  by  the 
phase  coupling  between  the  writing  beams.  The  phase  of  the  grating  varies  linearly  with 
depth  into  the  crystal,  so  that  the  grating  planes  are  rotated  to  form  slanted  gratings. 
Maximum  diffraction  efficiency  will  occur  at  a  read  angle  different  from  that  during 
writing.  Comparing  with  Eq.  (2b),  the  amplitude  niw  of  the  photorcfractive  index  grating 
given  in  Eq.  (15a)  must  be  replaced  by  2n10  to  yield  the  results  of  the  uniform  phase 
grating. 


Distance  z  (cm) 

Fig.  6a.  Nonuniform  amplitude  2nicff(z)/ni\y  of  photorcfractive  grating  for  different 
values  of  phase  shift  <J)  of  index  grating  from  intensity  pattern.  The  beam  ratio  P  during 
writing  is  equal  to  0.01,  and  2nni\yAwcos0wi)  =  20  cm"1. 
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Fig.  6b.  Nonuniform  phase  A\|/Cff(z)  of  photorcfractivc  grating  for  different  values  of 
phase  shift  <)>  of  index  grating  from  intensity  pattern.  The  other  parameters  arc  identical  to 
Fig.  6a. 


Similarly,  if  0  —  rt/2,  the  phase  coupling  constant  is  equal  to  zero  and  y  is 
equal  to  its  maximum  value  of  V  =  iTnqAwCOsGwi-  In  this  case,  Eq.  (14)  simplifies  to 

n(x,z)  =  nb  -  niw  - - - ^  sinKx  .  (15b) 

Vpcrz/2  +  C-^ 

Hence  the  photorcfractivc  index  grating  in  the  above  equation  has  a  uniform  phase  and  an 
amplitude  that  varies  with  the  depth  in  the  grating  and  is  spatially  shifted  by  tt/2  from  the 
intensity  pattern  of  Eq.  (9).  The  nonuniformity  2nicff(z)/niw  of  the  grating  amplitude  is 
a  function  of  the  input  beam  ratio  P  and  the  strength  of  the  beam  coupling  V  and  is 
illustrated  in  Fig.  7  for  r=20cm"1,  a  coupling  value  easily  achieved  in  single  crystal 
BaTiC>3.  The  grating  amplitude  is  maximum  when  the  two  writing  beams  have  equal 
intensity;  for  P  =  1  ease,  this  occurs  at  the  input  plane  itself,  while  for  the  P  =  0.01  case, 
the  two  writing  beams  have  equal  intensity  only  after  they  have  propagated  some  distance 
into  the  crystal  and  the  weak  probe  has  been  amplified  at  the  expense  of  the  strong  pump 
beam  via  photorefractive  two-beam  coupling.  For  p  =  100,  the  probe  is  stronger  than  the 
pump  to  begin  with,  and  two-beam  coupling  effects  serve  to  further  amplify  the  probe  in 
the  medium.  Hence  the  optimum  beam  ratio  is  never  achieved  and  the  grating  amplitude 
decreases  monotonically  in  the  medium.  It  is  important  to  study  the  effects  of  this 
departure  from  uniformity  of  photorcfractivc  index  gratings  on  the  diffraction  efficiency  of 
a  weak  reading  beam.  As  we  shall  show  in  the  next  section,  maximum  diffraction 
efficiency  occurs  for  p  =  0.01,  when  the  grating  peaks  near  the  middle  of  the  interaction 
length.12 
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Fig.  7.  Nonuniform  amplitude  2nldy(z)/n  of  photorcfractivc  grating  for  different  values 
of  beam  ratio  p  during  writing.  The  two-beam  coupling  constant  T  is  equal  to  20  cm*1. 


fCI.  Readout  of  n/2  Shifted  Photorcfractivc  Gratings  with  a  Weak  Incoherent  Beam 

The  basic  interaction  geometry  is  illustrated  in  Fig.  8.  Aw  and  Bvvare  the  two 
coherent  writing  beams  of  wavelength  that  intersect  within  a  photorcfractive  crystal  of 
thickness  L.  The  angle  between  the  two  beams  is  20WI  in  the  medium.  The 
photorcfractivc  grating  induced  by  this  pair  of  beams  is  read  by  an  incoherent  third  beam 
Ar  of  wavelength  XR  that  is  incident  at  a  different  angle  (0RI  in  the  medium)  for  Bragg 
matched  readout.  A  fourth  beam  BR  is  generated  by  the  diffraction  of  the  read  beam  off 
this  grating.  The  theoretical  problem  of  interest  is  to  find  an  expression  for  the  diffraction 
efficiency  of  the  photorcfractive  index  grating  when  (J)  =  ti/2,  as  given  by  Eq.  (15b).  The 
intensity  of  the  read  beam  IAR  is  assumed  to  be  weaker  than  the  writing  beam  intensities 

IAW,  W> so  l^at  readout  process  docs  not  perturb  the  grating  being  read.  Also,  Bragg 
condition  is  satisfied  during  readout,  so  that  k AR  -  kBR  =  kAW  -  kBW  =  K.  By 
following  the  methods  used  in  the  previous  sections,  we  arrive  at  the  coupled  wave 
equations  for  the  amplitudes  of  the  read  beam  AR  and  the  diffracted  beam  BR 
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Fig.  8.  Schematics  of  simultaneous  writing  and  readout  of  photorefractive  index  grating. 
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Since  the  read  beam  and  its  diffracted  component  can  have  a  polarization  different  from 
that  of  the  writing  beams,  rcff  at  readout  can  be  different  from  that  during  writing.  Such  a 

possibility  is  allowed  for  by  distinguishing  between  the  grating  amplitudes  nx  during 

writing  and  readout,  even  though  both  the  processes  involve  the  same  space  charge 
grating  with  wave  vector  K.  Using  boundary  conditions  AR(0)  =  Aro>  Br(0)  =  0,  we 

obtain  the  following  expression  for  the  diffraction  efficiency  (=  1BR(L)/ARQ12)  of  the 
hologram5 


» 


T[  =  exp^-  ^^jsin2  ^(tan'Rc^Vp)  -  tan'RVp))] 


(17) 


where  k  =  rcn1RARcos0RI  is  the  coupling  constant  during  readout,  k  and 
T  =  2rcnlw/Xwcos0WI  arc  related  but  can  be  quite  different  in  value  depending  on  the 
polarization  states  of  the  writing  beams  and  that  of  the  reading  beam.  Note  that  k  (r)  is 
the  amplitude  (intensity)  coupling  constant,  so  that  2k  =  F  for  identical  polarization  and 
wavelength  during  writing  and  readout.  Then  the  expression  (17)  for  diffraction  efficiency 
simplifies  to 


ocrL  ]  P  [exp(rL/2)  -  l]2  , 
cosGri  1 1+P  '  1  +  P  cxp(fL) 


(18) 


which  is  identical  to  the  result  obtained  by  Kukhtarev  et  al12  for  a  lossless  medium 
(ccr  =  0).  In  the  limit  of  no  coupling  during  writing  (r—>0,  k  remains  finite),  Eq.  (17) 
reduces  to  the  Kogelnik  result  given  by  Eq.  (8)  if  we  use  the  fact  that  n1R  is  equal  to 
2n10  of  the  fixed,  uniform  grating;  this  limit  can  be  achieved  in  noncubic  photorefractive 

crystals  by  writing  a  grating  that  can  be  read  out  only  by  anisotropic  diffraction.  In 
particular,  we  examine  the  diffraction  efficiency  of  holograms  written  in  BaTiC>3  for  the 
following  two  cases.  In  the  first  case,  the  writing  beams  at  =  514.5  nm  are 
extraordinary-polarized,  resulting  in  strong  two-beam  coupling  (rc=20cnr1).  In  such  a 
case,  the  index  modulation  varies  significantly  over  the  thickness  of  the  crystal  because 
one  beam  is  strongly  amplified  as  the  other  strongly  depletes.  For  the  second  case,  the 
writing  beams  arc  ordinary-polarized  for  which  the  beam  coupling  constant  is  much 
smaller  (r0=lcm'1).  In  both  cases,  the  grating  is  read  out  by  a  weak  beam  at 
XR=  633  nm  that  is  extraordinary-polarizcd  (K=rc/2=10cm'1).  Plots  of  the  diffraction 
efficiency  r|  as  a  function  of  medium  thickness  L,  as  given  by  Eq.  (17),  are  shown  in 
Fig.  9,  where  the  beam  ratio  is  taken  to  be  0.1  (absorption  is  neglected,  for  simplicity). 
Also  plotted  for  comparison  is  the  diffraction  efficiency  of  a  uniform  hologram  (with  no 
coupling,  corresponding  to  Eq.  (8)),  with  the  same  reading  beam  coupling  constant  k. 
Note  that  for  the  case  of  polarization  asymmetry  during  writing/rcading,  the  argument  of 
the  sine  function  in  Eq.  (17)  can  be  quite  large,  allowing  q  to  exhibit  one  complete  cycle. 
Figure  10  is  a  plot  of  the  diffraction  efficiency  q  as  a  function  of  beam  ratio  P  for 
L  =  0.5  cm.  The  values  of  the  coupling  constants  arc  the  same  as  in  Fig.  9.  Note  that 
writing  with  ordinary-polarized  beams  gives  results  similar  to  the  Kogelnik  case  due  to 
the  smaller  two-beam  coupling  effects.  Corresponding  to  this  choice  of  parameter  values, 
the  maximum  diffraction  efficiency  occurs  at  an  input  intensity  ratio  of  0.01  for 
extraordinary-polarizcd  writing  beams,  while  the  maxima  occur  at  P  -0.1  and  8  for 
ordinary-polarized  writing  beams. 
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Fig.  9.  Diffraction  efficiency  r|  as  a  function  of  medium  thickness  L.  Input  beam  ratio  P 
is  equal  to  0.1,  and  absorption  effects  arc  neglected  (a  =  0). 

(i)  T  -  2k  =  20  cm’1  (E-write  /  E-rcad  curve) 

(ii)  T=  1  cm  1 , 2k  =  20  cm  1  (O-writc  /  E-rcad  curve) 

(iii)  T=  0,  2k  =  20  cm'1  (Kogclnik  ease) 
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Fig.  10.  Diffraction  efficiency  q  as  a  function  of  beam  ratio  P  during  writing.  Grating 
thickness  L  is  equal  to  0.5  cm.  All  other  parameters  arc  the  same  as  in  Fig.  9. 


IDl.  Readout  of  Photorcfractivc  Gratings  with  a  Strong  Incoherent  Beam 


r4 


The  basic  interaction  geometry  is  identical  to  that  shown  in  Fig.  8,  and  the 
incoherent  read  beam  intensity  I  R  can  be  comparable  to  or  larger  than  the  intensities  of 

the  writing  beams  IAW,  IBW.  Hence  the  new  grating  formed  by  the  read  beam  and  its 
coherently  generated  diffracted  beam  must  be  taken  into  account,  besides  the  grating  being 
read.  In  this  section,  the  diffraction  efficiency  of  a  photorcfractivc  grating  is  examined  for 
arbitrary  phase  shifts  of  the  index  grating  from  the  intensity  pattern.  All  the  four  waves 
are  assumed  to  be  plane  waves  with  identical  linear  polarization,  so  that  the  grating 
amplitude  is  equal  for  the  writing  and  readout  beams.  The  total  intensity  is  the  sum  of 

the  intensity  of  all  the  four  beams  and  the  interference  terms  arising  from  the  two  pairs  of 
coherent  beams 

1(r)  =  2  2  {IAi'2  +  IB/  +  [  Aj  cxP(iKj‘r)  +  c-c-  •  O9) 

J 


where  the  grating  vectors  arc  defined  by  K ■  =  kAj  -  kUj,  j  =  W,  R.  If  Bragg  condition  is 
satisfied  for  the  two  coherent  beam  pairs:  Kw  =  KR  =  K,  then  the  gratings  formed  by  the 
two  beam  pairs  combine  to  form  a  single  grating  that  couples  all  the  beams 


n  = 


[c* 


AWBW 


iawi2+  ibw 


ar  br 


+  ia/+  ib/ 


cxp(iK*r  )  +  c.c.]  . 


(20) 


We  can  derive  the  following  coupled  wave  equations  by  using  the  methods  developed  in 
the  previous  sections 


dAj  AwBw +  ar  br  o 

dz  “  7  IAWIZ+  IBwf+  IA/+  IBrIz  i 

*  AwB  w  +  ArBr 

(k  "  7  IAWIZ+  1BwI’2+  IAr1z+  IBrIz  } 

where  j  =  W,  R,  and  y  is  the  complex  field  coupling  coefficient  defined  as 
y=  iicnje^/C^cosO!).  In  arriving  at  these  equations,  we  have  ignored  the  difference 
>  between  the  wavelengths  (Xw  «  XR  =  X)  and  angles  (0WI  =  0RI  =  0^  of  the  writing  and 

readout  beams,  so  that  the  coupling  constant  is  the  same  for  both  pairs  of  beams.  For 
Bragg-matched  readout,  the  frequencies  of  the  various  beam  pairs  can  differ  by  a  few  tens 
of  gigahertz  and  their  angles  differ  by  a  few  milliradians  for  typical  values  of  1  Jim 
grating  period  and  1mm  thickness  of  the  photorcfractivc  grating.  The  coupled  wave 
equations  (21)  can  be  solved  by  using  methods  developed  in  Refs.  13-15:  the  constants  of 


integration  associated  with  these  equations  are  identified;  two  new  variables  f  =  AR/AW 
and  g  =  BR/BW  are  defined  and  their  equations  arc  solved.  The  intensities  of  the  four 
beams  can  then  be  expressed  as  the  following  functions  of  z9 


IAr(z)I2 


Cw  lf(z)g(z)l2-cRlf(z)l2 
lg(z)l2  -  lf(z)l2 


IBr(z)I2 


cR  lg(z)l2-cw  lg(z)f(z)l2 
lg(z)l2-  lf(z)l2 


IAw(z)l2 


_  CW  l§(Z)^  CR 

■  lg(z)l2-  lf(z)l2 


IBw(z)l2 


cR  cw  lf(z)l 
lg(z)l2-lf(z)l2 


(22) 


where 


f(z) 


* 

sf0  +  (af0  +  2b)  tanh  kz 
s  -  (a  -  2bf0)  tanhK'z 


g(z) 


sgp  ~  (ggp  +  2  b)  tanh  k*  z 
s  +  (a  -  2bg0)  tanh  k*z 


(23) 


and  cR  =  IAr(0)I2+  IBr(0)I2,  cw  =  IAw(0)l2+  IBw(0)l2,  a  =  cR-cw>  b  =  AR(0)Aw(0)  + 

BR(0)iy0),  s2=  a2+  41b!2.  k  =  sy/4I0,  f0  =  Ar(0)/Aw(0),  g0  =  BR(0)/Bw(0),  I0  = 
[IAr(0)I2+  IBr(0)I2+IAw(0)I2+IBw(0)I2]/2.  The  above  equations  arc  general  solutions  for 
forward  four-wave  mixing.  For  the  problem  at  hand,  we  have  the  boundary  conditions 
Br(0)  =  0,  and  AR(0),  Aw(0),  and  Bw(0)  arc  known  quantities.  The  diffraction  efficiency  is 
then  defined  as  T|  h  IBr(L)12/  IAr(0)I2,  where  IBR(L)I2  can  be  calculated  by  evaluating  the 
second  equation  of  (22)  at  z  =  L.  From  the  expressions  for  f(L)  and  g(L),  it  is  obvious 
that  T|  will  be  a  function  of  tanh  kL,  where  k  is  a  complex  quantity  for  any  general  value 
of  the  spatial  phase  shift  (]).  Using  mathematical  identities,  tanh  kL  can  be  expressed  in 
terms  of  tanh  and  tan  ^L,  where  the  subscripts  r  and  i  refer  to  the  real  and  imaginary 
parts  of  a  complex  quantity,  respectively.  Hence  we  can  expect  T|  to  exhibit  both 
oscillatory  and  exponential  behavior  as  a  function  of  Ik!L. 


fi 


[i]  4)  =  n  /  2 

y  is  a  real,  negative  quantity  equal  to  -  yr,  with  yr=  jrn^XcosGj).  The 
expression  forq  simplifies  to9 

16  IAW(0)  IBW(0)  taph  (OAs 

11  “  sz(l  +  tanh\L)  -4s  tanh  KrL  [IAR(0)  +  IAW(0)  -  IBW(0)]  * 


where  =  sy/4I0,  IAj(0)  =  IAj(0)I2/2,  etc.  The  expression  for  q  involves  only  hyperbolic 
functions  of  k^,  and  is  asymmetric  with  respect  to  change  in  sign  of  Kr,  which  implies 
unequal  diffraction  efficiencies  for  readout  from  the  two  input  ports.  A  change  in  the  sign 
of  the  coupling  constant  is  equivalent  to  readout  from  the  other  input  port.  This  is  true 
for  both  equal  and  different  intensity  writing  beams,  except  that  the  latter  case  would 
involve  an  interchange  of  the  unequal  intensities  of  the  writing  beams  besides  the  change 
in  read  beam  intensity  to  describe  readout  from  the  other  input  port. 


For  a  weak  reading  beam  [IAR(0)  «  1aw(0)»  IBW(0)],  the  input  intensities 
dependent  factor  s/(2I0)  that  is  present  in  Kr  is  equal  to  unity,  so  that  Kr  =  y J2.  Hence 
Eq.  (24)  simplifies  to 


_  _g_  [  cxp(yrL)  -  1]  2 
^  -  1+p  1  +  P  cxp(2yrL) 


(25) 


where  p  is  the  input  beam  ratio  of  the  writing  beams  defined  as:  P  =  IBW(0)  /  IAW(0).  This 
result  is  identical  to  the  diffraction  efficiency  obtained  in  Eq.  (18)  if  we  neglect  absorption 
and  use  T  =  2yr.  The  same  expression  for  q  was  obtained  by  Kukhtarcv  el  al:1  for  a 
weak  read  beam  and  n/2  shifted  photorcfractivc  grating.  For  small  coupling  strengths, 
y^  «  1,  and  Eq.  (25)  reduces  to  q  ~  Pyr2L2  /  [(1+p)2].  This  expression  is  identical  to 
Kogelnik’s  result  (Eq.  (8))  for  small  arguments  of  the  sine  function  with  nj  =  2n10. 


Figure  11  is  a  plot  of  the  diffraction  efficiency  as  a  function  of  the  read  beam 
intensity,  as  given  by  Eq.  (24).  The  equal  intensity  of  the  write  beams  is  taken  to  be 
unity:  IAW(0)  =  IBW(0)  =  1.  Note  that  the  modulation  depth  can  still  be  smaller  than 
unity  due  to  the  read  beam  intensity  in  the  denominator  of  the  expression  for  m.  For 
small  coupling  strengths  (yrL  ~  1)  and  weak  read  beam  (IAR(0)  < 
diffraction  efficiency  is  independent  of  the  read  beam  intensity.  Increase  of  the  read  beam 
+'  intensity  above  the  write  beam  intensities  serves  to  erase  the  grating  being  read,  and  the 

diffraction  efficiency  falls  to  zero.  For  strong  coupling,  the  behavior  of  q  is  different — the 
diffraction  efficiency  is  -50%  for  weak  read  beams,  and  increases  as  the  read  beam 
becomes  stronger  than  the  write  beams,  reaches  a  broad  maximum  and  eventually  falls  to 
*  zero  when  IAR(0)  »  IAW(0),  IBW(0).  Enhancement  of  the  diffraction  efficiency  occurs 
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Fig.  11.  Diffraction  efficiency  r\  as  a  function  of  read  beam  intensity  IAR(0)  for  various 
coupling  strengths  yL.  The  normalized,  equal  intensities  of  the  write  beams  are  unity: 

IAW(°)  =  IBW(0)S=  I-" 


when  the  new  grating  formed  by  the  read  beam  and  its  diffracted  component  is  in  phase 
with  the  grating  being  read.  Figure  12  is  the  same  set  of  curves  as  in  Fig.  11  but  for 
negative  coupling  constant  yr,  which  corresponds  to  readout  from  the  other  input  port.  In 
this  case,  the  two  gratings  are  out  of  phase,  so  that  there  is  no  enhancement  of  the 
diffraction  efficiency  for  any  intensity  of  the  read  beam.  As  predicted  by  Eq.  (24),  the 
maximum  diffraction  efficiency  is  only  50%  at  large  coupling  strengths.  Figure  13  is  a 
plot  of  the  diffraction  efficiency  T|  versus  the  normalized  coupling  strength  yju  for 
negligible  absorption  and  equal  intensity  of  write  beams,  IAW(0)  =  IBW(0)  =  1.  The 
asymmetry  in  T|  is  negligible  for  a  weak  read  beam  and  becomes  appreciable  only  when 

Iar(0)  >  W(0)>  W(Q)* 
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Fig.  12.  Diffraction  efficiency  rj  as  a  function  of  read  beam  intensity  IAR(0)  for  various 
negative  coupling  strengths  yL.  The  normalized,  equal  intensities  of  the  write  beams  are 
unity:  IAy/(0)  =  =  ^  • 


Fig.  13.  Diffraction  efficiency  T|  as  a  function  of  coupling  strength  y^L  for  various  read 
beam  intensities  IAR(0).  The  normalized,  equal  intensities  of  the  write  beams  are  unity: 
^aw(O)  =  ^bw(O)  ~ 


For  the  case  of  <J)  =  0,  which  corresponds  to  a  local  index  grating  that  follows  the 
intensity  pattern,  y  (and  hence  k)  is  a  purely  imaginary  number  that  we  take  to  be  equal 
to  iy(ii^).  The  expression  for  diffraction  efficiency  simplifies  to9 


4Iaw(0)Ibw(^) 


[IAR(0)  +  IAW(0)  -  IBW(0)]‘  +  4Iaw(0)Ibw(0) 


sin  k;L 


(26) 


Apart  from  the  multiplicative  factor  before  the  trignomctric  function  and  the  read  beam 
intensity  dependence  of  K;,  this  result  is  reminiscent  of  Kogclnik’s  expression  (Eq.  (7)) 
for  diffraction  efficiency  of  uniform  gratings.  For  a  weak  reading  beam  [IAR(0)  «  IAW(0), 
IBW(0)]  and  writing  beams  of  comparable  intensity  [IAW(0)  =  IBW(0)],  the  multiplicative 
factor  and  the  input  intensities  dependent  factor  s/(2I0)  that  is  present  in  k;  are  both  equal 
to  unity.  Hence  Eq.  (26)  reduces  to:  r|  =  sin2(yiL/2),  where  y;  =  rcn^Xcose,)  is  the 
amplitude  coupling  constant  in  Eqs.  (21);  this  expression  is  identical  to  Kogclnik’s 
formula  (Eq.  (7))  with  n,  s  2n10.  If  the  writing  beams  have  different  intensities,  then 

Eq.  (26)  simplifies  to 


T|  = 


,4g_ 


a  +  P  r 


sin ’(f> 


(27) 


which  is  identical  to  the  expression  obtained  by  Kukhtarcv  ct  al.n  for  replay  of  a  local 
photorefractive  grating  with  a  weak  read  beam.  For  small  coupling  strengths,  yxL  «  1, 
and  Eq.  (27)  reduces  to  T|  «  py^L^Kl+P)2],  which  is  identical  to  Kogelnik’s 

expression4  for  small  arguments  of  the  trignomctric  function.  Increasing  the  read  beam 
intensity  serves  to  reduce  the  multiplicative  factor  in  Eq.  (26)  until  it  is  zero  for  a  strong 
read  beam  [IAR(0)  »  IAW(0),  IBW(0)L  implying  erasure  of  the  grating  from  readout.  The 
s/2I0  factor  that  is  present  in  is  not  very  sensitive  to  changes  in  read  beam  intensity 
and  decreases  slightly  from  unity  as  the  intensity  of  the  weak  read  beam  is  increased.  Once 
IAR(0)  >  ^aw(O)’  W(0)>  factor  increases  and  is  again  equal  to  unity  for  a  strong  read. 
Note  that  T]  is  symmetric  with  respect  to  change  in  sign  of  Kit  which  implies  that  there  is 
no  nonreciprocity  in  i\  for  local  photoinduccd  gratings,  even  when  the  read  beam  is 
intense. 


For  the  general  case  of  a  nonzero  spatial  phase  shift  between  photorefractive 
grating  and  fringe  pattern  and  a  read  beam  intensity  comparable  to  that  of  the  writing 
beams,  the  parameter  dependence  of  diffraction  efficiency  is  studied  numerically.  Figure  14 
shows  the  nonlinear  behavior  of  rj  as  a  function  of  the  read  beam  intensity  Iar(0)  for 
various  phase  shifts  <J).  Note  that  the  nonlinearity  is  more  pronounced  and  that  higher 
values  of  T|  are  possible  when  §  =  n  /  6.  However,  photorefractive  gratings  formed  at 
<{>  =  7t  /  2  are  more  difficult  to  erase.  The  curve  for  (J)  =  7t/2  is  identical  to  that  shown  in 


Fig.  11.  For  (J)  =  0  case,  both  the  multiplicative  factor  and  s/2I0  in  ka  decrease  as  the 
reading  beam  intensity  is  increased,  and  the  decrease  of  the  first  function  is  more  rapid. 
The  coupling  strength  y/L  is  equal  to  10  in  the  calculations,  so  that  for  an  initial  read 
beam  intensity  of  0.1,  K;L  is  equal  to  2.38.  This  is  considerably  more  than  rr/2  at  which 
a  sine  function  peaks,  so  even  though  the  argument  of  the  sine  function  is  decreasing 
from  the  initial  value  of  2.38  with  increasing  read  beam  intensity,  r\  still  increases  with 
Iar(0)  until  the  factor  s/2I0  in  ^reduces  to  2 jr/tyL)  =  0.63.  Beyond  that  value,  r|  starts 
decreasing  mainly  due  to  the  rapidly  falling  values  of  the  multiplicative  factor  with 
Iar(0).  The  curve  for  §  =  n/6  in  Fig.  14  represents  a  combination  of  effects  present  in 
((>  =  7t/2  and  $  =  0  cases.  Figure  15  shows  the  nonreciprocal  behavior  of  r[  with  respect 
to  direction  of  readout  for  various  phase  shifts  (J)  and  a  read  beam  intensity  equal  to  that  of 
the  writing  beams.  For  equal  intensity  writing  beams,  a  change  in  the  sign  of  the 
coupling  constant  is  equivalent  to  readout  from  the  other  input  port.  Note  that  <J)  =  0 
(<))  =  n/2)  gives  oscillatory  (exponential)  behavior  of  r\  as  a  function  of  ly|L,  while 
<J>  =  7r/6  exhibits  a  combination  of  both  behaviors,  with  decaying  oscillations.  The 
nonreciprocity  is  enhanced  for  <J>  =  rt/6,  while  <J)  =  0  case  is  symmetric  with  respect  to 
readout  from  either  input  port  even  when  the  read  beam  is  as  intense  as  the  write  beams. 
Figure  16  is  a  plot  of  r\  versus  lyiL  for  a  weak  reading  beam  (IAR(0)  =  0.01).  The 
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Fig.  14.  Diffraction  efficiency  ri  as  a  function  of  read  beam  intensity  IAR(0)  for  various 
values  of  spatial  phase  shift  <(>.  The  modulus  of  the  coupling  strength,  lylL,  is  chosen  to 
be  10,  and  the  writing  beams  arc  taken  to  be  equal;  i.c.,  1AW(0)  =  I1)W(0)  =  1. 
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Fig.  15.  Diffraction  efficiency  r\  as  a  function  of  coupling  strength  lyiL  for  various  values 
of  spatial  phase  shift  <J).  The  read  beam  is  taken  to  be  equal  to  the  writing  beams;  i.e., 
Iar(0)  =  Iaw(0)  =  W0)=1. 
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Fig.  16.  Same  as  in  Fig.  15,  but  for  IAR(0)  =  0.01. 


diffraction  efficiency  is  now  a  symmetric  function  of  coupling  constant  for  all  values  of 
<j),  and  the  oscillatory  behavior  for  §  =  0  reaches  a  peak  value  of  unity,  as  expected  from 
Eq.  (26) — for  the  choice  of  parameters  in  this  figure,  the  slant  of  the  gratings  is  zero, 
leading  to  a  Bragg  matched  readout  with  100%  maximum  diffraction  efficiency. 


III.  SUMMARY 


Kogelnik’s  formula  for  diffraction  efficiency  of  a  volume  phase  grating  is  derived 
for  Bragg  matched  readout.  The  differences  between  a  photorcfractive  index  grating  and  a 
thick  phase  grating  of  uniform  amplitude  and  phase  arc  outlined.  For  a  weak  read  beam 
and  polarization  asymmetry  during  writing  and  readout,  the  diffraction  efficiency  of  a 
photorefractivc  grating  is  independent  of  the  read  beam  intensity;  it  is  a  function  of 
grating  thickness,  coupling  constants  during  writing  and  readout,  and  input  intensity  ratio 
of  writing  beams.  For  a  strong  read  beam,  the  diffraction  efficiency  is  a  nonlinear  function 
of  the  read  beam  intensity.  The  diffraction  efficiency  is  also  nonreciprocal  with  respect  to 
readout  from  the  two  input  ports.  The  nonlinear  and  nonreciprocal  properties  of  diffraction 
efficiency  are  enhanced  when  the  phase  shift  has  a  value  lying  between  0  and  n/2. 
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